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ABSTRACT 


AN ANALYSIS OF CARRIER PHASE JITTER IN AN MPSK RECEIVER 

UTILIZING MAP ESTIMATION 
BY 

BRIAN THOMAS KOPP, B.S., M.S. 


Doctor of Philosophy, Engineering 
Specialization in Electrical Engineering 
New Mexico State University 
Las Cruces, New Mexico, 1994 
Dr. William P. Osborne, Chair 


The use of 8 and 16 PSK TCM to support satellite 
communications in an effort to achieve more bandwidth efficiency 
in a power-limited channel has been proposed. This project 
addresses the problem of carrier phase jitter in an M-PSK receiver 
utilizing the high SNR approximation to the maximum aposteriori 
estimation of carrier phase. In particular, numerical solutions to 
the 8 and 16 PSK self-noise and phase detector gain in the carrier 
tracking loop are presented. The effect of changing SNR on the 
loop noise bandwidth is also discussed. These data are then used 
to compute variance of phase error as a function of SNR. 
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Simulation and hardware data are used to verify these 
calculations. The results show that there is a threshold in the 
variance of phase error versus SNR curves that is a strong 
function of SNR and a weak function of loop bandwidth. The NI- 
PS K variance thresholds occur at SNR's in the range of practical 
interest for the use of 8 and 16-PSK TCM. This suggests that phase 
error variance is an important consideration in the design of these 
systems. 
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Chapter 1 

INTRODUCTION 

"Observe due measure, for right timing is in all things 
the most important factor. " 

Hesiod, 

poet & philosopher, 700 B.C. 


The information age is moving forward at an impressive 
speed. In recent years, global telecommunications traffic has 
maintained a growth rate of no less than 15 % compounded 
annually [ 1 ], This, coupled with the increased demands for 
scientific data from space and atmospheric based remote sensing 
and experimentation, has created a strong interest in increasing 
the capacity and performance of available communication satellite 
resources. Toward this end, New Mexico State University's space 
telemetering and telecommunications research center has been 
pursuing research projects whose thrust is the investigation of 
more efficient methods of utilzing an existing satellite channel. 

One area of particular interest has been the study of spectrally 
efficient modulation formats. 

A current modulation standard for satellite communications 
is phase shift keying. Specifically, binary phase shift keying 
(BPSK) and quadrature phase shift keying (QPSK) are used 
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extensively. For example, Intelsat earth station standards A, B, C, 

E, and F have provisions for the use of QPSK in single channel per 
carrier (SCPC) modes and time division multiple access (TDMA) 
burst communication modes [2]. Further, the National Aeronautics 
and Space Administration's (NASA) most extensive satellite 
communication system, the tracking and data relay satellite 
(TDRS) system, uses BPSK and QPSK for all of its user services and 
TDRS tracking, telemetry, and command links with the exception 
of one space shuttle service [3]. 

By modulating with BPSK it is possible to convey one bit of 
information in a given channel through the transmission of one of 
two phases of a sinusoidal carrier. With QPSK, four phases are 
used and thus two bits can be conveyed in that same channel. It is 
possible to extend this process through the use of eight phases to 
send three bits (8PSK) or, going even further, using sixteen phases 
to convey four bits (16PSK). Both extensions can utilize the same 
channel as the BPSK and QPSK signals. However, if the signal 
power level is maintained there will be a corresponding increase 
in bit error rate (BER) as the distance between adjacent 
constellation points drops with the conveyance of more phases. In 
satellite systems, where both available channel bandwidth and 
available power are usually at a premium it is desirable to 
consider the use of 8PSK and 16PSK. For example, if it is desired to 
increase the performance of a QPSK link without changing the 
signalling speed, 2/3 rate 8PSK trellis coded modulation (TCM) 
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could possibly be employed. TCM utilizes specific constellation 
signal set expansions and convolutional encoding to achieve a 
decrease in BER over uncoded QPSK. If a two-fold increase in data 
throughput is needed and loss in performance can be tolerated it 
may be possible to consider 16PSK. The economical key to making 
these extensions in an existing satellite channel is that the channel 
must be a relay or "bent pipe" channel. If the channel involves a 
regenerative transponder onboard the spacecraft this obviously 
won't work without modification to the satellite. It is much 
cheaper if the extensions involve engineering changes in the 
ground communications equipment only. 

With that in mind much research is being conducted in the 
area of 8PSK and 16PSK [4] - [14]. Specifically the 8PSK and 
16PSK TCM schemes have been paid particular attention in the 
literature [4] - [12]. The current available research examines 
higher order PSK systems, and specifically TCM systems, all 
chiefly in terms of bit error rate performance. One important area 
in the consideration of the use of these coherent schemes that has 
not been fully addressed in the literature is the synchronization 
performance of the receiver. Specifically, the quantity of carrier 
phase variance, or jitter, in the reference generated in the 
receiver for demodulation has not been calculated for popular 
receiver designs, at various signal to noise ratios (SNR), under 
tracking or "locked" conditions. This characteristic is an extremely 
important one since jitter in the carrier translates into uncertainty 


3 


of the location of the received PSK constellation points. If this 
jitter is large enough it can cause a degradation in BER 
performance. Further, it will set a minimum range of operational 
SNR below which the jitter will be so severe that phaselock in a 
feedback-type receiver will be unattainable and/or 
unmaintainable. This is of particular importance in 8 PSK and 
16PSK where the constellation points are much closer together 
than BPSK and QPSK. In the literature it is possible to find 
research that does analyze the impact of certain amounts of 
carrier phase jitter on the BER performance of 8PSK and 16PSK 
receivers [5]. However, a quantitative analysis of how much jitter 
is to be expected is generally unavailable. 

Since the quantity of jitter depends on the type of receiver 
used, the first task of this project was to decide which receiver to 
analyze. The PSK signal being received is a double-sideband 
suppressed carrier (DSB-SC) signal. The lack of a discrete carrier 
component coupled with discrete phase transistions conveying the 
data makes coherent demodulation of PSK signals a challenging 
task. There are several receiver designs that were considered for 
use in this study and they can be grouped roughly into two major 
types: forward acting designs and feed-back designs. The decision 
of which receiver design to use was based on several factors. First, 
a practical receiver design was a criterion since it was intended 
from the start to build an M-ary PSK (MPSK) link to verify 
theoretical and simulation results for BPSK, QPSK, 8PSK, and 
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16PSK. Second, the results of this study would be of great value to 
NASA, the principal benefactor of the telemetering and 
telecommunication center, so choosing a type of receiver that they 
use and are familiar with, particularly at the TDRS White Sands 
ground terminal (WSGT) twenty-five miles northeast of NMSU, 
was also a consideration. At the WSGT a type of loop called a 
modified Costas loop is used for the demodulation of both BPSK 
and QPSK. This feedback loop design is based on the loop first 
proposed by Costas [15] for the demodulation of DSB-SC signals. 
The basic modified Costas loop is shown in Figure 1-1 and for 
comparison the original Costas design is included in Figure 1-2. 
The third consideration in the decision process was that the 
choosen design needs to be backwards comparable with published 
data, if possible, to ensure that the research data is accurate. In 
other words, a baseline study with BPSK and QPSK using the same 
design would allow for a comparison with published data on how 
much jitter should occur. When the calculated, simulated, or 
experimental data agrees closely with published data for BPSK 
and QPSK then there is verification that the design being used is 
accurate and extensions to 8PSK and 16PSK can be made. This 
third consideration inferred that a design that is flexible and can 
test BPSK through 16PSK with minimal modification is desirable 
and that the choosen design be popular enough to have already 
received attention in the literature in so far as BPSK and QPSK 
carrier phase jitter are concerned. 
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Figure 1-1. The modified Costas loop. 


DSB-AM signal 



Figure 1-2. The Costas loop. 
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To facilitate making a choice, several receiver designs in 
each of the categories mentioned above were examined. Two 
forward-acting designs were examined. The times-N loop, shown 
in Figure 1-3, is a forward-acting type loop in which the received 
MPSK signal is multiplied by itself N times. The N in the loop's 
name refer's to the M in MPSK. By performing this multiplication 
the data is effectively multiplied out resulting in a spectral 
component at a frequency of N times the carrier that can then be 
tracked by a standard phaselock loop (or even filtered by a 
bandpass filter). Unfortunately, as N, i.e., M increases, the 
frequency at which the carrier phase tracking loop must operate 
increases. For example, a 370 mega-Hertz (MHz) 16PSK signal 
would require a "times 16" loop. The quiescent frequency of such 
a phaselock loop would be approximately 6 Gigahertz (Ghz). 


r(t) = sin((O c t+0:(t)) 



Figure 1-3. The "Times N" loop. 
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Implementing a carrier tracking loop at such a frequency in 
firmware or hardware can be prohibitory because of the 
difficulties involved with working at such a high frequency. 

Another forward-acting design, the feed-forward carrier 
synchronization system, is an approach that contains no feedback 
loop [16]. The system estimates carrier phase by calculating its 
arctangent over 2N+1 samples where N is varied depending on the 
signal to noise ratio. The feed-forward system is shown in Figure 
1-4. This system is used in burst communication applications since 
the possibility of hang-up during acquisition, while not eliminated, 
can be significantly reduced. 1 In evaluating this loop in terms of 
the decision criteria it is noted that there is a significant amount 
of digital hardware and 2 pairs of phase detectors present in the 
loop. While this alone is not a significant deterrent when 
compared to some of the other designs, and in particular the 
design finally chosen, it is noted that this design is substantially 
different from that used by NASA at the WSGT. It contains no 
feedback loop, as does the modified Costas loop, and its unique 
ability to acquire a received signal quickly is not of great 
importance at the WSGT where a particular communication link 
between the ground station and a user spacecraft, called a "sho", 


1 Hang-up is a term for a condition that can occur in PLL’s. When the intial 
reference phase is 180 out from the lock point a phase error of zero will 
occur at the output of most phase detectors. However this point is an 
unstable equilibrium point for the PLL and thus noise will cause the phase 
to move toward the lock point. This process can be slow to start and the 
result is a long acquistion time or "hang-up" before the PLL locks. 
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may last many minutes and result in the transmission of many 
gigabits of data. 


r(t) -yfl sin(co c t+0 j(t)) 



Data 


Figure 1-4. The feed-forward loop. 

The second type of PSK DSB-SC receiver design that was 
considered is the feedback loop design. This loop has several 
variations. Costas loops, demodulation/remodulation loops, 
reverse modulation loops, and maximum aposteriori phase (MAP) 
estimation loops are all forms of feedback loops. All of these 
feedback loops involve the use of some form of PLL to provide a 
coherent reference for demodulation, except in the reverse 
modulation loop where a bandpass filter can be substituted [17, 
pg. 177]. There are several feedback designs that fall into a special 
category called decision-directed loops. A loop can be considered a 
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decision-directed loop if at some point in the receiver the 
transmitted baseband data is effectively recovered and used in 
the demodulation process. Examples of decision-directed loops are 
the QPSK Costas loop with hard limiters and the MPSK high-SNR 
approximation to the MAP estimation loop. 

A feedback loop that is a type of Costas loop and that 
showed early promise as a possible candidate for the research into 
carrier phase jitter is a loop referred to as the Leclert and 
Vandamme variation [18]. It is shown in Figure 1-5. In this loop 
an error signal that drives a VCO is calculated from several hard- 
limiter processors, each related to the recovered data. The hard 
limiters provide the inputs to dual polar quantizers which 
estimate the transmitted phase. These estimates are then used 
along with the hard limiter outputs to calculate an error signal. 

Two more types of feedback loops are the demodulation- 
remodulation loop [17, Fig. 6.6, pg. 176] and the reverse 
modulation loop [17, Fig. 6.7, pg.177]. In the former, the recovered 
data is modulated onto the VCO reference and then phase 
compared with the received signal. If good decisions are being 
made then the only phase difference between the two signals will 
be due to a carrier phase offset. In the reverse modulation loop 
recovered baseband data is "reverse" modulated on to the 
received signal. The resultant signal will have a spectral 
component at the carrier frequency that can be tracked. Further, 
in reverse modulation loops, like feed-forward systems, the 
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Delay 
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Figure 1-5. 8-PSK Leclert and Vandamme variation. 















chance of hang-up during acquistion can be minimized by 
replacing the PLL with a bandpass filter. 

The fourth type of feedback loop that was considered is the 
MAP estimation loop. This loop uses an optimum aposteriori phase 
estimation calculation, based on the recovered data, to determine 
what the most likely transmitted signal carrier phase was. This 
complex estimate is then used to calculate an error signal to drive 
a VCO which provides the coherent reference. 

In the early stages of the receiver design selection process 
the feedback loops received the most consideration. All of the four 
feedback designs mentioned above are relatively adaptable for 
construction in an MPSK system. In contrast, the times-N loop is 
adaptable to MPSK environments for construction but only at a 
considerable cost. To support an M of 2, 4, 8, or 16 would, in 
effect, require that four receivers be constructed if the task were 
completed with analog circuitry. A digital implementation would 
require substantial circuitry as well, and some form of wide 
bandwidth numerically controlled oscillator to achieve the range 
of from twice-the-carrier-frequency to 16 times the carrier 
frequency. This lack of adaptability and the fact that times-N 
loops aren’t employed at the WSGT were sufficient drawbacks to 
warrant withdrawal of this type of loop from consideration. 
Although easily adaptable to the MPSK environment, the feed- 
forward systems distinct variation from the feedback technique 
used at the WSGT coupled with significantly large overall 
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hardware requirements, compared to several of the feedback 
loops, effectively eliminated it from consideration as well. What 
remained at this point in the selection process was to decide 
which feedback loop to use for the research on carrier phase jitter. 

Of the four feedback loops two stand out for their overall 
simplicity of implementation. These are the Leclert and 
Vandamme variation and the MAP estimation loop. The other two, 
the demodulation-remodulation loop and the reverse modulation 
loop, both require increased circuit complexity to implement when 
compared to the first two loops and thus were eliminated from 
consideration. Having narrowed the decision to the 
Leclert/Vandamme and MAP estimation loops, it was decided that 
simulations should be performed to compare the performance of 
the two loops. While the methodology for conducting these 
phaselock loop simulations will be discussed in detail later in this 
dissertation, it should be noted here that the simulation source 
code and results are presented in Appendix A. The results indicate 
that the performance of the Leclert and Vandamme variation was 
less satisfactory than that of the MAP estimation loop. This 
performance difference is attributable to the Leclert and 
Vandamme’ s discrete processing algorithm involving only 
hardlimiters which results in a loop that performs like a "bang- 
bang" servo control. It should be stressed that the difference in 
performance was minimal. It should also be noted that the Leclert 
and Vandamme variation requires a slight increase in digital 
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circuitry compared to the MAP estimation loop. Neither of these 
reasons was significant enough to dismiss the Leclert and 
Vandamme variation by themselves. However, this loop design is 
relatively new and little research has been published on it in the 
literature. In contrast, research into the use of MAP estimation for 
demodulation was first investigated in the early 1950's [19] and 
much has been written about it in the literature since [20]-[23], 

[24, chap. 5], [25]-[28], [29, chap 11], [30]-[36]. This interest in the 
literature for MAP estimation was the final reason that lead to the 
selection of the MAP estimation loop design over the Leclert and 
Vandamme variation. The MAP estimation loop is relatively easy 
to model and construct in an MPSK environment. It is very similar 
to the modified Costas loop employed at the WSGT (in fact at high- 
SNR's they are mathematically equivalent), and there is a large 
enough source of literature research to verify that the models 
developed and tested, and the hardware constructed, were done 
properly. Yuen, at The Jet Propulsion Laboratory, [37, chap 5] and 
Hinedi and Lindsey [35] have both published data related to the 
jitter that proved useful in verifying the results of this study for 
BPSK and QPSK and thus allowed the extensions to 8PSK and 
16PSK to be considered valid. 

One restriction was placed on the MAP estimation loop to 
simplify its modelling and to make the data from the study more 
useful to NASA at the WSGT. The MAP estimation technique uses 
an algorithm that provides an estimate of the phase of the 
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received signal. This algorithm is mathematically complicated and 
involves the use of hyperbolic trigonometric functions. Before the 
advent of digital signal processing, implementing this algorithm 
with analog circuitry would have required a prohibitive amount of 
electronics. However, if it is assumed the received signal exists in 
a high-SNR environment then the MAP estimation algorithm 
reduces to an equation that is easily implemented. Further the 
analysis of the loop is significantly simplified as well. A similar 
reduction in complexity occurs if a low SNR environment is 
assumed. However, it is the high-SNR approximation that was 
particularly attractive to this study. The modified Costas loop used 
at the WSGT is equivalent to the QPSK Costas loop with hard 
limiters and thus equivalent to the high-SNR approximation of the 
MAP estimation loop. Further the above mentioned published 
results on BPSK and QPSK carrier phase jitter are applicable to a 
comparison with the original MAP estimation loop only at high- 
SNR's. If a feedback loop that employs the high-SNR 
approximation to the MAP estimation loop is used, the results are 
completely comparable. Insofar as presenting these results to 
NASA was concerned it is possible to inform them what they 
should expect in terms of carrier phase jitter if they continue to 
use the same type of receiver design for 8PSK and 16PSK. 

Having concluded the selection process and chosen a 
receiver design that utilizes the high-SNR approximation to the 
MAP estimation of carrier phase, it was possible to begin the 
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analysis of carrier phase jitter. The analysis begins with 
understanding jitter and what it is, specifically, in feedback loops. 
To accomplish this one of the most basic of feedback loops, the 
phaselock loop, can be used. The next chapter of this dissertation 
presents, as background, a brief discussion of the basic phaselock 
loop: how it works and how its jitter can be analyzed. Chapter 3 
presents the MAP estimation loop as a phaselock loop for the 
carrier recovery of PSK DSB-SC signals. The approximation for 
high-SNR is described as well. In Chapter 4 the high-SNR MAP 
estimation loop is reduced to a baseband model to simplify 
analysis and simulation. It is here that the equation that describes 
carrier phase jitter in the high-SNR approximation loop is 
presented. Although this equation has been previously solved 
analytically for BPSK and QPSK it is shown that for 8PSK and 
1 6PSK an analytical solution for several of the components of the 
jitter do not exist and another approach is necessary. In the next 
chapter (Chapter 4) numerical solutions are obtained for the 
components for 8PSK and 16PSK. The accuracy of the numerical 
techniques is verifyed using the BPSK and QPSK analytic solutions. 
The 8PSK and 16PSK data could then be tentatively assumed to be 
correct. Through the use of simulation techniques and a laboratory 
hardware design for the high-SNR approximation loop the 
numerical data could be verified more assuredly. All of the results 
for carrier phase jitter are presented together in Chapter 6. 
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As the reader progresses several other important factors 
and characteristics will come to light as well. For example the loop 
noise bandwidth exhibits some interesting phenomena that are 
not previously discussed in the literature and are of great value to 
the MPSK receiver designer. This bandwidth changes significantly 
with SNR, narrowing as the SNR drops. It is attributable to the 
dropping phase detector gain as the SNR decreases. Smaller design 
loop damping factors can minimize this effect. Further, as the SNR 
decreases, a threshold region is reached where the trade-off for 
carrier phase variance that exists between SNR and loop 
bandwidth is no longer linear. In this region, if it is desired to 
maintain a certain level of jitter while dropping the SNR by 3 
decibels (dB), it is not possible to merely narrow the loop 
bandwidth by 3 dB. The loop bandwidth would require 
significantly more reduction to maintain the same level of jitter at 
half the SNR. This nonlinear exchange does not exist at higher 
SNR's where maintaining the same amount of jitter while dropping 
the SNR 3 dB can be achieved by narrowing the loop bandwidth 3 
dB. 

Another interesting, and alarming, characteristic that was 
revealed in the study of the 8PSK and 16PSK carrier tracking 
loops concerns the minimum SNR that will maintain carrier lock in 
the receiver. As will be shown, the carrier phase variance is large 
for 8PSK and 16PSK at low SNR's that would be considered 
operational for BPSK and QPSK. At these SNR's 8PSK and 16PSK 
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show difficulty in maintaining phaselock. While this data was not 
obtained through the mathematical use of acquistion models, it 

was observed in both the simulators and the hardware at similar 
SNR's. 

All of these results will be presented in the subsequent 
chapters where appropriate and summarized with the jitter data 
in the results chapter (Chapter 6). Chapter 7 is the conclusions 
chapter in which results are discussed and suggestions offered for 
what to do next, now that the MSPK jitter has been quantified. 
Following the conclusions, there are several appendices. Appendix 
A presents the Leclert and Vandamme simulation comparison 
with the high-SNR approximation to the MAP estimation loop. 
Appendices B and C conain the optimum MAP estimator 
derivations for 8PSK and 16PSK. The data for two of the 
components of phase error variance are in the next two 
Appendices. Appendix D contains the equivalent noise variance 
data and Appendix E contains the phase detector gain data. The 
variance of the phase error data is tabulated in three Appendices. 
Appendix F has the theoretical results for BPSK and QPSK as well 
as the numerical results for 8PSK and 16PSK. The simulation 
results are in Appendix G and the hardware results are in 
Appendix H. 
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Chapter 2 

THE PHASELOCK LOOP 

"Here and elsewhere we shall not obtain the best 
insight into things until we actually see them growing 
from the beginning ... " 

Aristotle 

philosopher, 384 - 322 B. C. 

2.1 The Basi r Phaselock Loop 

Investigating the MAP estimation loop carrier phase jitter 
requires a model of the loop's steady state operation. The jitter 
that it is desired to measure occurs when the loop has acquired 
and is locked. To build the mathematical model from which the 
jitter can be extracted, and to extract the jitter itself, several 
mathematical procedures must be performed. To understand 
these procedures and their significance it is prudent to discuss 
first their purpose in the most basic of feedback synchronizers: 
the phaselock loop. 

The in-depth analysis of PLLls is available in many sources 
[24] [29] [38] [39]- The discussion presented here is meant only to 
cover those procedures that are required to understand carrier 
phase jitter and how it can be evaluated. The basic PLL is used to 
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phaselock two sinusoidal signals that are 90° apart in phase. The 
quadrature relationship is necessary when a multiplier type phase 
detector (the most common for PLL's) is used. The first of the two 
sinusoidal signals, the received signal, has an unknown phase 
component and is denoted 

r(t) = sin (co c f + d;(t)). (2-1) 

The second signal is the sinusoid generated in the receiver 

v{t) = *j2A v cos((o c t + 0 0 (tj). (2-2) 

The constant amplitude coefficients are chosen for convenience. 
Note that the phase arguments, 0,(t) and G 0 (t), are considered to 

have spectrums significantly below their common carrier 
frequency of co c , i.e., they are varying slowly [39, pg. 21]. For the 

locally generated signal to phaselock to, or track, the received 
signal, information regarding which signal is lagging and which 
signal is leading is required. An error signal which conveys this 
information and provides equal amounts of absolute error, 
whether a lead or lag condition exists, can be achieved with a 
device that yields an odd function of the phase difference 
between the signals. The device that provides the error signal in a 
PLL is called a phase detector. The multiplier-type phase detector 
multiplies the two sinusoids together and lowpass filters the 
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product to remove the double frequency term. The output from 
such a device with inputs (2-1) and (2-2) and a gain of K d (which 

has dimensions of volts~ x ) is the odd function error signal 

<*0 = K d ^F r \ sin(e,(<) -»„(/)). (2-3 ) 

Besides the phase detector, two other basic components are also 
part of the PLL. They are the linear loop filter and the voltage 
controlled oscillator (VCO). All three components are shown in the 
PLL of Figure 2-1. The function of the loop filter varys with the 
application of the PLL but in general its output can be represented 
by convolving its impulse response, /(/), with its input (2-3) as 

x(t) = j f(t - X)K d 4F r \ sinf^A) - e o {X))dl . ( 2-4) 

o 



Figure 2-1. The phaselock loop. 
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The VCO is a unique device in that its output is a sinusoid 
whose frequency varies in proportion to the voltage of its input 
signal. When the VCO's input is zero its output has a constant 
frequency o) c . When the VCO's input isjr(f), the output of the loop 

filter, its output frequency, will be 


c o'(t)=co c + K 0 x(t ). (2-5) 

Note that co'(t) is thus the instanteous frequency of the locally 
generated sinusoidal signal since it is the VCO's output that 
provides the reference sinusoid in Figure 2-1. Since the 
instantaneous frequency and overall phase of the reference 
sinusoid are related by the expression 

4*M s(0 '( t ) ( 2 - 6 ) 

dt v ’ 

and since the overall phase of the reference in (2-2) is 

<i> o {t) = (0 c t + d o {t) (2-7) 

the deviation of the VCO fromtu,., namely K a x{t), can be expressed 
in terms of the phase component of the reference sinusoid, 0 o {t). 
Using (2-5 )-(2-7), this relationship is 

^ = *•„*(<)• ( 2 - 8 ) 

dt 
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Combining (2-3), (2-4), and (2-8) it is possible to state the 
integro-differential equation that describes the operation of the 
PLL as 


“ K ° K * V*VA.J/(< - A)sin(e,(A)- e o (X))dX . (2-9) 

at 0 

An extremely important quantity that has not been 
mathematically introduced yet and that is needed to further 
explore the integro-differential equation of (2-9) is the phase 
error 


e e (t)=e i (t)-0 o (t)^ (2-10) 

This term represents the phase difference between the received 
sinusoid and the locally generated sinusoid. Since it is the purpose 
of the PLL to track the received signal's phase and to do so does 
not necessarily require absolute knowledge of the signal's phase, 
it is equivalent to define the function of the PLL as that of 
maintaining as small a phase error as possible between the 

received and locally generated sinusoids. Introducing the phase 
error, 6 e (t ) into (2-9) yields 

A j/(! - i) sm (e,(X))di . (2-11) 
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Equation (2-9), or (2-11), lead directly to the baseband (low 
frequency) model of Figure 2-2. Note that the VCO is modeled as 
an integrator corresponding to the relationship between the VCO 
input frequency deviation and its output phase change. This 
model is useful for two reasons in particular. It provides direct 
access to the phase error and it has eliminated the need to keep 
track of the frequency, a> c . While the latter reason may seem 

insignificant, if this frequency is included when the loop is 
simulated (as in a complete time simulation of the loop) valuable 
computing resources will be consumed in unnecessarily high 
sampling rates. This model demonstrates that it is not necessary 
to simulate the loop in such a manner. Rather, a baseband 
simulation will suffice. 


Loop Filter: f(t) 



Figure 2-2. The nonlinear baseband model. 
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This nonlinear model is certainly useful for studying 
transient response behavior and acquistion phenomena but its 
nonlinear nature makes it unattractive for studying steady-state 
characteristics where a linear approximation will suffice and thus 
allow the use of some very powerful analysis tools. The carrier 
phase jitter we are interested in analyzing can be obtained 
through the use of a linear approximation to the model in Figure 
2-2. The assumption that must be made to linearize the loop 
model is that the phase error is small. This assumption is 
predicated by the notion that the loop has acquired and is 
tracking the received signal. The loop is in steady state. If the 
phase error is small then 

sin(0 e (f))»0,(O* (2-12) 


This approximation provides the necessary change to convert the 
nonlinear model of Figure 2-2 into the linear model of Figure 2-3. 
Given that the model is now a linear one it is possible to describe 
its transfer function. To do so consider the further modified model 
of Figure 2-4 where now the Laplace transforms of the loop filter 
and the VCO have been introduced. Further, the assumption has 
been made that the Laplace transforms of the phase components 
of the received signal and the reference signal exist. It is now 
possible to construct the transfer function for the loop. While 
there are several transfer functions available in this control loop 
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the desired one for assisting in the calculation of the carrier phase 
jitter is the relationship between the received signal phase, 0 ,( 5 ), 
and the VCO output phase, e o (s), namely 



Figure 2-3. The linear baseband model. 



Figure 2-4. The Laplace transformed linear baseband model. 
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From Figure 2-4 it is clear that 0„(s) can be expressed as 


e c ( s )=e r (s)4p;A v K d F(s)^. 


(2-14) 


Using (2-10) and (2-14) to construct the desired transfer function 
yields 


4P r \K d K,,F( S ) 

«*) l + 4F r KK d K 0 ^- 

s 

To develop the transfer function further requires that the loop 
filter be selected. A commonly used loop filter is one that creates 
a second order control loop [38, pg. 9]. Such a loop is 
unconditionally stable (A third order loop is not). Further, if the 
loop filter is constructed using a perfect integrator the loop can 
track out static frequency differences between the received signal 
and the locally generated signal (a first order loop cannot). To 
achieve a second order control loop that can track out static 
frequency differences requires a loop filter transfer function of 

F(*) = l + - (2-16) 
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where the constant a is the constant of integration for the loop 
filter integrator. Inserting (2-16) into the last part of (2-15) yields 
the transfer function 

H(s)= ^ 

s 2 +sJF r A v K d K 0 + a JF r A v K d K o ' (2 " 17) 


If the phase detector gain, VCO gain, and VCO magnitude 
coefficient are combined to form a new constant K, the transfer 
function becomes 


m- 

s z +s^P r K + aJP r K 


(2-18) 


where K = K d K a J\. This transfer function now has the form of the 
classical second order control loop transfer function 


H(s) 


s2C(Q n +0)l 
s 2 + s2£a) n + co l 


(2-19) 


where £ is the damping factor of the loop and o)„ is the natural 

frequency of the loop. The relationships between these two 
classical loop design parameters and the gain terms, ^l\K and a, 
are 
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fi\K = 2Z(Q n , 


( 2 - 20 ) 


JF r Ka = (0 2 n , 


or conversely 


_> BEE 

2 } 


and 


( 2 - 21 ) 


0} n =^fJp^Ka. 

Note that throughout this analysis the signal power, P r , has 
appeared. It is an important component of the gain of the loop. If 
the signal strength changes or fluctuates during operation, the 
loop performance will change. 

One last characteristic of the PLL that will be introduced is 
its loop noise bandwidth. Practically speaking, the PLL is a 
bandpass filter [17, pg. 170]. Its input is a sinusoid in the presence 
of channel noise and therefore has a certain SNR associated with it 
(assuming some form of predetection filtering). Its output is 
(when tracking properly) a sinusoid of exactly the same frequency 
with significantly increased SNR. The baseband model from which 
the transfer function is derived has a lowpass characteristic, thus 
the name baseband. Its noise bandwidth reflects how wide this 
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lowpass characteristic is by equating an equivalent ideal 
"brickwall" lowpass bandwidth through the use of one of the 
integrals [38, pg.30] [39, pg. 124] 


h = 


1 

2tt|//(0)| 2 


\\H{(ofdco 



J| 


( 2 - 22 ) 


It is important to note that the result of calculating the loop 
bandwidth with either part of (2-22) is a number expressed in 
Hertz. For the second order loop of interest this bandwidth is well 
known and is expressed as [38, pg. 31] 


B, = co r 




4C. 


(2-23) 


using the classical design parameters or as Viterbi presents it [24, 
pg. 36] 


'L ~ ~ 


(2-24) 


when the gain terms are used. It is important to note here that 
equations (2-22) - (2-24) all yield the same result, i.e., the 
double-sided loop noise bandwidth. 
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2.2 The Phaselock Loop With Noise: Calculating litte r 


Having analyzed the PLL and obtained a transfer function 
from its linearized model it is now useful to introduce noise into 
the PLL and examine how the loop models change. Introducing 
noise into the loop will also allow for calculation of the carrier 
phase jitter since it is, afterall, noise that causes jitter to occur. 
Before proceeding it is reiterated that this simple PLL 
presentation is readily available in the literature [24] [29] [38] [39] 
and is a compilation of the information available there. Consider 
now the composite signal 

r(t ) = -j2P r sin (a c t + 0,(r)) + w(t) (2-25) 

where the noise term, w(t), is additive white Gaussian Noise with a 
two-sided power spectral density of watts/Hz. At the front 

end of the PLL there is a predetection bandpass filter, centered at 
^ Hz and with a bandwidth of B IF Hz. At the output of this filter 

the composite signal has become 

r(t) = ^|2P' r sin( (o c t + 0, (t )) + n(t) (2-26) 


where 


n(t ) = V2/ij (r)cos(fi) c f) - V2n 2 (0 s i n ( 6, c f ) 


(2-27) 
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is a narrowband noise process whose power spectrum is shown in 
Figure 5a. The noise components of n(t), namely, ni (t) and n 2 (t), are 
zero-mean, independent, white Gaussian noise processes [29, pg. 
28]. Figure 5b shows the baseband power spectrum associated 
with both «,(?) and n 2 (t). The power in n(t) is calculated 

P n =2B lF ^ = B IF N 0 . (2-28) 


S w (f) 

* 
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Figure 2-5a. The narrowband noise power spectrum. 
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Figure 2-5b. The baseband noise power spectrum. 
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Since the power in the signal is P r the SNR at the input to the PLL 
is 


SNR r (2-29) 

N,B, f 

The composite signal at the input to the PLL can now be expressed 
r(t) = -JlF r sin(G) c r + 0,(0) + V2n, (f)cos(fi> c r) - 'J2n 2 {t)sin(a c t) . (2-30) 

When this signal is applied to the phase detector along with (2-2), 
the VCO signal, the result is the noisey error signal 

e(t) = K d JF r A v sin{6 i (t)-e o (t)) 

+ A' rf A v ni(O cos (-0 o (O) • (2-31) 

~K d A v n 1 { Osin(-0 o (O) 

Noting that the cosine function is even and the sine is odd (2-31) 
becomes 

e(t) = K d ^[P~ r \sin(e i (t)-d 0 (t)) 

+ ^A v n,(Ocos(0 o (O) • ( 2 ‘ 32 ) 

+ K d A v n 2 ( Osin(0 o (O) 

The signal amplitude coefficient, ^ , is sometimes incorporated 
in the noise components of (2-32) so that when the variance of 
the phase error is calculated it can be interpreted in terms of the 
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SNR at the input to the PLL [38, pp. 26-27][39, pp. 106-107]. This 
incorporation is accomplished by the manipulation 


e{t) = K d Jp' r A v 


sin(0, (0 - e 0 (0) + ^=cos(0 o (0) + sin(0 o ( t )) 


(2-33) 


The noise components in the brackets of (2-33) are often referred 
to collectively as the equivalent noise of the loop, stated [39, pg. 
106] 


n e (t) = ^=cos(0 o (O) + ^=sin(0 o (O). 


(2-34) 


The error signal is thus 

e(t) = K d ^A v [sin{6 i (t)-d 0 (t)) + n e (t)]. (2-35) 

The output from the loop filter will now be 

*(0 = J fit- A)A: rfA /P^A v [sin(0,(A) - 0 O (A)) + n e (X)]dX . (2-36) 

o 

Therefore, the integro-differential equation for the loop can be 
expressed as 
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de,(<) <* 3(0 

dt dt 


K c K d .JP r \\ f(t - A)[sin(^(A)) + n e (X)]dA . 
0 


(2-37) 


This equation leads to the nonlinear model of Figure 2-6. 


0i(t) 



Figure 2-6. The nonlinear baseband model that incorporates noise. 

The next step is to linearize the model and develop its 
transfer function. Once this is accomplished the jitter can be 
calculated as the variance of the phase error. To calculate the 
variance requires that the statistics of the equivalent noise be 
stationary. Therefore, before proceeding with the development of 
the linear model the statistics of the equivalent noise must be 
examined. In (2-34) the equivalent noise is seen to be strictly 
dependent on the VCO phase, e o (t). However, it is noted that the 


35 



loop bandwidth is, in general, much narrower than the bandwidth 
of the input noise for carrier tracking loops such as this [24, pg. 
31]. Thus the main lobe of the sine function that is the 
autocorrelation function of either of the noise components, n,(r) 

and n 2 {t), will be narrow, relative to the time r = — In other 

B l 

words the loop, which has a lowpass-type response, responds 

slowly to changes in the input noise and while the VCO phase is 

dependent on noise that occurred on the order of — seconds ago, 

b l 

it is not dependent on noise that occurred at the input in the 
range of — seconds ago. This makes esoteric sense too when it is 

DjF 

considered that the function of the PLL is to filter out the noise at 
its input and thus produce a sinusoid that has less noise, i.e., jitter 
associated with it. As for the equivalent noise, it can now be 
stated that the VCO phase, in effect, does not vary during the 
small time intervals over which the noise is correlated. Then for 

the purpose of evaluating the statistics of the equivalent noise it 
can be written 

n e (t) = ^cos(6 0 ) + ^sin(d 0 ) (2-38) 

where 9 0 is a constant. The autocorrelation function of the 
equivalent noise is 
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(2-39) 


— rtj (/+!)«! (0cos 2 (e o ) 

P r 


E\n e (t+x)n e (t)]=E 


+-^-n 2 (t+x)n 2 {t)sm 2 (d 0 ) 

+— n, (/ + r)n2(Ocos(0 o )sin(0 o ) 

P r 

+— n 2 (^+ T)ni(O sin ( 0 o ) cos ( 0 o) 


Combining like terms and using the linearity property of the 
expected value function results in 

E[n e (t + X )n e (t)] = —{£[«! (t + x) ni (t)}cos 2 {e o )+E{n 2 (t + T)n 2 (0]sin 2 (fl o )} 
*r 

+ — cos(0 o )sin(0 o ){£[nj(f+ x)n 2 (t)] + E[n 2 (t + r)ni(r)]} 

E r 


(2-40) 

However, since n x (t) and n 2 (t) are independent and have zero mean 
their crosscorrelation functions, E[n,(f+ x)n 2 {t)] or E[n 2 {t + x)n x {t)] 
equate to zero as can be seen for the former crosscorrelation term 

E[n { ( t + r)n 2 (r )] = E[n x ( t + x)\E[n 2 ( 0 ] = 0 • (2-41) 

Having equated the crosscorrelation terms of «,(/) and n 2 {t) to 
zero the autocorrelation function of the equivalent noise becomes 
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E[n e {t + T )n e (0] = J- {E[n x (t + l)n x (r)]cos 2 (d 0 ) + £[n 2 ( t + r)n 2 (r)]sin 2 (0 o )} . 


(2-42) 


Since n, (t) and n 2 {t) are identically distributed their 
autocorrelation functions are equivalent and thus (2-42) reduces 
to 


E[n e (t+ x)n e {t)] = ^-E[n x (t+ T)« 1 (O]{cos 2 (0 o ) + sin 2 (0 o )} = -^-£[n 1 (r+ t)«,(0]. 


(2-43) 


Noting that n x {t), like n 2 (t), is stationary, its autocorrelation is thus 
a function only of the time difference, t. Therefore, it is possible 
to write the autocorrelation of the equivalent noise as 


*n e ( T ) = 



(2-44) 


Since the autocorrelation function is the Fourier transform pair of 
the power spectral density of n e (t) it is possible to state 




for |/|<^. 


(2-45) 
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Having identified the equivalent noise as a stationary Gaussian 
process and having obtained its power spectral density, it is 
possible to obtain the variance of the phase error by analyzing the 
output power spectral density of the linear loop model in response 
to the equivalent noise. 

Linearizing the model of Figure 2-6 leads to the model of 
Figure 2-7. If the input phase, 0,(r), is held constant and equal to 

zero then the output of the linear system, i.e., the VCO phase, will 
be due only to the equivalent noise input. The transfer function 
from the equivalent noise input to the VCO phase output is seen in 
Figure 2-7 to be identical to that of (2-18). For a linear system the 
output power spectrum can be obtained using 

•%.(/Htf(/)l\(/)- (2-46) 



Figure 2-7. The linear baseband model that incorporates noise. 
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Noting that with 0, (0- 0 the phase error is the negative of the VCO 

phase and so the output power spectral density of the VCO phase 
is also that of the phase error. This is stated 

So e (f) = Hffs ne (f)- (2-47) 

Since the equivalent noise has zero mean and therefore d e {t) does, 
too, the variance of d e (t) is its autocorrelation function evaluated 
at zero 



T=0 


\Se t (f)e^df 


-oo — OO 


(2-48) 

Next, consider that since the power spectral density of the 
equivalent noise, as described by (2-45), has a constant value 
over the entire frequency range that |//(^)| is essentially non-zero, 
it is possible to express (2-48) as 

B !F 

(2-49) 

r JB !f 
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where (2-45) has been used for the power spectral density of the 
equivalent noise. Using the definition for the loop noise bandwidth 
given in (2-22), and noting that the loop noise bandwidth is much 
smaller than the IF bandwidth, the variance of the phase error 
can be expressed as 

< 2 - 50 > 


Since the loop transfer function of (2-18) has unity DC gain the 
jitter becomes 

.A&A (2-51) 


Using the input SNR to the PLL, described by (2-29), the variance 
of the phase error of the PLL becomes 



1 


B, 


2 SNR r B [F 


(2-52) 


A final, useful though intangible quantity that is often considered 
in understanding the performance of a PLL is the SNR in the loop. 
It is usually described in terms of the variance of the phase error 

as 


SNR l = -^- = 2SNR r ^- 

b l 


(2-53) 
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The quantity is considered intangible because there is no signal 
that exists in the loop that exhibits this SNR. Its function is as a 
reference for measuring performance and it will be encountered 
again in the descriptions of the performance of the high-SNR MAP 
loop. 

As has been demonstrated jitter in a PLL is represented by 
the variance of the phase of the VCO's signal due to noise at the 
PLL input. This jitter is dependent on the SNR of the received 
signal and on both the IF and loop bandwidth. All of these factors 
will appear again, in one form or another, when the jitter in the 
high-SNR MAP approximation loop is analyzed. The procedure will 
be very similar as well. The linear baseband model for the loop 
will be developed and its transfer function used to obtain the 
jitter of the VCO's phase due to a noisy input. 
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Chapter 3 

MAP ESTIMATION DERIVATIONS FOR CARRIER PHASE JITTER 


"It is truth very certain that when it is not in our 
power to determine what is true, we ought to follow 
what is most probable. " 

Rene Descartes 

mathematician & philosopher, 1596 - 1650 


3.1 Introduction 


In a coherent MPSK system we have the transmitted signal 



(3-1) 


of duration T seconds and symbol energy E, where 6 m (t) is the 

phase shift key modulation. For M = 2, 0 m (t) can take on the values 

(0, n). For M >2 the modulation can take on the values -—* 1 ^ for 

M 

m = 0,l...(M — l). The channel adds white Gaussian noise (AWGN) to 
the signal and shifts the carrier phase by an unknown amount, 0,. 

It is considered here that the function of the coherent receiver is 
to detect the data in the received PSK signal in the presence of 
noise. The received signal (after predetection filtering) is denoted 

r (0 = cos {(O c t - G m (t) + 0, ) + n(t). ( 3-2 ) 
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As was done with the PLL of Chapter 2, the noise component, n(t), 
will be assumed to be narrowband noise having passed through 
the predetection bandpass filter 1 before entering the receiver. For 
MPSK, the optimum receiver can be achieved with the use of two 
correlators [40, pps. 257-259]. These correlators require a 
coherent reference. In order to obtain this reference it is 
necessary for the receiver to estimate what the unknown phase 
component, 0,, is. There are several criteria that can be used to 
achieve the optimum estimate of 0, [24, pg. 126]. However the 

estimation technique of interest here is the maximum aposteriori 
criterion. This approach involves selecting 0, in such a manner 
that the aposteriori probabilty of 0, , given that r(t) was received, 

is maximized [26]. The a posteriori probability density of the 
phase 0 ( , given received signal r{t), is denoted p(0,|/?) . 2 It can be 

expressed in terms of the channel transition probability density 
and the received signal and carrier phase probability densities 
using Bayes theorem 


n(0\R\- P ^ di ^ di) 
p[m ~ 


(3-3) 


where p(/?|0,) is the transitional probability density of receiving 
the signal R given a carrier signal with phase 0 ( . We are 


1 Note that the bandwidth of the predetection filter is wide compared to the 
symbol rate. 

2 Note that time designator has been dropped for mathematical 
convenience. 
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interested in maximizing this quantity to obtain the best estimate 
of the carrier phase and will therefore be computing its 
derivative. However, before this is done the natural log function is 
introduced to separate the terms in (3-3). This will not affect the 
derivative since the natural log function is a monotonically 
increasing function. Applying the log function results in 

ln(p(e,|R)) = ln(p(/i|e i )) + ln(p(e,))-ln(p(/i)) (3-4) 

Now (3-4) can be maximized by taking its derivative with respect 
to the carrier phase and setting it equal to zero. We add the 
further constraint that the phase of the optimum receiver's 
coherent reference will equal 0, when the maximum a posteriori 

phase estimate is obtained. This will help formulate the estimator 
into a control loop. 

Noting that the probability of receiving R is not a function of 
the carrier phase (3-4) can be written as 

9i=9 0 

Since the carrier phase 0, is uniformly distributed it is possible to 
state 


01n(p(0,|/?)) 

dO t 


00 , 


01n(p(0,)) 


= 0. (3-5) 
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Be 




|0,=0 O 


de, 


= o. 


0,=e o 


(3-6) 


It is now possible to derive the maximum aposteriori phase (MAP) 
estimator by using p(/?|0,), the likelihood function. 


3.2 MAP Derivation For BPS K 

To proceed with obtaining p(F\e i ) the M=2 case (BPSK) signal 

set will be considered first. Once its MAP estimator is obtained the 

subsequent estimators for M=4,8, and 16 will be presented. To 
compute p(/?| 0, ) begin with the received signal 


r{t) = 



di) + n(t). 


(3-7) 


The MPSK signal can be expressed in two dimensions using the 
orthonormal basis 



which results in r(/) = r, (*)%(*) + r 2 (/)'P 2 (/) where 
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(3-9) 



and 



Performing the integration for r, and r 2 yields 


rj = 4 e cos(0, - 0 m ) + ni - 

r 2 = -4 e sin(0, - 0 m ) + n 2 = s 2 + n 2 . 


(3-10) 


To continue, the likelihood function can be obtained from the 
conditional density of R given both the modulation angle 6 m and 
the carrier phase 0, . Noting that this conditional density can be 
expressed in terms of the orthonormal basis components r { and r 2 
it is expressed 

P (F\e i ,d m ) = p^di^pfald^dn) (3-11) 

The conditonal densities on the right of the equality in (3-11) are 
the noise probability densities of and n 2 expressed in the 
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received signal space with means ^ and s 2 respectively. The 
conditional density of (3-11) can therefore be expressed 3 





( r i- j i ) 2 

2 a 2 


V5 


KG 


( r 2 ~ s 2 ) 2 
2 a 2 


(3-12) 


To obtain the likelihood function the modulation angle random 
variable is integrated out of the conditional density in (3-12) as 
follows 


p( R \ e i)= lp{ R KO m )p(O m )dd m . (3-13) 


For BPSK the modulation angle probability density is 

p( e m) = j${O m )+^S(d m -x). (3-14) 

Substituting (3-14) into (3-13) and applying the sifting property 
of the delta function the likelihood function becomes 

p( R \ e i) = \p( R \QnQ m =Q) + ^p(R\Q i ,O m = n). (3-15) 


3 Note that the noise samples out of the quardrature 
Gaussian and have variance <7 2 = 

2 E 


arm correlators will be 


48 



To insert (3-12) into (3-15) the quantities r l -s l and r 2 -s 2 must 
first be analyzed for both values of 6 m . For 0 m =0 

rj - s, = r, - sfE CO s(6>, - 0 ) = r x -^E 008(61,) ( 3- 1 6a) 


and 


r 2 ~ * 2 = r 2 + sin(0, - 0) = r 2 + ~JE sin(0,). 


(3-16b) 


For d m = K 


r l -s ] =r l -^cos(d i -7t) = r l +y[Ecos(e i ) (3-17a) 

and 

r 2 -s 2 =r 2 + *J~E sin(0, -n) = r 2 - 4 e sin(0 i ) - (3-17b) 


Now using (3-16) and (3-17), it is possible to insert (3-12) into (3- 
15) and obtain the likelihood function 
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( r \ — 4 E cos( 0, )) 2 ( r 2 +4E sin ( 0/ )) 2 


p m=^\ 


2<J Z 


2 a z 


(n +V£ cos(0; )) 2 (r 2 “V£ sin(0j )) 2 

e 2<r 2 ^ 2cr 2 


( 3 - 18 ) 


Completing the square and combining exponential terms yields 




(V 2 -2r\ V£cos(0 ( - )+£cos 2 (0j)+r 2 + 2 r 2 V£sin(0; )+£’sin 2 (0/ )j 

2cr 2 

” __ 2 _ 


e 

+- 


|r 2 + 2 rj V£cos( 0 f )+£cos 2 ( 0 , )+r| -2r2 V£sin( 0 ; )+£sin 2 ( 0 , )j 
— • — — - ■ — > 

2 


( 3 - 19 ) 


By reducing, the conditional density becomes 
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pW«.) = ^ 


[ ( r ? +r 2 + ^) (-2rj V£cos(0,)+2r2 V£ r sin(0 ! )) 

j 2a 2 


2tt<7 


2 a x 


( r ? +r 2 +£ ') ( 2 ri-y/£cos( 0 ;)- 2 r 2 V£sin(0,)) j 

2cx 2 ,, 2CT 2 


( 3 - 20 ) 


Rewriting the exponents in the form of the hyperbolic cosine 
produces 




(rf+rl+E) 


2 o 1 


(2rj 4 E cos( 6j )- 2 r 2 4e sin (0; )) 


2<t z 


(2rj V£ cos( $ 1 )— 2/*2 V£ sin(0, )) 


2 cr 


( 3 - 21 ) 


Inserting the hyperbolic cosine function results in 


p(R\8i) = 



(r 2 + r 2 +£ ) 

2 <t 2 


( 

cosh 

V 


(lr\ V E cos( dj ) - 2 r 2 V# sin( 0, )) 
2 c 5 


A 


( 3 - 22 ) 
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Substituing in r, and r 2 from ( 3 - 9 ) the conditional density 
becomes 


p(*| 0 ,) = 


K+i+f) 


'1 2 
* " 2 c 2 


27rcr 


cosh 


r ry 

V£ cos(0,)J r ( f )-]j— cos {(Q c t)dt 


VjB sin( 0 ( )j r(f)J— sin(ft) c r)cfr 
n V 7 " 


( 3 - 23 ) 


Combining the integrals yields 


p(R|e,) = 


(rf+r|t£) 


2;tc7 


2 ( 7 ^ 


x cosh 


2 £ 


- jr(d(cos(d) c f)cos(0 l )-sin(© c f)sin(0 1 ))<ir 


( 3 - 24 ) 


Applying a trigonometric identity results in 
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1 


f 


p(K|e,)= 


M 


2 +rj+E ) 


2a 


2 KG 


cosh 


v 


2 E 


~ 2 ~ j r(t)cos(o) c t + 6 t )dt 


(3-25) 


To continue with deriving the MAP estimator in (3-6) the log 
likelihood function, 


ln[p(*|0 f )] = ln 


27TCT 2 


(rf+rl+E) 

la 1 


cosh 


2 E 


- J r(t)cos[(O c t + 0, )dt 


(3-26) 

is used. The MAP estimator is obtained by taking the derivative of 
(3-26) with respect to 0, . To compute this derivative we note the 

multiplicative property of the logarithm function and the 
following application of the chain rule 

<9{ln[cosh(«)]} = 1 0[cosh(«)] _ sinh(«) d(u) = f > d(u) 

d6 l cosh(«) 00, cosh(w) 00, 00, 


(3-27) 


Taking the derivative of (3-26) results in 
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[ r \ +r 2 +E ) 

<?ln 

^ r 9vt2 

2 kg 1 

ddi 


- tanhl 


2 E 


-JV(r)cos(Gy + 0 f )<# 



x "^T" J r (0 sin {(O c t + 0, )dt 


(3-28) 

Next, note that the quantity rf + rf, the square of the magnitude of 
the received signal, is not a function of0, and therefore the first 

term of the sum in (3-28) equates to zero. To complete the MAP 
estimator derivation set the right side of (3-28) equal to zero and 
insert the necessary condition for the optimum estimator 
mentioned in (3-6) regarding the control of the receiver reference 
(i.e., e i = 0J. The resulting MAP estimator can be used in a control 

loop where a loop filter and voltage controlled oscillator can be 
used to drive the MAP estimator equation, or more appropriately 

the error signal, to zero [26]. To conclude, the BPSK MAP estimator 
equation is [27] 


2 E 


- tanh 


-\r(t)cos((D c t + 0 o )dt 


x]r{t)sin(G) c t + e o )dt = 0. (3-29) 
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When it is employed in a control loop to drive an error signal to 
zero the result is the loop of Figure 3-1. 



Figure 3-1. The optimum MAP estimator loop for BPSK carrier 

tracking. 


Historically, the implementation of complicated functions 
like the hyperbolic tangent function was a costly procedure in 
terms of the large quantity of analog hardware it required. As a 
result alternative approximations were considered. In the case of 
the hyperbolic tangent function two approximations have been 
considered for implementation in the loop of Figure 3-1. The first 
is accurate when the function argument is large. In that case the 
hyperbolic tangent function may be replaced by a hard limiter. 
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This is known as the high-SNR approximation. When the function 
argument is small the hyperbolic tangent function can be replaced 
by a "through connection" since the function has a linear unity 
gain response to small arguments. The loop of interest in this 
study is the high-SNR approximation to the MAP estimation 
technique and it is shown in Figure 3-2 for BPSK. 



Figure 3-2. The high SNR approximation loop for BPSK carrier 

tracking. 


3.3 MAP Derivation For OPSK 

To obtain the optimum MAP estimator, and subsequently 
the high-SNR approximation for QPSK the derivation presented for 
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BPSK must be altered starting with the probability density 
function for the modulation angle,p(0 OT ). Once this is done the 

likelihood function for QPSK can be used to arrive at the 
appropriate estimator. Therefore to begin, the QPSK modulation 
angle probability density function is 




(3-30) 


Once again using the sifting property to integrate out the 
modulation from the conditional density of R given both the 
modulation angle d m and the carrier phase 0, results in the 

likelihood function 



(3-31) 


To proceed (3-12), the conditional noise probability density, is 
inserted into (3-31) for each of the four values of 0 m . To do this 

the two quantities r x - s, and r 2 - s 2 must first be analyzed for the 
four values of 0 m . For 0 = — 

m 4 

r x -s x =r x -4E cos^ 0,-^1 (3-32a) 
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and 


r 2 - s 2 = r 2 + V£ sinf 0, - — . 

V 4) 


For e m 


3k 

T 


r \ ~ - 4e cos^0 ( - ^ 


and 


r 1 -s 1 = r 2 +4Es\n\6 i ~y 


For e m 


K 

4 


r, -j, =r, 




and 


r 2 ~ s 2 ~ r 2 + sin^ 0 ( - + j . 


(3-32b) 


(3-33a) 


(3-33b) 


(3-34a) 


(3-34b) 


For e m 


3 7T 

T 


58 



n- J l = r l 




(3-35a) 


and 


r 2 -s 2 = r 2 




(3-35b) 


Now using (3-32) - (3-35) it is possible to substitute (3-12) into 
(3-31) and doing so yields 


p(«l»/)=5^r 


2 2 

1 -4e cos^ Si ( r 2 +V£ sin^ j 

„ 2a 2 


2 a 


-4E cos [Ot j \j 2 +JE sin^i 


2 a 


2 a 


,-Vfcos^/—]] ^+VEsin^i~J 


2 a L 


2a L 


VEcos^O,— j^+VEsin^fy — 


2cr 


2 <T 


(3-36) 
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To proceed with minimum confusion, each of the terms in (3-36) 
corresponding to a particular modulation angle will be considered 
separately. Thus, 

|q-VIcos^e,+-^jj |r 2 +VEsin^0,+-^jj 

(A) = — je 2 ° 2 e (3-37) 


corresponds to the first term in (3-36). Completing the square of 
the exponential arguments in (3-37) and then combining them 
results in 

r^-2r\ -Je cos^0, +—^ j+£cos2^0j+-^ 

+r| +2r 2 V £ sin^ dj +— ^j+£sin2 +— ^ 

W-;=r* • < 3 - 38 > 


By combining and separating like terms (3-38) becomes 


;+£-2ri Vfcos 0 1 +-^j+2r 2 V£sin^0,-+-^ 


(A) = 


la 1 


&7t<7 


(3-39) 


Applying trigonometric identities to expand the sinusoidal terms 
forms 
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By combining the integrals and distributing the carrier 
components within the single integrand, (3-42) becomes 



Next, using the trigonometric identities 
cos(*)cos(}>) - sin(*)sin(>>) = cos(x + y) 


(3-43) 


and 


(3-44) 


cos(x)sin(;y) + sin(*)cos(;y) = sin(jc + y) 

the integral in (3-43) is separated into two new integrals. This 
new form is 


<A) = 


rf+rj+E 


8 KG 1 


E ‘ , , . , [e t 

— jr{t)cos(m c t+S i )dt Jr(r)sin(cu (: f+0 / )dt 

0 » / n 


2 


2 c L 


(3-45) 
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At this point the second term of (3-36), the term corresponding to 
a modulation angle offl m =~, is considered. From (3-36) this 

second term is 

^r|-V£cos^0,+^jj ^2+V£sin^0,+^j 

<B> = -i, « ^ e ^ . (3-46) 

' 1 O rrsr 2 ' 


Using the same approach as was done for the first term, (3-46) is 
first resolved into 


(B) = 


r 2 +r 2 +E 

la 2 


8 no 


-2r\ 


cos(0/)cos^j-sin(0 / )sin^j 


xe 


2a z 


2 r 2 V£^sin(^ )cos^^ j+cos(0; )sin^^ j 


xe 


2 a 


(3-47) 


V2 , . V2 


Making use of the fact that cos^^ j and 

second term becomes 


sin 


=4± the 
V 4 J 2 


rf+r^+E -rp/2E[-cos(^)-- sin(^)] r 2 V^[-sin(e,)4-cos(ei)] 


(B) = 


1 9 cA 

- e la e 


2 o L 


2 a 


8/T(7 


(3-48) 
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Inserting the equations for r, and r 2 , given in (3-9) converts (3- 
48) into 


(B) = 


8n<j 


rf+rl+E 

2 a 1 


l r ( t )-\l~cos(co c t)c/t 


V2£[-cos(0, )-sin(6>, )] 


x e 


2 . 


LO 


MOy-sin(t» c /)<fe 


2 a 

V2£[-sin(0j )+cos(0, )] 


Xe 


2a l 


(3-49) 


Once again combining the integrals and distributing the carrier 
component inside the integrand yields 


(B) = 


r? +ri +E 

1 "^7^ 

T e 10 

8 kct 


— I r(f )[-cos(fl> c r)cos(0, )-cos(<u £ ./)sin(0, )+sin(fl> c r)sin(0 i )-sin(a> c r)cos(0, )]c/i 


xe 


2a A 


(3-50) 


Using the same trigonometric identities as with the first term and 
separating the integral in the same manner yields 


(B) = 


1 

8 no 2 



[e t 


" 

r e t 1 

r~ 4 -r\ +£ 

- 2 J-MO cos (^+0/V/ 
’ 1 0 

M. 


|-M')sin(<o c r+0,)<* 


2 a 1 e 2 a 1 2 a 1 
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(3-51) 


The third term of (3-36) is 

(C) = — ~~T^ 2? e 2a 2 . (3-52) 

N 8 no 2 

Following the procedure established for the first two terms 


(0 = 


r \ +r 2 


%7tG 


-2 r x 4E 


cos(0/ )cos|^~ j-sin(0j )sin^— ^ j 


xe 


2<t z 


(3-53) 


2r2 4 e 


sin(0 { )cos^— ^ j+cos(0| )sin^-~ j 


x e 


2 <7 l 


Here it is noted that cosf-™) = ^ and that sinf ^ 

Therefore (3-53) becomes 




2 2 
r \ +r 2 


(c) = 


8ttct 


2 a- 


+E -r\ V2£[cos(0 i )+sin(^ )] r 2 V2£ r [sin(0 f )-cos(0,- )] 
e 2<j2 e 


2 a 


(3-54) 
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Once again inserting the equations for^ and r 2 , given in (3-9) 
resolves (3-54) into 


<C> = 


8tt<t 


2 2 

r 1 + r 2 

~ 2a 2 


lr(t)^J-cos(to c t)dt 


V2£[cos(0 { )+sin(0 i )j 


x e 


xe 


2o z 


T [2 

MO^-sin {(O c t)dt 
.0 v 1 


V2£ r [sin(0j)-cos(0 / )] 


2a 


(3-55) 


Once again combining the integrals and distributing the carrier 
component inside the integrand yields 


<0 = 


l 

T e 

8/TCT 2 


rf+rj+E 

2cr 2 


- jr(O[cos(w c r)cos(0, )+cos(« c f)sin(0, )-sin(a) c /)sin(e, )+sin(to c f)cos(0, )\di 


Xe 


2a 1 


(3-56) 

Using the same trigonometric identities as with the first two term 
and separating the integral in the same manner yields 
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<C>= 


rf+rj+E 

1 2^ 

-e 2<J e 


2^jr(t)cos(a> c t+dj)dr 2^ Jr(Osin(© c /+0, )dt 


la 


2a 


8 no 


The final term in (3-36) 


(3-57) 


(D) = 


. 2a 2 


%KO 


2a L 


(3-58) 


resolves into 


(D) = 


8 no 


r\+r\+E 

2a 2 


x e 


xe 


-2r\ %/£ 

cos(0/)cos 

(-?: 

j-sin(0/)sin 

(-t) 

] 

2r 2 4E 

sin(0/ )cosj^ 

la 1 

-t) 

fcos(0;)sin| 

:-t)] 

(3-59) 


2a* 


Noting that cos( ~\ = = “> ( 3 " 59 > becomes 


r 2 +rf+£ — T| V2£[-cos(8,)+sin(6j)] r 2 V2£[-sin(ej)-cos(e,)] 


(D) = 


2cr 


2<r z 


2<J 


87r<x 


(3-60) 
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Inserting the equations for r x and r 2 , given in (3-9) resolves (3- 
60) into 


(D) = 


&7T<7 


r \ +r 2 
2a 1 


J r (t\l~cos(co c t)dt 


V f 2£ r [-cos(0 / )+sin(^)] 


Xe 


2a L 


T 

I 

.0 


jr(0jysm(a> c t)dt 


V2£[-sin(0, )-cos(0; )] 


xe 


2 a L 


(3-61) 


Once again by combining the integrals and distributing the carrier 
components within the integrand (3-61) becomes 


(D) = 


1 

87tcr 2 


r 2 +r 2 +E 
2 a 2 


xe 



T 

j r ( / )[ - cos( g> c / ) cos( 6 1 )+cos(co c t)s\n(Gi )+sin(fi> c r)sin(0/ )+sin(<y c r)cos(0 r - )]dt 
0 


2 a 2 


(3-62) 

As with the previous terms, trigonometric identities are used to 
separate the integral in (3-62) into two new integrals. This new 
form is 
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f{x,y) = e x e y +e x e y + e x e y + e~ x e y 


( 3 - 65 ) 


where 


X = 




r{t)cos((O c t + d^dt 


and 


(3-66) 



This sum can be expressed as 

f{ x *y) = [f x +e x }[e y + £ _:y ) = 4cosh(*)cosh(;y). (3-67) 


Applying this new form of the sum in (3-64) yields 


p(fl|0,) = 


1 


r \ +r 2 

2a 1 


2 no 

xcosh 

x cosh 


1 f E t 

~ 2 -W ~ MO cos(a) c t + 0, )dt 

0 v / o 

1 f e t 

— 2 -Jt ^ r (0 sin («c r + 9i)dt 

O Vi 0 


(3-68) 


Next the log likelihood function is obtained and can be stated 
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ln[p(fl|0 { )] = ln 


iKO 2 


r\+r\+E 
2 a 1 


+ ln< cosh 


+ \n < cosh 


-r M ! r(t)sin(a c t + d^dt 
<y V 1 n 


(3-69) 


As with BPSK, the next step is to take the derivative of (3-69) 
with respect to 0,. Once again it is noted that the first term in the 

log likelihood function will equate to zero in the derivative. The 
derivative of the QPSK log likelihood function is stated 




T 

J r(f )sin(fi) c f+ Oi)dt 
o 


T 

x J r(t) cos (co c t + 6 ( )dt 
o 


- tanhl 


^ r(t)cos[(O c t + 0,*)A 
a o 


r 

x Jr(f)sin(<tf c f + O^dt 
o 


(3-70) 

To complete the MAP estimator derivation for QPSK set the right 
side of (3-70) equal to zero and insert the necessary condition for 
the optimum estimator mentioned in (3-6) regarding the control 
of the receiver reference, i.e., d i = e o , just as was done for BPSK. 

The resulting MAP estimator is [32] 
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tanh 


2 E 


-\r(t)sm((o c t + d 0 )dt 


x J r(r)cos(© c f + 6 0 )dt 
o 


-tanh 


2 E 


■\r{t)cos{(O c t + B 0 )dt 


• (3-71) 


x J r(/)sin(© c t + 6 a )dt = 0 


When a loop filter and VCO are used in conjuction with the above 
control signal the result is the carrier tracking loop of Figure 3-3. 
To obtain the high-SNR approximation for QPSK the use of 
hardlimiters is once again employed in replacing the hyperbolic 
tangent functions. The resulting loop is shown in Figure 3-4. 

3.4 MAP Derivation For 8PSK 


The rigorous derivation for the MAP estimator for 8PSK is 
presented in Appendix B. The optimum MAP estimator for 8PSK is 
stated here as 

Jtstanhfayj-ciytanhfax) | 5 2 xtanh(5 2 )')-c 2 >>tanh(c 2 x) (3-7?) 

1 + cosh(c 2 *)cosh(s 2 y) ^ cosh(c 1 x)cosh(5 1 _y) 
cost^q-^cosh^y) cosh(c 2 Jt)cosh(.S 2 .y) 


where 
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Figure 3-3. The optimum MAP estimator loop for QPSK carrier 

tracking. 
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Hard Limiter 


Figure 3-4. The high SNR approximation loop for QPSK carrier 

tracking. 

The operation of this estimator uses the scaled outputs from 
quadrature arm correlators, as is done for BPSK and QPSK, i.e., (3- 
74) and (3-75). Once again the obvious obstacles to 
implementation that (3-72) possesses motivates a simplified 
solution, namely the high-SNR approximation. To obtain this 
approximation hardlimiters are used in place of the hyperbolic 
tangents for large arguments. Further, instead of hyperbolic 
cosines, the following approximation is used [31] 

cosh(A) = |A| » 0. (3-76) 
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This follows since the SNR is incorporated in the calculation of x 
and y as a scale factor. The first step in proceeding with the 
evaluation of the high-SNR approximation is to examine the 
impact of using the hyperbolic cosine approximation of (3-76). 
Inserting (3-76) into (3-72) reduces the estimator equation to 

5,jctanh(5 1 >')-c 1 >'tanh(c ] x) _ s 2 *tanh(s 2 ;y)-C 2 ;ytanh(c 2 x) _ n 

Jc 2 x\ e \s 2 y\ + i g M g M " U - 

l+ e \°iA e M 1 + e M e M 

To consider the impact of using hardlimiters in place of the 
hyperbolic tangents it is noted that the inputs to the hardlimiters 
are the scaled correlator output coordinates (scaled by the SNR) 
describing where the received 8PSK symbol was detected in the 2 
dimensional constellation space. If the received symbol is detected 
in the first quadrant both hardlimiters will output a +1 and the 
estimator equation will then reduce to 


s\x-c\y 

e \ciA e \ s iy\ 

h 


six-c 2 y 

1 + 


= 0, x >0, y>0. 


(3-78) 


Correspondingly, if the received symbol is in the second quadrant 
the hardlimiter corresponding to the inphase correlator will 
output a -1 while the quadrature correlator will output a +1. The 
resulting estimator is 
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(3-79) 


£l£jt£lZ | J 2* + e 2? 

I e h x \ e M 

e \ c i x \ e M 1 + e h x \ e hy\ 


= 0, x < 0, y > 0. 


It follows that for the third quadrant where x < 0 and y < 0 the 
estimator will become 


-Si* + c^y -s 2 x + c 2 y 

e \ c 2*\ e hy\ e M e M 

1 + 1 + 


= 0, j: < 0, < 0 


(3-80) 


and in the fourth quadrant where * > 0 and j < 0 the estimator is 


Z£ifz£iZ | ~ s 2 x ~ c 2y 
e h*\ e foy\ J c \ x \ e M 

1 + 1 + 


= 0, x > 0, j > 0. 


(3-81) 


The four estimator equations of (3-78) to (3-81) correspond to 
using the high-SNR approximation for each of the four two- 
dimensional quadrants. To proceed further each 8PSK sector must 
be addressed individually in terms of the received 8PSK symbols. 
To accomplish this consider first the received 8PSK sector wherein 
* > 0, >> > 0, and * > y. Since this sector is in the first quadrant, the 
reduced estimator equation of (3-78) is used. It can be reduced 
further to 


\ + e { c 2 x+s 2y- c \ x ~ s \y) 


+ 


s_2X - c 2 y 

l + e i c \ x+s \y- c 2 x -s2 y) 


= 0, jc > 0, y > 0. 


(3-82) 
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To reduce the denominators of (3-82) consider again that this is a 
high-SNR approximation and that the outputs from the correlators 
will be approximately the actual transmitted data scaled by a 
large factor that incorporates the SNR. For reception in the first 
sector this is stated 


JC 



=m 


(3-83) 


and 


y ~ fj. sin 





(3-84) 


Inserting (3-83) and (3-84) into the denominators of (3-82) 
results in 


s\x-c x y 

\ + e ^( c 2 c \ +* 2 *i -qq -qq ) 


+ 


s 2 x-c 2 y 

j + e ^(qq +qq - c 2 c \~ s 2 s \ ) 


= 0, x > 0, y > 0, x > y . 


(3-85) 


This reduces to 


1+ ^(2c 2 c,-1) 


s 2 x-c 2 y 
1 + g Ml-2c 2 q) 


= 0, x > 0, y > 0, jc > y. 


(3-86) 
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Since 0 < 2 c 2 q < 1 and is a large number the denominator of the 
first term of the sum in (3-86) will reduce to unity and the 
denominator of the second term will become very large. The 
result will be that the estimator equation of (3-86) will reduce to 

5 i x - c \y = 0, * > 0, y > 0, * > _y. (3-87) 

Cross multiplying off the large scale factor in x and y and noting 
that in the first sector s^=Q and c, = / are the data estimates of 
what were transmitted, the estimator equation can be considered 
as the error signal in the loop that is to be driven to zero and 
expressed as 


e(t) = IQ-QI,I>0,Q>0,I>Q (3-88) 

where I and Q_are the outputs from actual correlators in the 
constructed carrier tracking loop and in the case of no noise at all 
represent vector coordinates to a point on the unit circle. Note that 
(3-88) has been obtained only for the case that the received 
symbol is detected in the first sector, i.e.,jt>0, >>o, and *>>>. In 
fact, this error signal is the same for any received 8PSK symbol 
designated by quadrature correlator outputs I and Q_and utilizing 
received sector locations to identify the data estimates. To show 
this is true the case for the second sector, namely when * > o , >» > 0, 
and y > x , is now considered. 
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Starting from (3-82) and using the the approximations 


X 


=> fl cos 



= [ic 2 


(3-89) 


and 


y ~ H sml — | = fis 2 


(3-90) 


in the denominator forms the estimator 


S±x - Cjy + Six-Cjy 

\ + e ^i c 2 c 2 +S 2 S 2 ~ c \ c 2 ~ s l s 2 ) i ^ e ^( c l c 2 +s l s 2- c 2 c 2~ s 2 s 2) 


= 0, x > 0, y > 0, y > x . 


(3-90) 


This reduces to 


jjx-qy ( s 2 x-c 2 y 

1 + ^( ,_2c i c 2) i + e ^( 2c i c 2-0 


= 0, x > 0, y > 0, y > x. 


(3-91) 


Once again since 0 < 2c 2 q < l and / j . is a large number the estimator 
is reducable, in this case, to 


s 2 x — c 2 y = 0, x > 0, y > 0, y > x . 


(3-92) 
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A loop error signal is obtained from this estimator, as before, by 
cross multiplying off the large scale factor in x and y and noting 
that in the second sector s 2 = Q and c 2 = / are the data estimates of 
what was transmitted. This error signal is then 


e(t) = IQ-QI,I>0,Q>0,Q>I. 


(3-99) 


where I and Q_are defined as for the first sector. Indeed it is 
possible to show that in all of the sectors the same error signal 
will result when the high-SNR approximations used above are 
assumed. This will now be demonstrated in an abridged fashion 
for the remaining six sectors rotating in a counter clockwise 
fashion as was started for the first two sectors. 

In the third sector where the estimator of (3-79) is used 
and, x < 0 , y > 0, and |x| < y , it is assumed that 


X ~ //cos 



-jl cos 


37T^ 

T > 


= -^2 


(3-100) 


and 


y ~ H sin 


' 5n 



= 2 - 


The estimator of (3-79) then reduces to 


(3-101) 
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(3-102) 


£l£+fiy | S2X + C21 
x+e Li{\-2c x c 2 ) 1 + /(2qc 2 -l) 


= 0, x < 0, y > 0, |x| < y 


and finally to 


S 2 X + C 2 y — 0, x < 0, y > 0, x < y • 


(3-103) 


In this sector s 2 = Q and -c 2 = / and so the error signal becomes 


e(t) = IQ-QI,I<0,Q>0,\I\<Q. 


(3-104) 


In the fourth sector where the estimator of (3-79) is again 
used and, jc < 0 , y > 0, and |x| > y > it is assumed that 


X = M cos (y') = ~ fl cos (f ] - -/iCl 


(3-105) 


and 


y»Aisin^ = ^sin^|) = Ais 1 . 


(3-106) 


The estimator of (3-79) then reduces to 


+ = 0, X < 0, y > 0, 1*1 > 0 

1 + ^(2c 2 c,- 1) 1 + e M(l-2c 2 q) 


(3-107) 


and finally to 
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ijX + c x y = 0, x < 0, y > 0, |je| > y . 


( 3 - 108 ) 


In this sector s x = Q and -c x = / and so the error signal becomes 


e(t) = IQ-QI,I<0,Q>0,\I\>Q. 


(3-109) 


In the fifth sector where the estimator of (3-80) is used and, 
x < 0 , y < 0, and * < y , it is assumed that 


X~fl cos 


9k 

— = -a cos 
8 1 


(D- 


-H c \ 


(3-110) 


and 


y » fi sin 


'9k\ 

Cs) 



(3-111) 


The estimator of (3-80) then reduces to 


-ftjr + qy | -S2X + C2I 
l + e M 2c 2 c l-0 i + e K l ~ 2c 2C\) 


= 0, x < 0, y < 0, x < y 


(3-112) 


and finally to 


— S]* + c x y = 0, jc < 0, y < 0, x < y . 


(3-113) 


In this sector -s x = Q and -c x = / and so the error signal becomes 
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e{t) = IQ - QI, I < 0, Q < 0 , / < Q . 


( 3 - 114 ) 


In the sixth sector where the estimator of (3-80) is again 
used and, x < 0, y < 0, and x > y, it is assumed that 


X 


(Un 

“/*COs(-g- 



= -^ 2 


(3-115) 


and 


y~fi sin 


8 , 


"^TtJ 


= -fJ-s 2 . 


(3-116) 


The estimator of (3-80) then reduces to 


-sjx + qy | -y 2 -y + c 2 y 

l + ^( 1_2c 2 c l) 1 + e M(2c 2 q-l) 


= 0, x < 0, y < 0, X > y 


(3-117) 


and finally to 

- 5 2 x + c 2 >' = 0,x<0,y<0,x>.v. (3-118) 

In this sector -s 2 = Q and -c 2 = / and so the error signal becomes 
e(t) = IQ- Q1,I < 0, Q < 0 , 1 > Q. 
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(3-119 



In the seventh sector where the estimator of (3-81) is used 
and, je > 0, y < 0, and x<\y\, it is assumed that 


f 13tt 

x~fi cosl — — | = jl cos 


(t> 


^2 


(3-120) 


and 


. ( 13 ^^ . (- 2)71 

y = ^sin — \ = fisin 


V 8 


8 


= -/** 2- 


(3-121) 


The estimator of (3-81) then reduces to 


-s^x-c^y -S2X-£il 

l + e ^( 1_2c i c 2) i + ^( 2c i c 2->) 


= 0, x > 0, y < 0, x < |y| 


(3-122) 


and finally to 


-s 2 x - c 2 y = 0, x > 0, y < 0, x < |y| . (3-123) 

In this sector -s 2 = Q and c 2 = i and so the error signal becomes 

e{t) = IQ-QI,I>0,Q<0,I<\Q\. (3-124) 

In the eighth sector where the estimator of (3-81) is again 
used and, *>0, y<0, and *>|y|, it is assumed that 
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(3-125) 


x ~ n cos 




=m 


and 


y~H sin 



= -^i* 


(3-126) 


The estimator of (3-81) then reduces to 


-s,* -c } y -J 2 x-C 2 y 
1 + e M( 2 q^2-!) l + e M(l-2qc 2 ) 


= 0, JC > 0, y < 0, x > |y| 


(3-127) 


and finally to 

-s^x -c x y = 0, JC > 0, y < 0, x > \y\. (3-128) 

In this sector -s, = Q and c, = / and so the error signal becomes 

e{t ) = IQ - QI t I > 0, Q < 0, 1 > |G| . (3-129) 

Since all of the error signals for all eight sectors are identical the 
limitations on / and Qare dropped and the error signal becomes 
simply 


e(t) = IQ-QI. 


(3-130) 
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This is the error signal that is required for implementing an 8PSK 
carrier tracking loop that utilizes the high-SNR approximation to 
MAP estimation of carrier phase. 

3.5 MAP Derivation For 16PSK 

The rigorous derivation for the MAP estimator for 16PSK is 
presented in Appendix C. It is stated here as 


—C\y tanh(cjj:) + Sjjc tantivy) 

+ 

+ 

+ 

, ^ cosh(c 2 Jt)cosh(s 2 ;y) ^ cosh^j^cosh^y) ^ cosh(.s 2 ;t)cosh(c 2 y) 

cosh(c 1 j:)cosh(j' 1 ^) cosh(c 1 Jc)cosh(5 1 y) cosh(c ] x)cosh(5 1 >') 
-c 2 ytanh(c 2 Jc) + s 2 x tanh(s 2 y) 

cosh(c 1 x)cosh(5 1 y) | ^ cosh(j 1 Jc)cosh(cj>') ^ cosh(.s 2 ;c)cosh(c 2 ;y) 

cosh(c 2 jr)cosh( 5 2 _y) cosh(c 2 ;c)cosh(j 2 >’) cosh(c 2 x)cosh(j , 2 >') 

-5 1 >'tanh(^ 1 x) + c 1 .xtanh(c 1 ;y) 

cosh(c 1 .x)cosh(j 1 ;y) ^ cosh(c 2 Jc)cosh(.s 2 y) ^ ^ cosh(j 2 j:)cosh(c 2 y) 

cosh(j 1 x)cosh(c ] >’) cosh(^x)cosh(cjy) cosh(5 ] jc)cosh(c 1 >') 

-s 2 y tanh(^ 2 x) + c 2 xtanh(c 2 y) 

cosh(c 1 x)cosh(5 ] >') ^ cosh(c 2 x)cosh(s 2 _y) ^ cosh(5 ] x)cosh(c 1 y) ^ 

cosh(s 2 ;t)cosh(c 2 y) cosh(,s 2 Jt)cosh(c 2 y) cosh(.s 2 Jt)cosh(c 2 ;y) 


(3-131) 


This estimator equation has the same form as that of (3-72) for 
8PSK. In fact the procedure for analyzing the high-SNR case is 
exactly the same. Each of the 16 sectors is considered individually 
and the error signals turn out to be identical. The result is that the 
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error signs! that is required for implementing a 16PSK carrier 
tracking loop that utilizes the high-SNR approximation to MAP 
estimation of carrier phase is the same as that of 8PSK. In fact the 
error signal is exactly the same for any M-ary PSK signal set. In 
the case of QPSK a review of the high-SNR carrier tracking loop of 
Figure 3-4 shows the error signal of (3-130) is exactly what is 
present at the input to the loop filter. For BPSK consider that the 
high-SNR loop of Figure 3-2 has the error signal 


Since the estimate of what was transmitted on the quadrature 
channel will always equal zero, i.e., 0 = 0, for BPSK the error signal 
of (3-132) may be generalized to that of (3-130). So it is therefore 
possible to conclude that for any M-ary PSK scheme the error 
signal in a carrier tracking loop utilizing the high-SNR 
approximation to MAP estimation has the error signal of equation 
(3-130). This general MPSK high-SNR loop is shown in Figure 3-5. 
Note that the decision device that generates the data estimates is 
based on determining which sector the received 2 dimensional 
symbol was detected in, just as the mathematical discussion of the 
high-SNR approximation used. 
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Figure 3-5. The high SNR approximation loop for MPSK carrier 

tracking. 
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Chapter 4 

DERIVING THE VARIANCE OF THE PHASE ERROR FOR THE HIGH- 
SNR MSPK CARRIER TRACKING LOOP 

"Mathematics is an experimental science and 
definitions do not come first, but later on. " 

Oliver Heaviside 

physicist & electrician, 1850 - 1925 
4.1 Introduction 

In the first chapter the need for quantifying the variance of 
the phase error in the MPSK carrier tracking loop that uses the 
high-SNR approximation to MAP estimation was discussed. In 
Chapter 2 the mathematical derivation involved in calculating this 
variance, or jitter, for a basic PLL was presented. Further, before 
discussing the variance of the more complicated MPSK carrier 
tracking loop it was necessary to discuss the loop itself. This was 
the purpose of Chapter 3. Now having established a mathematical 
foundation of understanding for both the basic idea of jitter and 
the tracking loop of interest, it is possible to proceed with the 
subject of this chapter: obtaining an equation that describes the 
actual variance of the phase error in a high-SNR MPSK carrier 
tracking loop. The procedure is straight forward and follows that 
of Chapter 2. Starting with the integro-differential equation of the 
loop, a linear baseband model is constructed from which an 
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appropriate transfer function is obtained. This transfer function is 
then used to relate the noise in the loop to the jitter or carrier 

phase variance. 

4 ? The Variance Of The Phase Error 

To begin, consider the high-SNR MPSK carrier tracking loop 
presented in Figure 3-5 and shown here again in Figure 4-1 with 
one modification in particular. A gain term on the quadrature arm 
integrators has been introduced. This term allows for the 
normalization of the average energy associated with the I and Q. 
components. It reflects the physical assumption that the receiver 
is utilizing perfect automatic gain control of the received signal. 
When the error signal in the loop is discussed shortly, the purpose 

of the gain term will become more clear. 

The output of the loop filter in Figure 4-1, stated in terms of 

its error signal input, e(t), is 

x(t) = J f{t - A)e(A)dA . 
o 

Using (2-8), which describes the VCO output phase in terms of its 
input, we can form the expression 

<4 ' 2) 

dt n 
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Figure 4-1. The high SNR approximation loop for MPSK carrier 
tracking, modified for error signal analysis. 

Again describing the difference between the received signal 
carrier phase and the VCO phase as the phase error changes (4-2) 
into 


dejj) _ dejt) 

dt dt 


K 0 \f{t-X)e{X)dX. 

o 


(4-3) 


The integro-differential equation cannot be described further 
without exploring the error signal. So far the error signal has been 
presented as a function of the quadrature arm correlator outputs 
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and the data estimates. It must now be described in terms of the 
phase error. 

First, note that the error signal, derived in Chapter 3, carries 
w.th it a time delay of one symbol period associated with the 
quadrature arm integrators. Strictly speaking this should be 

incorporated in the error signal [29, pg. 65] and doing so for this 
loop results in 


e(t) = e p 7 | IQ-Ql] 


(4-4) 


where pis the Heaviside operator and r is the symbol period over 
which the integration takes place. Insofar as the loop bandwidth 
of this type of loop is much less than! the delay component 

is essentially unity over the frequencies of interest and can be 
disregarded [29, pg. 65], It is of interest to note here that while 
this delay does introduce questions as to the performance of the 
loop in the nonlinear sense, (acquisition and mean-time between 
slip characteristics) it does not impact the steady state 
performance analysis of phase jitter when it is neglected. 

To relate the quadrature arm outputs, /and Q, to the phase 
error consider first the output from each of the quadrature mixers 


1 f 2E 

~2y~7 r C0S (^e - 9 m ) + n(t)cos( 6 ) c t + 0 0 ) 
1 l~2 E 

+ 2 V ~T~ cos ( 26 V + @ e — 6 m ) 


(4-5) 
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and 




2 N 
1 


H[-sin(0 e - e„)]+ "(')sinK' + «») 
T 1 

Hsin(2G) c f+^- 0 m) 


( 4 - 6 ) 


a e-e is the phase error and all the phase time 

where e e = 0* - 1S> p mitnut of 

and y(,) 

w m be effectively filtered out. If the gain term is 


G = 


1 


_ 2 _ 

TE 


( 4 - 7 ) 


then the output from the integrators will be 


/ = cos(0^ -6m) + 


( 4 - 8 ) 


and 


Q = -sin(d e -O m ) + N Q 


( 4 - 9 ) 


a m are zero mean Gaussian random 
where the noise terms N, and N fi are 

■ Wes Because the narrowband noise n(r) is necessar y 
variables, because , , Ht h — it can be 

wideband compared with the IF data bandw,dth, 
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considered as a white noise process in so far as the integrators are 
concerned. Doing so allows the following calculation of the noise 
variance at the output of either integrator 




TT 

j f G 2 w(t)cos((o c t + 6 0 )w(u)cos(a> c u+ 9 0 )dtdu . 
oo 


(4-10) 


Bringing the expected value function inside the double integral to 
the noise processes results in 

TT 

°N, =G 2 \\E[w{t)w{u)\cos{(O c t + 0 0 )cos{(o c u + 0 0 )dtdu. (4-11) 

00 


Noting that for a white Gaussian process 
£[w(r)w(w)] = Rywit'U) = -S(t- u) 


(4-12) 


and then applying the sifting property to the outside integral 
converts (4-11) into 


°n, =G 2 ^r I cos ( ft> c r + 9 0 ) cos( a c t + e o )dt . 


(4-13) 


Multiplying out the squared cosine term yields 


° 2 n, ~ G 2 -^-j- + -cos(2co c t + 29 0 )dt 


(4-14) 
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and finally 


2 _ r 2^1 = JJ^I = *<>. (4-15) 

° N ' 2 2 ET 2 2 2 E 

Thus the variance of the noise components out of the quadrature 
integrators is the reciprocal of twice the symbol energy-to-noise 
ratio in the channel. This calculation can be demonstrated 
spectrally by considering the response of the integrator to the 
narrowband noise in the frequency domain. The integrator’s 
impulse response (truely, each is an integrate and dump function 
with a zero-order-hold designed to latch the dumped value over 
the next T second symbol interval) is a flat pulse of duration T 
seconds. It follows that the magnitude of its transfer function is a 

sine function, centered at zero, and with first nulls at the symbol 

rate — Hz. The modulated narrowband noise, which, once again, is 
’ T 

wide compared to the symbol rate, is flat over the significant 
frequencies of interest in the integrator transfer function. It is 
therefore possible to consider the noise output of the integrator as 
being due to its transfer function’s effect on a white noise 
spectrum. This allows for accepting the above calculations 
resulting in the noise variance of (4-15). 

Having identified the components of (4-8) and (4-9) they 
can now be considered when substituted into the error signal. The 
error signal is now stated as 
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e(t) = cos(0 e - 6 m )Q + N,Q - (-sin (6 e - O m ))/ - N Q I . 


(4-16) 


Writing the error signal as a control signal with an additive 
disturbance identifies this loop's equivalent noise term. Doing so 
yields 


e(t) = cos(9 e -6 m )Q + sin(0 e -0 m )l + N e (t ) 


(4-17) 


where 


N e (t) = NfQ- N q 1 (4-18) 

is the equivalent noise in this loop. That is, the IF noise appears 
inside the loop, in the error signal as noise weighted by the data 
decisions. The error signal of (4-17) can be reduced further by 
considering that for this loop, wherein the loop bandwidth is much 
smaller than the symbol rate, the phase error changes very slowly 
compared to a single symbol interval. Since it is essentially 
constant over several of these intervals it is appropriate to 
consider the average of the control function in analyzing the error 
signal. Therefore stating the error signal as its expected value (the 
expected value of the equivalent noise is zero) plus a disturbance 
results in 
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e(t) = E[e(t)] + N e (t). 


( 4 - 19 ) 


The expected value of the error signal is the effective phase 
detector characteristic of the loop [39, pg 116] and is denoted 

PD(e e ) = E[e(t)] (4-20) 

Therefore the error signal itself can be stated as 

e(t) = PD(0 e ) + N e {t). (4-21) 

Recall from Chapter 2 that the output from a multiplier type 
phase detector in the absense of noise is essentially the sine of the 
phase error. The error signal of (4-17) can be shown to reduce to 
exactly the sine of the phase error in the noiseless case if it is 
assumed that good data decisions are being made. Consider the 
noiseless error signal 


e'(t ) = cos(d e -8 m )Q + sin(0 e -0 m )l. (4-22) 

The data estimates can be rewritten as / = cos(0 m ) and <2 = sin(0 m ) 
where e m is the estimate of the transmitted angle. With 
trigonometric manipulation (4-22) becomes 
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(4-23) 


e'(t) = cos(0 e )[cos( 9 m ) sin( 0 m ) - sin ( 6 m ) cos( G m )] 

+ sin(0 e )[sin(0 m )sin(0 m ) + cos(0 m )cos(0 m )] 

Using the I and Q. designators, (4-23) becomes 

e\t) = cos(e e )[l m Q m -Q m I m ] + sin(0 e )[Q m Q m +I m I m \. (4-24) 

Incorporating the assumption that good data decisions are being 
made reduces (4-24) to the expected phase detector characteristic, 
namely 

e'(t) = sin(0 e ). (4-25) 

Returning now to the phase detector characteristic of (4-20), 
and again noting that the expected value of the additive 
equivalent noise is zero, it is possible to describe the phase 
detector characteristic as 

PD(6 e ) = E[e{t)\ = E[cos(0, ~ 0 m )G + sin(0 e -0 M )/]. (4-26) 

Again using trigonometric identities it is possible to manipulate 
(4-26) and reduce it to 


/>D(0 e ) = E[sin(0 ( 


'e + 



(4-27) 
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This last form of the phase detector characteristic is exactly the 
form that will be needed in the next Chapter when the gain of the 
phase detector is calculated. 

This analysis of the error signal was begun to obtain a form 
that included the phase error directly. The form of interest at this 
point is (4-21). Inserting it into (4-7) yields 

d 6 e( ll = - K 0 \ f{t - X){pr{e e (tj\ + N e (t)}dX (4-28) 

dt dt q 

where the time designator for the phase variables has been 
reinstated for uniformity. As in Chapter 2 this form of the integro- 
differential equation that describes the loops operation leads 
directly to a baseband model, in this case shown in Figure 4-2. To 
proceed with linearizing the baseband model of the loop it is first 
necessary to examine the phase detector characteristic more 
closely. It is convenient to express the phase detector 
characteristic in the form [39, pg. 269] 

PD{9 e ) = oc SNR sin(0 e ) (4 29) 

where oc^nr Is the effective gain of the phase detector. It 
incorporates the impact of making poor data decisions in such a 
manner as to become smaller as more incorrect decisions are 
made, i.e., as the SNR goes down. Note that when the SNR is high 
and good data decisions are being made this gain term equates to 
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unity resulting in a phase detector characteristic that is once again 
the sine of the phase error. This gain term manifests itself as the 
slope of the phase detector characteristic measured at zero phase 
error. Mathematically it is expressed as [39, pg. 269] 


a SNR ~ 


_ dPD{6 e ) 


d6„ 


e e =o 


(4-30) 


N-(t) 


e o (0 


TVeM 

PD(0 e (t)) ■ 

/TV e(t) 

J 

xU 


vco 


Loop Filter: f(t) 

'I 

KqJxCiOcIu 

0 


Jf(t-X)e(A,)dX 

0 



' X(t) 



Figure 4-2. The nonlinear baseband model. 

To linearize the loop it is assumed that the phase error is 
small and therefore (4-29) reduces to 

PD{6 e )~a SNR 6 e . (4-31) 

The linear baseband model of Figure 4-3 results. Note also that 
the employed loop filter is the same as before so that the loop will 
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be a second order control loop. From the linear model it is possible 
to develop the second order transfer function, just as was done in 
Chapter 2. For this loop the combined DC gain is a SNR K 0 and the 

transfer function that describes the relationship between the 
received signal phase and the VCO output phase is 

s (XsnrK 0 {s + a) (4-32) 

s 2 + &snrK 0 ( s + a ) 


Ne(t) 


Sift) 

J 

i 

^“SNR 

A 

e(t) 

U>" 

VCO 

Loop Filter: f(t) 


e c (t) 


K 0 Jx(u)du 

0 


jf(t-A,)e(A,)dA, 

M 


x(t) 

0 



Figure 4-3. The linear baseband model. 


The transfer function of equation (4-32) can be equated to the 
classical second order control loop transfer function and its design 

parameters through 


H(s) 


2 W + 

s 1 +2£ '(o' n s + (o' n 2 


(4-33) 
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where 


*°n ~ ^ a SNR 'JKqO- - V &SNR a n 


(4-34) 


and 


C' - ^ a SNR tJ ° ~ ^j a SNR^‘ 


(4-35) 


To proceed with identifying the carrier phase jitter the transfer 
function that describes the relationship between the equivalent 
noise input and the VCO phase output is needed. From Figure 4-3 
this relationship is 


H W = ^(s + a) 

" U s 2 +0WU* + a) 


— 

a SNR 


(4-36) 


To obtain the variance of the phase error the output power 
spectral density of the linear loop model is analysed in response 
to the equivalent noise. Doing so is once again preceeded by the 
notion that the equivalent noise is stationary. In fact, it is 
cyclostationary [35] with autocorrelation function 


RN e (*) = °N e 



M<7. 


(4-37) 
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To see that this is true consider that the equivalent noise itself can 

be represented as a pulse amplitude modulated waveform with 
period T. Further, its amplitude has zero mean and variance . 

Recalling the approach used in Chapter 2 to describe the 
variance of the phase error, the output power spectral density for 
the phase error of the loop will be 

s 9 ,(/)= K (/)| 2 **,(/)• (4 ' 38> 

And it follows that the variance can be described as 


a}, = (4_39) 


The power spectral density of the equivalent noise is obtained 
from its autocorrelation function (4-37) and is stated 

S Ne (/) = On, t sinc2 ( / T ) • (4-40) 

In the case of the PLL of Chapter two the fact that the IF 
bandwidth was large compared to the loop bandwidth allowed for 
a simplification in the calculation of the variance, namely that the 
power spectral density of the equivalent noise is flat over the 
range of lowpass frequencies for which the square of the 
magnitude of the loop transfer function is essentially non zero. 
The same holds true here. The first null of the sine-squared 
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function of (4-40) is at a much higher frequency than the highest 
significant frequency in the square of the magnitude of the loop 

transfer function. Therefore the variance of (4-39) can be reduced 
to 




// , 

T \\H n {ffdf 

-r 


(4-41) 


where /' is an arbitrary frequency slightly greater than the most 
significant frequency to be considered in the square of the 
magnitude of the loop transfer function. Using the transfer 
function relationship of (4-36), the variance becomes 



< t r 


a 


SNR -f' 


Jl 


(4-42) 


Using the same reasoning as in Chapter 2 the integral in (4-42) is 
equated to the loop bandwidth, b l , and the variance is therefore 



<?nJB l 

a SNR 


(4-43) 


And finally, since the symbol period, T, is the reciprocal of the 
symbol rate, SR, the variance of the phase error in an MPSK loop 
utilizing the high-SNR approximation to MAP estimation may be 
stated 


104 


(4-44) 



Bl 

SR a SNR 


Making any comparisons between this equation and that of the 
variance for the basic PLL, namely, (2-52) will be withheld until 
after the components of (4-44) are further explored in the next 

Chapter. 

With the exception of the symbol rate, all of the components 
that make up the variance equation of (4-44) vary with SNR. 

What is more, how they vary with SNR is dependent on which M- 
ary PSK scheme is being used. The first of these, the loop 
bandwidth, is directly dependent on the phase detector gain. As a 
result the loop bandwidth can become substantially narrower as 
the SNR decreases. The phase detector gain itself decreases with 
decreasing SNR. The amount it decreases at a particular SNR is 
greater as M goes up. The third component, the equivalent noise 
variance, increases with decreasing SNR. However, at a particular 
SNR , as M is increased, the variance of the equivalent noise can 
decrease (only at small SNR’s). It never increases with M. 

All of these changing components lead to a dynamic and 
complex relationship between SNR and the phase error variance. 
To proceed further with understanding the variance requires an 
indepth analysis of the above mentioned three components. This 
is the subject of the next chapter. 
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Chapter 5 

THE COMPONENTS OF PHASE ERROR VARIANCE 


"To talk sense is to talk quantities. " 

Alfred North Whitehead 
mathematician & philosopher, 1861 - 1947 

5.1 Introduction 

In Chapter 4 the variance of the phase error in the MPSK 
carrier tracking loop that uses the high-SNR approximation to 
MAP estimation was derived. It is restated here for convenience 
as 



b l 

SR a SNR 


(5-1) 


The three components of this variance that change with SNR were 
also introduced. They are the loop bandwidth, b l , the phase 
detector gain, a SNR , and the variance of the equivalent noise, ol . 

It is the subject of this Chapter to discuss these three components 
in detail. As was mentioned at the end of Chapter 4, each of them 
has a dynamic response to changing SNR and to selected 
modulation type. Their combined effect on the variance of the 
phase error will be presented as results in Chapter 6. 
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5 , 7 The Equivale n t Noise Varian ce 


The first component of the phase error variance that will be 
discussed is the variance of the equivalent noise. This component 
consists of two characteristics as will be demonstrated. The first of 
these reflects the presence of channel noise in the IF receiver 
front-end. The second relates the impact of making bad data 
decisions on the noise in the loop itself. This second characteristic 
is particularly intriguing since it manifests itself in such a way as 
to slightly lessen the overall variance of the equivalent noise and 
therefore the variance of the phase error. To begin the analysis of 
the equivalent noise variance recall the equation for the 
equivalent noise given in (4-18) restated here as 

; (5-2) 

N'(t) = N,Q-N Q I. 

In calculating the variance of the equivalent noise it is first noted 
that the mean of the equivalent noise is 

e[n.(0]=£[n,gMV] = o ’ (5 ' 3) 

since the statistics of N,Q and of N a I are identical. The equivalent 
noise variance is therefore 

a\ = E^b) 2 ] = fiK' 2 + w ?e 2 - ™,N e io\- (5-4) 
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To proceed further, each modulation type is considered in turn. 
For BPSK, (5-4) can be reduced because the squared estimate on 
the I channel is always unity and the Q, channel estimate is always 
zero. For BPSK the variance becomes 

<= £ K1 =<t "c (5_5) 

which is simply the variance of the noise out of the quadrature 
arm integrator. Since noise out of each of the integrators has 
equivalent statistics, equation (5-5) can be equated to equation 
(4-15) and the result is 



For BPSK the variance of the equivalent noise in the loop is exactly 
the variance of the noise out of the quadrature arm integrators. 
For QPSK the variance of (5-4) is expanded as 

o\ = E[N e (t) 2 ] = E[Nli 2 \ + E[NjQ 2 ]-2E[N,N Q iQ\. (5-7) 

In the case of QPSK the data estimates on either channel are 
independent of their orthogonal noise components. In other 
words, / is independent of n q and Q is independent of N,. This 

permits the manipulation of (5-7) wherein 
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< = e[w«(<) 2 ] = £[/ 2 ]£[^] + £[e 2 ]£[jV, 2 ]-2£[w,/]£[Ar s e]. 


(5-8) 

Since the statitics of N r l and n q Q are equivalent (5-8) becomes 
<? N , = £[w e (<) 2 ] = £[/ 2 ]£[^]+£[g 2 ]£[w, 2 ]-2{£(w,7]} 2 . (5-9) 

Note once again that the quadrature arm integrator output noise 
statistics are equivalent and the result is 

CT 2 Ni =£[w,( I ) 2 ] = £[w, 2 ]{£[/ 2 ] + £(fi 2 ]}-2{£[w,/]} 2 . (5-10) 

This becomes 

<y\ = £[w,(<) 2 ] = £jw, 2 ]£[i 2 +S 2 ]- < 5 - n) 

Using (4-15) and noting that the vector magnitude of the data 
estimates is always unity, the QPSK equivalent noise variance is 

< =eK« 2 ] = ^-2H'V‘]} 2 - (5_12) 
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What remains is to calculate fifty/]. This is accomplished by 

equating the expected value of a function of two random variables 
through 

oo oo 

£ [^]= J ! N M N l)f(L)dN,dI m (5-13) 


where the two probability denisty functions are that of the noise, 


/("/) = 




(5-14) 


and that of the transmitted I channel data, 





1 * 
_l_ A 

r, + ^ 2 ) 

r 2j 

T (J 

2 

l" 2 J 


(5-15) 


To examine the integration of equation (5-13) it is noted that/ 
can take on two values depending on which of the two possible l m 
values were transmitted and depending on whether or not ty 
caused an error to occur. The result is that there are four integrals 
corresponding to four combinations of possible transmitted data 
and possible estimated data that must be examined. This is 
expressed by integrating out the transmitted data from (5-13) 
resulting in 
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V 




J_V2 

2 2 


+ ' 


1 -V2 

2 2 


j Nif(Ni)dNj 

V2 

2 

fV2 


/• + • 


1 

2~ 2~~ 


J2 

fN t f(N,)dNi 


(5-16) 




1 V2 

2 2 


jiV;/(iV/)^/ 
V2 
l 2 


The first and second integrals correspond to a transmission of 

/ = ^2l while the third and fourth integrals correspond to sending 
m 2 

i in the second integral the estimate has been made 

m 2 ‘ 

incorrectly because the noise is sufficiently large in the negative 
direction. In the fourth integral the estimate has been made 
incorrectly because the noise is sufficiently large in the positive 
direction. By inserting (5-14) into (5-16) and performing the 
integration the expected value of the product of the I channel 

noise and estimate becomes 



Returning this result to (5-12), which describes the QPSK 
equivalent noise variance yields 
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2 


2 \^Le =%l 1 _ 2 , 

4n 2 E k 


Figure 5-1 shows the equivalent noise variance for BPSK and QPSK 
using (5-6) and (5-18). Note that at low SNR's the variance for 
QPSK is less than that of BPSK. This effect is known as "self noise" 
and in decision directed loops such as this it refers to the effect 
that making bad decisions (which occur more frequently at low 
SNR s) has on the noise inside the loop. This will be demonstrated 
graphically when 8 and 16PSK self noise are discussed. For QPSK 
though, the self noise is prevalent only at SNR’s below about 5 dB. 
In fact, at higher SNR's the BPSK solution for equivalent noise 
variance can be used as an approximation. 

To calculate the equivalent noise variance for 8 and 16PSK it 
is first determined that no reduction in the complexity of (5-4) is 
possible, as was done with BPSK and QPSK. The data estimates and 
corresponding orthogonal noise components are interdependent. 
For example, a noise sample on the I channel could move the 
received data point vector across a sector boundry and cause an 
error on the estimate of what was transmitted on the Q channel. 
This is shown for the 8PSK constellation in Figure 5-2. The hash 
marks on the vertical axis denote the possible transmitted 
quadrature channel data. 
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SNR (E/N ) in dB 

O 


Figure 5-1. BPSK and QPSK equivalent noise variance. 
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Therefore to calculate the variance of the equivalent noise 
for 8PSK and 16PSK a direct approach was utilized. The variance 
was numerically calculated by evaluating the triple integral 


=fi[w,(() 2 ]= J S j s (N„N e ,e m )f(N„N a ,e m )dN,dN Q de m 


(5-19) 


where 


g(N„N e ,e m ) = [N,Q-N a i } 2 (5-20) 

is the square of the equivalent noise and f(N h N Q ,6 m ) is the joint 

density over the noise components and the discrete transmitted 
modulation angle. Each of the noise components and the 
transmitted modulation angle are all independent of each other 
and thus the joint density may be expressed as 

f(N„N a ,e m )=f(N,)f(N e )f{e m ). (5-21) 

If 8PSK modulation is assumed and it is assumed that there are 
equally likely transmitted symbols, then integrating out the 
modulation from (5-19) results in 
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(5-22) 


< = I 1 II 

i — 0 — oo — oo 

Since the double integral of (5-22) will yield the same result for 
each of the eight modulation points it is possible to reduce (5-22) 
to 


oJ,-J jg(N„N Q .0 m =O)f(N,)f(N Q )dN,dN e (5-23) 


where the zero-angle transmitted modulation point has been 
assumed for mathematical convenience. To proceed a change of 
variables is performed to move the joint Gaussian probability 
"mountain" out to the modulation point corresponding to d m =0. 

This allows the double integral to be performed over sectors that 
correspond to the decision regions for the data estimator. 1 The 
new double integral is 


7 {fi2j«2j 

<y 2 N e = z \ J J g( x - y. e rn = 0)f(y)f(x - 1 )dydx > 
j =0 [Pij a\j 


(5-24) 


where the point (jc = i,y = 0) corresponds to the modulation point on 
the positive I axis. The limits of integration correspond to the 
values of x and y that are in each of the eight decision regions. To 


!The thick dashed lines and the axes in Figure 5.2 form the boundries of the 
decision regions. 
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proceed with performing this numerical integration it is helpful to 
view the components of the integrand as 3 -dimensional surfaces. 
Figure 5-3 shows the function g(x-l,y,e m =o), the square of the 

equivalent noise, for 8PSK. For comparison, the same function is 
plotted for BPSK, QPSK, and 16PSK in Figures 5-4, 5-5, and 5-6, 
respectively. For an SNR of 10 dB the joint Gaussian probability 
mountain, centered at the (x = 1,^ = 0) modulation point, is shown in 

Figure 5-7. The entire integrand of the double integral is shown 
for 8PSK at an SNR of 10 dB in Figure 5-8. The variance of the 
equivalent noise is the volume under the integrand surface. What 
the 3 -dimensional view contributes here is assistance in setting 
the limits of integration for the numerical calculation of (5-24). 
Before the use of surface analysis was employed it was not fully 
understood just how localized the surface prominences are. To test 
the technique the numerical approach was first applied to QPSK. 
This numerical data are shown in Figure 5-1. The source code for 
all M-ary PSK equivalent noise variance calculations and the data 
are given in Appendix D. Having verified the functionality of the 
approach for QPSK the technique was applied to 8PSK and 16PSK. 
The resulting data are shown in Figure 5-9 along with the closed 
form solutions to BPSK and QPSK. Note that at low SNR's the 
variance drops with increasing M. To understand how this is 
possible consider again the surface plots of the square of the 
equivalent noise. For BPSK this surface is a parabolic one, 
increasing with both positive and negative Y. For the other three, 
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the surface is discontinuous and, in particular, has channelized 
minima in regions that correspond to the increasing parabolic 
surface of BPSK. As M increases the number of minima increases 
and thus the volume under the square of the equivalent noise 
surface decreases. Correspondingly, the variance itself decreases. 
At high-SNR's, when the Gaussian probababilty mountain is very 
narrow, only the continuous, parabolic region (for all M) that 
contains the modulation point(* = l.r - 0) will provide significant 
contributions to the variance calculation. Thus, all of the MPSK 
equivalent noise variance calculations converge at high-SNR's. 



Figure 5-3. The 8PSK squared equivalent noise. 
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Figure 5-7. The joint 



ssian probability mountain at 10 dB SNR. 
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Variance of the Equivalent Noise 






s 2 ThP Gain O f Iha Itee Detector 


The next component of interest in the calculation of the 
variance of the phase error is the gain associated with the phase 
detector. In Chapter 4 the gain was shown to be a function of SNR 
through (4-27). Further, it was presented as the slope of the phase 
detector characteristic, measured at a phase error of zero, as 
shown in (4-30). Together these representations will allow the 
computation of the gain of the phase detector for 8PSK and 16PSK. 
As with the variance of the equivalent noise for 8PSK and 16PSK, 
the interdependence between the data estimate on the I channel 
with the noise on the a channel and the interdependence between 
the data estimate on the dchannel with the noise on the I channel 
prohibit the closed form solution of the phase detector gain. 
Before discussing the solutions for 8PSK and 16PSK, the solutions 

for BPSK and QPSK will be presented. 

Using (4-26) and (4-30) the phase detector gain is stated as 


a SNR ~ 


<9E[cos(0 e - 0 w )Q + sin(0 e - flj/] 1 


dd„ 


e e =o 


(5-25) 


Expanding the trigonometric functions yields 
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a SNR ~ 


(5-26) 


{cos( 9 e ) cos( e m ) + sin( 0 e ) sin ( G m )}q 
+{sin(fl,)cos(0„ ) - cos(0, )sin(»„, )}/ 

30, 


l0 e =O 


Writing the trigonometric functions of the transmitted modulation 
angle as I and Qdata components and combining like terms forms 


a{cos(0 e )£[/ M G-/a„]+sin(0 e )£[a I e+/ JB /]} 

a SNR - ^ 


90, 


(5-27) 


0,=O 


Distributing the derivative and applying its multiplicative 
property results in 


a SNR ~ ^ ImQ ~ I On 


9cos (gj 
90, 


+ cos(0 e ) 


3^i m Q-ig,\ 


0 e = 0 


90, 


8 e = 0 


+ E^Q m Q + I m I 


9sin(0j 


90, 


+ sin(0 e ) 


9E\Q m Q+lJ 


e ,= o 


90, 


0^=0 


(5-28) 


The first and fourth terms of (5-28) equate to zero leaving 
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a SNR ~ 



9E[iQ„] 

+£ M„,.„ +£ M 

d0 e 

de. 


n 

o 

o 

II 

Q5 


(5-29) 

For BPSK, all but the last term of (5-29) equate to zero because 
the data sent on the Q, channel and the estimate of what was sent 
on the Q_ channel are both always zero. The last term is the 
expected value of a function of random variables and is stated 

£[/„/]= I 1 Uf(U)f(N,)dN,dI m . (5-30) 


It is noted that the phase error has been set to zero before 
beginning the integration. Expressing the estimate of the I channel 
data as a function of the noise on the I channel and the particular 
transmitted data on the I channel makes it possible to integrate 
out the data and form 

= = l)f(N,)dN, -i =-l)f(N,)dN, . 

^ -oo -« 


(5-31) 
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Now, considering the two possible values for the data estimate 
transforms (5-31) into 


J OO | J 

«[V] = 2 jf(N,)dN,+ i )f(N l )dN l -±]f(N,)dN l . 

-1 -CO L -OO ^ 1 


(5-32) 


Using (5-14) and the "Q" function, which is related to the normal 
probability function and defined by 


Q(y)=j 


1 -^r 


V27T 


dx, 


(5-33) 


it is possible to express (5-32) as 



\ 

-1 

J 


-2* 


V °N, J 


+ f« 


'-i' 


1 a f>, ! 




V^/y 


(5-34) 


This reduces to 



( i > 

f . \ 


r \ 

Q 

— 

-Q 

l 

= 1-20 

1 


k g n, ) 


J 


K 0 "! ) 


(5-35) 


Again stating that for BPSK, (5-35) is the only non-zero term in 
the expanded equation for the phase detector gain, (5-29), and 
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noting the "Q!' function's relationship to the standard "erf 
function, namely 



1 -erf] 


V 


_y_ 

V2 


(5-36) 


it is possible to state 


a SNRgpsK 



(5-37) 


Inserting the expression for the variance of the I channel noise 
yields [37, pg. 91] 


a SNRgpsK 



(5-38) 


For QPSK, all of the terms of (5-29) are non-zero, Further, the last 
two terms are equivalent and can be stated 

1 1 (5_39) 

— OO — oo 


Following the same procedure as for BPSK, next the modulation is 
integrated out yielding 
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«M = ~ J s[n,.L 


2 2 __ 

1V2 7 




2 

V2 


(5-40) 


f[Nj)dNj 


Now, considering the two possible values for the data estimate 
transforms (5-40) into 


E\U 


1 V2 V2 

2 2 2 




_V2 

1 V2 -V2 r 2 

2 2 2 


J/JtyJlW, 


_V2 
2 

| /(W;)<W( + iz^2| jV(Ar,)<w, 

T 


(5-41) 


Once again applying the "Q" function, the result is 


: M 4 e 


-V2 "1 

--Q 

f ^ 1 

+ -Q 

f — V 2 1 

-ic 

f V2 1 

l 2<T *7 J 

4 

l 2<7 ", ; 

4 

v 2cT "/y 

4 

v 2cT "/> 


. (5-42) 


This can be reduced to 


eM=|-2 


' V2 A 
( 2<T » ; 


(5-43) 


or using the "erf" function with the I channel noise variance 
substituted in as 
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(5-44) 



To complete the phase detector gain analysis for QPSK requires 
that the first two terms of (5-29) be examined. The procedure for 
calculating either one is similar. It is noted, however, that the 
phase error cannot be set to zero until after the derivative is 
evaluted. Therefore to proceed, consider the expected value 
function in the first term of (5-29), namely 

E[l m Q] = J I ]lm&{lm)f{Qm)f{H Q )dN Q dQ m dI m . (5*45) 


Next, integrate out the modulation and express the estimate of the 
Q_ channel data as a function of the noise and modulation data to 
form 




22 2 

1 1 V2 

22 2 


' , _V2 „ _VT| 

N o •> I m rs >Qm ^ 


— oo V 


*Q' l m 2 

V2 
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N Q ,I m =^-,Q m = 2 


1 1 -V2 7 , 
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1 1 -V2 
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' — s/2 ^2^ 


|/KH> 

|/K)<» e 




(5-46) 
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To continue the phase error, 6 e , must be considered. As was done 
previously, it is assumed that the phase error is small. Next 
consider that a phase error will rotate the received constellation 
in the decision space. Depending on whether the phase error is 
positive or negative and depending on what modulation data is 
transmitted the rotating effect can help or hinder the Q, channel 
data decision process. When the rotation assists the Q. channel 
decision, by moving the constellation point away from the I 
channel axis, it moves the point closer to the Q. channel axis and 
thus hinders the I channel decision. 

In terms of integrating n q , this phase error affects the 
boundry value in the decision space beyond which// G will cause 

an error. If the phase error has moved the constellation point 

away from the I channel axis, i.e., the Q, channel decision boundry, 
a greater value of n q is required to cause an error in the Q 

channel estimate. A positive phase error will cause a transmitted 
symbol in the first or third quadrant to appear closer to the I 
channel axis in the decision space. This is shown in Figure 5-10. In 
the second and fourth quadrants a positive phase error will mean 
the symbol appears farther from the I channel axis. 

Therefore, in the first quadrant the Q_ channel noise needs to 
be 


Nq < -sin 



\ 

) 


(5-47) 
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to cause an error in choosing Q. Correspondingly in the third 
quadrant an error will occur if 

If in the second quadrant the noise is 

N Q < _sin (x " = " Sin ( 4 + 

or in the fourth quadrant 

Nq > -*.)■ sin (f + ®' ) 

an error will occur. These are the four noise boundries required to 
expand (5-46) taking into account the phase error. When they are 
applied to (5-46) the result is 


(5-48) 


(5-49) 


(5-50) 
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Figure 5-10. The QPSK constellation with and without phase error. 








Applying the "Q!' function yields 




After combining complementary "CT functions (5-53) becomes 
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£ M]=H e 


(5-54) 


sin 



--+-QI 

4 2 1 



and finally 



(5-55) 


By a similar approach for the second term of (5-29) its 
expected value function becomes 



(5-56) 


Referring to (5-29) the QPSK phase detector gain is now stated 
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(5-58) 

To further reduce the partial derivatives in (5-58) it is necessary 
to take the derivative of an integral through the use of the 
relation 

^jf( x ,y)dx = j^[f( x ,y)}k + f(q,y)^-f(p,y)^. ( 5 - 59 ) 

Since the "Q" function integrand in (5-58) is not a function of the 
phase error, the first term in (5-59) will equate to zero. Further, 
since the upper limit on both "Q." function integrals is infinity, the 
second term of (5-59) will equate to zero. The result is that (5-58) 
becomes 
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(5-61) 


+ 1 - 2(2 


f V2 ' 


2<t 


tf / ) 


Evaluating (5-61) at a phase error of zero yields 


a SNR QPSK - 1 - 2 ( 2 | 
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(5-62) 


and finally 
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Using the "err' function and inserting the I channel noise variance 
(5-63) can also be stated [35] 


a SNRQPSK 


= erf 
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2N„ 


(5-64) 


This is the closed form solution of the phase detector gain for 
QPSK. The results for both the BPSK and QPSK closed form 
solutions are plotted in Figure 5-11. Note that the logarithm of the 
square of the gain is actually plotted. This will prove useful when 
the entire phase error variance equation is reconsidered in the 

next chapter. 
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Figure 5-11. BPSK and QPSK phase detector gain. 
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For 8PSK and 16PSK a closed form solution is unobtainable. 
A review of the components in (5-29) will reveal this. The data 
estimates, both / and Q, require that the expected value functions 
in each of the components of (5-29) be integrated over both N, 
and Nq because of the interdependence between each of the 

estimates and both noise channels. For each component a solution 
will arise involving the integration of a "Q" function. This cannot 
be solved in closed form. 

To circumvent this problem the gain is calculated directly 
from the phase detector characteristic itself. By considering a 
small phase error, at a specific SNR, and calculating the phase 
detector characteristic output, the gain can be obtained. This 
process uses the approximation that for small phase errors the 
phase detector characteristic has a linear response. By so 
calculating the detector's response to a small phase error and 
dividing that response by the error, the slope, i.e., the gain, is 
computed. 

To proceed, the phase detector characteristic of (4-27) is 
employed and restated here as 

pd( e e ) = £[sin(e c + e m -e m )]. ( 5 -65 ) 


The first step is to express the expected value function as a 
summation over both discrete random variables, d m and 6 m , 
yielding 
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(5-66) 


PD(d e )= X X sin (^ + ~ 6m ‘ ) p ( 0m ‘ ’ ^ e ) * 


i=0 j = o 


Next the mutual probability is decomposed resulting in 
PD(6 e ) = X X Sm ( 6 e + ~ 6m > I 0 " 1 ' ’ ) * 


(5-67) 


i=o y=o 


The expression of conditional probability in (5-67) describes the 
probability of receiving one of the M possible modulation angles 
when a particular modulation angle is transmitted. The apriori 
probability of the transmitted modulation angle in (5-67) equates 

to -1 when the transmitted symbols are all equally likely thus 

M 

producing 

PDi»e)=-b X S“"(e. + 9»; -8-.WS (5 ' 68) 

M i=0 7=0 v 


Because the constellation is circularly symmetric each element of 
the sum over the transmitted modulation angles will be the same. 
Since there are M of these elements it is possible to reduce (5-68) 
by fixing the transmitted modulation angle, e.g., at d m -0, and 

then replacing the outer summation by multiplying through by M. 
In this way (5-68) becomes 
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(5-69) 



What remains is to describe the conditional probability of 
receiving a particular modulation angle estimate given that0 m =o. 

This is accomplished by integrating the conditional density that 
describes the probability of receiving a particular phase, 7, given 
0 m =0, over the decision region of phase that corresponds to the 

received modulation angle. The density is configured to account 
for a phase error, 6 e , and is expressed [41] 


p(y \o m =o,e e ) = - 


E 

, *0 
2n 


, R En , „ , cos 2 ( y+ e e ) 

1 + J cos [y+6 e )e N ° 


N n 


xQ 


-Jjf cos (r+0 e ) 


(5-70) 


As was just mentioned, to compute each of the M conditional 
probabilities in (5-69) it is necessary to integrate the density of 
(5-70) over the decision region corresponding to each of the M 
possible received phases. 

Through the use of (5-69) it is possible to examine the phase 
detector characteristic for any M. Figure 5-12 shows the phase 
detector characteristics, using this approach, for M=2,4,8 and 16 
all at an SNR of 20 dB. For each M one detector cycle has been 
plotted. Note that all the characteristics have positive nearly 
equivalent slopes at zero phase error. This is because at 20 dB all 
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the PSK modes, with the exception of 16PSK, have a phase 
detector gain of unity. The gain for 16 PSK is slightly less than 
unity. Having equivalent gain follows when it is considered that at 
such a high-SNR the decision device in the loop is making near 
perfect decisions as to what was transmitted. As the SNR goes 
down the gain will go down. This is shown for each PSK mode in 
turn using Figures 5-13, 5-14, 5-15, and 5-16 for BPSK, QPSK, 

8PSK, and 16PSK, respectively. In each of these plots a family of 
phase detector characteristics is presented wherein each 
characteristic is plotted at a different SNR. As the SNR drops the 
maximum magnitude and gain drop. Note also that the maximum 
magnitude of the phase error that can be tolerated before the 
phase detector output swings the opposite way and causes a cycle 
slip, also decreases with decreasing SNR. This is distinguished by 
the movement of the phase error location, of the maximum 
magnitude of the characteristic, inward toward the lock point as 
the SNR goes down. Evaluating any of the characteristics in Figures 
5-13 through 5-16, at a very small phase error allows for the 
calculation of the gain of the phase detector characteristic, as 
mentioned before, simply by dividing the result by the phase 
error. For QPSK this data is plotted along with the closed form 
solutions in Figure 5-11. Once again this technique of comparing 
the numerical approach to the closed form solution was used for 
verification. In Figure 5-17 the 8PSK and 16PSK numerical results 
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are presented along with the BPSK and QPSK closed form solutions. 

All of the numerical calculations were conducted using the 
source code in Appendix E. The raw data are presented there as 
well. 



Figure 5-12. MPSK phase detector outputs at 20 dB SNR. 
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Figure 5-16. 16PSK phase detector characteristics at SNR's of 10, 

12, 14, 16, 18, and 20 dB. 





Figure 5-17. MPSK phase detector gain. 


4The Toop Bandwidth 


The third component of the phase error variance equation 
that will be discussed is the loop bandwidth. As was shown in 
Chapter 2, the transfer function of a classical second order 
phaselock loop can be expressed by using the gain terms in the 

loop as 

+ (5-71) 

H(I) - 


or by using the classical design parameters as 
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( 5 - 72 ) 


h(s)= 2 s2( “ , ; + 0> " 

s + s2L s co n + (X) n 


Further, its loop noise bandwidth can also be expressed in the 
same two ways. When the gain terms are used the double-sided 
loop noise bandwidth can be stated 


B = Vg£+f 


( 5 - 73 ) 


When the design parameters are used the result is 


B, = co 



( 5 - 74 ) 


just as described in Chapter 2. In Chapter 4 the transfer function 
for the second order high-SNR loop was derived and is restated 
here as 

H(s)= a sm K »( s + a ) ( 5 - 75 ) 

s 2 +a sm K„(s+a) (5 /5) 


It was then equated to the design parameter version of the 
transfer function by stating 


H(s) 


s 2 +2£'co' n s + a>' 2 


( 5 - 76 ) 
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where 


(On = V &SNR a n 


(5-77) 


and 


C' - V a SNR C 


(5-78) 


are the "dynamic" design parameters. If the loop is designed 
assuming an SNR operating condition that is high enough to make 
the phase detector gain equate to unity, then the actual natural 
frequency and damping factor of the loop will be exactly those of 
the classical design. If operation is then conducted at lower SNR's, 
where the phase detector gain is less than unity, the actual 
natural frequency and damping factor will decrease from the 
original design values. Among other phenomena, this has the 
effect of narrowing the loop noise bandwidth. Substituting the 
design parameters of (5-77) and (5-78) into the loop noise 
bandwidth equation of (5-74) yields 


This reduces to 
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B, = (£> r 


a smC + 


(5-80) 


At high-SNR's, when a SNR = 1 the loop bandwidth of (5-80) will 

equate to (5-74) just as expected. However at low SNR's when 
a sNR « 1 equation (5-80) becomes 


( 0 „ 

B,. .... »— s. 


How SNR 


4C 


(5-81) 


There is a transition from one loop bandwidth at high-SNR's to a 
narrower loop bandwidth at low SNR's. For all the PSK modes the 
transition occurs at higher SNR's as M increases. This is shown in 
Figures 5-18, 5-19, and 5-20 where the normalized bandwidth 
has been plotted for 3 different damping factors. Note that the 
damping factor plays a crucial role in determining the size of the 
transition. As £ increases the loop bandwidth will have a much 
larger transition resulting in a much narrower relative bandwidth 
at low SNR's. 

As has been shown, the loop bandwidth changes because the 
phase detector gain changes. Both the phase detector gain and the 
equivalent noise variance change because incorrect data decisions 
are being made. In the next Chapter the combined effects of all 
three of these components will be assimilated to calculate the 
variance of the phase error as a function of SNR. 
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Chapter 6 

PHASE ERROR VARIANCE RESULTS 

"Why Think ? Why not try the experiment ?” 
John Hunter 

anatomist & surgeon, 1 728 - 1 793 
6.1 Introduction 


In Chapter 5 all of the complex components of the phase 
error were investigated individually. In this Chapter they will be 
combined in the equation for the phase error variance of the 
MPSK carrier tracking loop that uses the high-SNR approximation 
to MAP estimation. Once again this very important equation is 
restated for convenience as 


2 _ b l < 

U 0 n • 

SR a SNR 


( 6 - 1 ) 


6.2 N-Phase Loss 

The first step in investigating the combined affect of each of 
the varying components of the right side of (6-1) is to reconsider 
the equivalent noise variance in light of the data obtained in 
Chapter 5. At high-SNR's the variance of the equivalent noise, for 
all M, can be approximated by the result for BPSK, namely, the 
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reciprocal of twice the signal-to-noise ratio. At low SNR’s, for M>2, 
the self-noise in the receiver lessens the equivalent noise 
variance. It is possible to express the equivalent noise at any SNR 
and for all M, in terms of the BPSK approximation and the self 

noise, as 


where | is the self noise and is a function of SNR. It follows that 
for BPSK the self noise term is unity for any SNR. As stated above, 
for M>2 at high-SNR's it will also equate to unity. At low SNR’s, 
when the self noise is significant, it can be calculated using the 
equivalent noise variance data in Chapter 5 and equation 6-2. 
Substituting this new expression for the equivalent noise 

variance into (6-1) yields 

2 (6-3) 

e ‘ 2 E SR (Xsnr 


The reciprocal of the third fraction of (6-3) has a special 
designation. It is known as the N-phase loss. 1 Its effect is to 
increase the variance as the SNR goes down and is expressed as 


1 riassicallv this term has used the name "squaring loss". Its origin comes 
fmm the thnes-Ntaop employed for BPSK where N - 2. In that loop the 
a eternal was sauared to remove the modulation and provide a 

known as squaring loss. 


155 


NL = 


(6-4) 


a SNR 

a • 

The variance of (6-3) now becomes 
ct 2 = K l b l J_ 

° e 2 E SR NL ' 


(6-5) 


Note that as NL gets smaller the variance goes up. This special 
term incorporates two competing effects for M>2. As the SNR goes 
down the self noise goes down. This should increase the N-phase 
loss and lessen the variance. However as the SNR goes down the 
phase detector gain goes down as well and this dominating effect 
lowers the N-phase loss overall, thus increasing the variance. In 
fact the N-phase loss has been approximated solely by the square 
of the phase detector gain for QPSK at reasonable SNR's of interest 
[35]. This approximation is valid for M>2 and is expressed as 

( 6 - 6 ) 

To see that this approximation is valid consider Figure 6-1 where 
both equations (6-4) and (6-6) have been plotted for M>2. Only 
for QPSK and 8PSK at low SNR's is the effect of self noise evident. 
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Figure 6-1. MPSK N-phase loss. 

63 The T <~>np Ban dwidth And The Symbol Ra te 

The second fraction of (6-5) is the ratio of the loop 
bandwidth to the symbol rate. As was mentioned in Chapter 4, for 
this type of loop this ratio is typically very small. To make the 
variance of the phase error data useful it is helpful to plot the 
data using the ratio of the high-SNR loop bandwidth to the symbol 
rate. That way loop designers can use the variance results at 
symbol rates of interest to them. This means that the variance can 

be expressed as 
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2 = Ho_.. HIGHSNR ^ 


B, 


2E 


B-snr 


SR NL 


(6-7) 


where fi SNR is the loop bandwidth compression factor. The 
normalized loop bandwidth data plotted in Figures 5-18 through 
5-20 are actually plots of the loop bandwidth compression factor 
for various MPSK using 3 different damping factors. While the 
compression factor does decrease the variance of the phase error, 
its effect is limited. It causes a maximum amount of compression 
of the bandwidth, at low SNR's, that is a function of the damping 
factor. It is expressed, using (5-74) and (5-80), as 


_ ( a SNR^C 2 + l) 
4f2+l) • 


( 6 - 8 ) 


When the phase detector gain becomes very small n SNR reaches a 
minimum value, corresponding to the maximum amount of 
compression to the loop bandwidth. It is of interest to note that 
once again, in terms of variance of the phase error, there is a 
trade-off between the slight beneficial effect of the compressing 
bandwidth, which lessens the variance with decreasing SNR, and 
the decreasing phase detector gain, i.e., N-phase loss, which 
increases the variance with decreasing SNR. The N-phase loss 
dominates in this trade-off since it decreases as the square of the 
phase detector gain and since the compression factor reaches a 
minimum value while the N-phase loss continues to decrease. 
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£4 The Variance Of The Phase Error 


Having isolated and considered the combined effects of self 
noise and phase detector gain as one term, i.e., N-phase loss, and 
having discussed the changing loop bandwidth-to-symbol rate 
ratio, it is now prudent to evaluate the variance of the phase 
error. Referring to equation (6-7) the variance can be plotted 
us ing the data obtained in Chapter 5. The variance data is 
presented for BPSK, QPSK, 8PSK and 16PSK in Figures 6-2, 6-3, 6- 
4, and 6-5 respectively. The data are presented for various high- 
SNR loop bandwidth-to-symbol rate ratios all with a damping 
factor of C = l and is tabulated in Appendix F. 



Figure 6-2. BPSK theoretical variance of phase error. 
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Figure 6-3. QPSK theoretical variance of phase error. 



Figure 6-4. 8PSK numerical variance of phase error. 
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Figure 6-5. 16PSK numerical variance of phase error. 

As expected the variance increases with diminishing SNR. At 
high-SNR's when the N-phase loss and loop compression factor are 
unity, the variance of the phase error can be approximated by 

2 _ Nn ^LhighSNR (6-9) 

highsnr 2 E SR 

At low SNR’s when the N-phase loss begins to decrease more 
rapidly the variance increases faster than the SNR decreases. This 
creates a thresholding effect where the variance curves begin to 
rise upward very quickly. There are several other characteristics 
of these curves that warrant disscussiori as well. 
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Notice first that at high-SNR's there is a linear trade-off 
available to the loop designer in pursuit of maintaining a specific 
variance of phase error. The trade-off is between loop bandwidth 
and SNR. This trade-off is evident in equation (6-9) as well. If it is 
desired to lower the operating SNR by 3 dB, while still maintaining 
the same variance, the designer need only narrow the loop 3 dB 
(assuming other system constraints permit this option). 

Unfortunately at low SNR's this trade-off is not linear. 
Narrowing the loop bandwidth 3 dB does not allow the designer to 
decrease the SNR 3 dB. This affect is demonstrated by considering 
the compression of the spacings between the variance curves in 
the thresholding region. Moving right from one variance curve to 
another along a fixed horizontal line of variance represents a 3 dB 
change in high-SNR loop bandwidth. As stated above, at high- 
SNR's this move also corresponds to a 3 dB drop in SNR. However, 
at low SNR's, specifically in the thresholding region, the move 
from one curve to another corresponds to less than a 3 dB 
decrease in SNR. The trade-off has diminished. 

Another interesting characteristic of the thresholding effect 
is that it is a strong function of the SNR and a relatively weak 
function of the loop bandwidth. This is attributed to the reliance 
of the N-phase loss on SNR. It is noted graphically by first 
considering again that at any SNR a 3 dB change in loop 
bandwidth corresponds to a 3 dB change in variance (moving 
vertically from one curve to another). Next consider that, in the 
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threholding region, a 3 dB change in SNR can correspond to a 
change in variance of several orders of magnitude. Thus the 

variance is strongly dependent on the SNR. 

A further, highly significant characteristic suggested by this 

data is that SNR’s, where power-limited channel operation has 
been conducted for QPSK, may be too low to support 8PSK and 
especially 16PSK operation. This is the case because the SNR at 
which the thresholding region occurs increases with M. As was 
mentioned in Chapter 1, there is strong interest in using coded 
8PSK and 16PSK systems to increase performance in limited 
channels. If the maximum available SNR is near or below the 
threshold region of the variance curves, then operation may not 
be possible because the variance does not meet operational 
specifications. 2 Further, as will be shown with the following 
simulation and hardware results, the loop may actually have 
trouble maintaining phaselock at SNR’s in the threshold region. 

AS Simula tion Results 

To verify that the numerical results for the 8PSK and 16PSK 
variance of the phase error were correct, simulations of the loop s 
operation were conducted. These simulations were conducted 
using the nonlinear baseband model of Figure 4-2. The code and 
raw data from the simulations are in Appendix G. The BPSK and 

2a variance of -35 db-radians 2 represents an approximate RMS jitter of 1 
while -25.6 db-radians 2 corresponds to about 3°. 
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QPSK simulation data match the theoretical results and thus 
confirm the simulation technique is accurate. The 8PSK and 16PSK 
simulation results are therefore assumed to be correct and indeed 
they confirm the numerical approach of calculating the variance is 
accurate when they are plotted together. 

The simulation data are shown in the following 12 plots. The 
first three plots, Figures 6-6, 6-7, and 6-8, are for BPSK at three 
different high-SNR loop bandwidth-to-symbol rate ratios. Note, 
once again, that the reference here to high-SNR loop bandwidth 
corresponds to the case where the phase detector gain is unity. 
This loop bandwidth can also be referred to as the "design" loop 
bandwidth in that the classical loop design parameters are the 
operational parameters only at high-SNR's. The selected loop 
bandwidths were chosen to match the available hardware 
configurations. The simulation data are plotted with a 95% 
confidence interval. 

Figures 6-9, 6-10, and 6-11 are for QPSK at the same design 
loop bandwidths as BPSK. Following these are Figures 6-12, 6-13, 
and 6-14 for 8PSK and Figures 6-15, 6-16, and 6-17 for 16PSK all 
reflecting the same bandwidths as that of BPSK and QPSK. 
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Figure 6-6 BPSK variance of phase error theoretical and 

simulation data. 
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Variance in dB-radians 



Figure 6-7 BPSK variance of phase error theoretical and 

simulation data. 
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Figure 6-8 BPSK variance of phase error theoretical and 

simulation data. 
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Figure 6-9 QPSK variance of phase error theoretical and 

simulation data. 





Variance in dB-radians 



6-10 QPSK variance of phase error theoretical and 
simulation data. 



Variance in dB-radians 



Figure 6-11 QPSK variance of phase error theoretical and 

simulation data. 






6-12 8PSK variance of phase error numerical and 
simulation data. 
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Figure 6-14 8PSK variance of phase error numerical and 

simulation data. 
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e 6-16 16PSK variance of phase error numerical and 
simulation data. 
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Figure 6-17 16PSK variance of phase error numerical and 

simulation data. 


6.6 Hardware Results 

As a third means of quantifying the variance of the phase 
error an MPSK high-SNR carrier tracking loop was constructed in 
hardware. The hardware has limits on the operational parameters 
it can support. It only has 64 loop bandwidth configurations. A 
specific loop bandwidth is selected by controlling the gains in the 
loop. Three of these bandwidths were selected over the range of 
operation of the hardware. They are the bandwidths used for the 
simulations mentioned in Section 6.5. The next 12 plots are the 
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hardware results for the variance of phase error. As with the 
simulations they are presented in sequence for each modulation 
scheme. The BPSK data are plotted in Figures 6-18, 6-19, and 6-20 
and the QPSK data in Figures 6-21, 6-22, and 6-23. Figures 6-24, 
6-25, and 6-26 are the 8PSK data and Figures 6-27, 6-28, and 6- 
29 are the 16PSK data. These data are tabulated in Appendix H 
along with a discussion of the hardware operation. 



Figure 6-18 BPSK variance of phase error theoretical and 

measured data. 
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Figure 6-19 BPSK variance of phase error theoretical and 

measured data. 
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Figure 6-20 BPSK variance of phase error theoretical and 

measured data. 
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Variance in dB-radians 



Figure 6-2 1 QPSK variance of phase error theoretical and 


measured data 
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Figure 6-22 QPSK variance of phase error theoretical and 


measured data. 
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Figure 6-24 8PSK variance of phase error numerical and 

measured data. 
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Figure 6-25 8PSK variance of phase error numerical and 

measured data. 
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Variance in dB-radians 



Figure 6-26 8PSK variance of phase error numerical and 

measured data. 
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Figure 6-27 16PSK variance of phase error numerical and 

measured data. 
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Figure 6-28 16PSK variance of phase error numerical and 

measured data. 




Figure 6-29 16PSK variance of phase error numerical and 

measured data. 
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ft 7 Discussion Of The Sim ulation And Hardware Results 


As the simulation and hardware data plotted in Figures 6-6 
through 6-29 indicate there is substantial verification of the 
numerical approach used to calculate the variance of the phase 
error. With the exception of a few samples, all of the simulation 
data is within a 95% confidence interval of the theoretical and 
numerical solutions of the phase error variance. Before the 
validity of the hardware results can be discussed, an interesting 
phenomena in that data must be addressed. This phenomena is 
the minimum measured phase jitter that is represented by the 

flair in the hardware variance curves. 

The presence of a phase jitter floor makes sense when it is 
considered that the hardware cannot obtain infinitely small 
variance of phase error data. There are several causes for this. 

First, the NCO uses direct digital synthesis to create its sinusoid 
outputs. As such it cannot exactly recreate the desired frequency. 
Related to this is the minimum phase error that will cause the 
minimum possible error signal to be created at the output of the 
EPROM's. This minimum error is approximately 0.5° measured 
along the constellation circle from the I or Q. axis. There is the 
added affect of the variance measuring device, i.e., the modulation 
analyzer. This device digitizes the I and Q. analog baseband data to 
creates samples that are then used to calculate the variance. 
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This variance noise floor in the hardware data dominates 
the results at high-SNR's only. At low SNR's the variance in the 
phase error data is accurately reflecting the mechanism of 
interest, i.e., the overall loss in performance due to making 
incorrect data decisions in the carrier tracking loop. When the 
MPSK theoretical and numerical data are "corrected" to account for 
the phase noise floor present in the hardware the results are 
accurate at all SNR's. 

One last noted characteristic that is present in the 8PSK and 
16PSK data from both the simulation and the hardware results is 
a minimum SNR below which phaselock cannot be maintained. 
While quantifying this phenomenon analytically requires an in- 
depth study of non-linear loop operation, some conclusions can be 
drawn from the simulation and hardware results since they are 
not affected by the constraints of linear analysis. 

In the plots of simulation data and hardware data for 8PSK 
and 16PSK, the lowest SNR at which a data point is plotted 
represents the minimum SNR for which phaselock can be 
maintained. The striking result in this data is that the SNR where 
phaselock terminates seems to be a weak function of loop 
bandwidth and a strong function of SNR. For both 8PSK and 16PSK 
the simulation and hardware data cover a 5.4 dB range of high- 
SNR loop bandwidth-to-symbol rate ratios (from 0.19% to 0.66%). 
However, the phaselock threshold moves at most 1.6 dB for the 
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8PSK hardware results and not at all for the 16 PSK simulation 
results. 
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Chapter 7 

CONCLUSIONS 

"Thinking makes it so. " 

William Shakespeare 
dramatist and poet, 1564-1616 


7.1. Discussion 

Quantifying the phase error variance of the carrier reference 
in an 8PSK or 16PSK receiver utilizing the high-SNR approximation 
to MAP estimation in a carrier tracking feedback loop was the 
motivating problem behind this research effort. This was 
accomplished by three different methods, including theoretical 
and numerical analysis, simulation, and laboratory hardware 
testing. All of these methods first required an understanding of 
the phase error variance and its origins. Chapter 2 is a discussion 
of what phase error variance is in a simple phaselock loop. This 
information alone is insufficient to understand the phase error 
variance in the phaselock loop of interest since it is a loop that 
tracks a DSB-SC MPSK signal as opposed to a CW carrier. The 
mechanism utilized to track the DSB-SC MPSK signal is an 
approximation to the optimum solution in the MAP sense. This 
approximation, i.e., the high SNR approximation, is a commonly 
employed and studied technique as outlined in Chapter 1. Chapter 
3 provides a discussion of the MAP estimator approach and the 
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high-SNR approximation approach to estimating carrier phase. 

Riter and Hon [26] have previously discussed the optimum MAP 
estimator for BPSK while Simon [32] discussed the optimum MAP 
estimator for QPSK and 8PSK. Simon also presented the high-SNR 
loop structure for MPSK. The 16PSK MAP estimator solution is 
presented here in this research effort following the approach of 
Simon. Using the integro-differential equation for the high-SNR 
loop, the MPSK variance of phase error was derived in Chapter 4. 
The technique utilized for this follows the classical approach (see 
Chapter 2) and was verified by the QPSK variance of phase error 

derivations of Lindsey and Hinedi [35]. 

Quantifying the variance of phase error for the BPSK and 
QPSK high SNR loop has been accomplished [37, pg. 91] [35] using 
the independence of the quadrature data estimates on their 
respective orthogonal noise components. Because of only recent 
increased interest in 8PSK and 16PSK and because the above- 
mentioned independence does not exist for these M>4 PSK 
schemes, the variance of phase error had not previously been 
quantified for M>4 in the high-SNR loop. The exact problem 
concerned the calculation of two components of the phase error 
variance, namely, the variance of the equivalent noise in the loop 
and the phase detector gain. A third previously uncalculated 
component of the phase error variance, the loop bandwidth, is a 
function of the phase detector gain. In Chapter 5 the calculation of 
all three of these components is accomplished. The approach is 


193 


verified using the previously published BPSK and QPSK results 
mentioned above. This theoretical analysis and set of numerical 
solutions represented the first method used to obtain the variance 
of phase error for 8PSK and 16PSK. The resultant data are shown 
in Chapter 6. 

Simulations of the high-SNR loop were conducted to verify 
the results of the numerical solutions for 8PSK and 16PSK. The 
nonlinear baseband model of the loop was utilized in these 
simulations and the method verified by first simulating BPSK and 
QPSK. This second method of obtaining the variance data confirms 
that the numerical solutions for the components of the phase error 
variance and the variance itself are accurate. 

The third method of obtaining variance of phase error data 
is in some ways the most important. An MPSK carrier tracking 
loop was constructed in hardware and used to obtain laboratory 
measurements of the variance of phase error. This is actual 
performance data as opposed to either analysis or simulation data. 
The hardware environment was designed so that the variance 
measurements reflected, as closely as possible, only the effects 
studied in this research. This was accomplished by removing other 
unknowns from the test configuration that would influence the 
variance such as symbol timing jitter, intersymbol interference, 
and received signal automatic gain control jitter. As the results of 
Chapter 6 and the discussion of Appendix H indicate, the 
hardware data verified both the numerical and simulation results. 
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Having quantified and verified the phase error variance for 
8PSK and 16PSK and having plotted it as a function of both SNR 
and loop bandwidth it was possible to draw certain conclusions. 
Chapter 6 contains a discussion of these observations and so they 
are only summarized here. 

The data show that at high SNR's the variance is directly 
related to both the SNR and loop bandwidth and can be obtained 
by the use of the approximation in (6-9). It should be noted that 
this form of the equation for variance closely follows the equation 
for variance of the basic PLL, namely (2-52). In fact equating the 
IF bandwidth to the symbol rate yields an exact match. This 
follows esoterically, too, when it is considered that at high SNR's 
the high-SNR MPSK loop makes good data decisions and produces 
an error signal that exactly matches that of the basic PLL. It is at 
low SNR's that the variance performance of the high-SNR MPSK 
loop is poorer than that of the basic PLL. 

At low SNR's the variance begins to increase more rapidly 
than the SNR decreases. Further the high-SNR trade-off between 
loop bandwidth and SNR is degraded. Another interesting 
phenomenon is that the SNR location at which the variance begins 
to increase rapidly, the threshold region, cannot be moved 
appreciablely by a judicious narrowing of the high-SNR design 
loop bandwidth. The data also demonstrate that the location of the 
thresholding region may cause difficultly in systems attempting to 
use M>4 PSK schemes if operation is power limited. For 8PSK this 
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threshold occurs at an SNR of about 8dB and for 16PSK it occurs at 
about 15dB. A case in point is that of NASA's TDRSS BPSK/QPSK 
links. Typically they operate with SNR's on the order of 10 dB. If 
8PSK is employed, as is being considered at this time, and the 
high-SNR loop structure is utilized, not only may the variance be 
excessive but the simulation and hardware data suggest that 
phaselock performance may be marginable. Further, this data 

suggests successful 16PSK operation would be unobtainable at 10 
dB SNR. 

Having considered some of these unique characteristics of 
the 8PSK and 16PSK variance of phase error data, quantified in 
this study, it is clear that coherent receiver design for these 
schemes must include a thorough analysis of the carrier tracking 
loop to be implemented. The high-SNR approximation is a popular 
choice, however, operation may not be possible if relatively high 
SNR conditions are not available. Designs involving the use of 

coded M-ary PSK schemes to increase performance in limited 
channels are a prime example. 

7.2. Where To Next? 

The use of the high-SNR approximation is just that: an 
approximation. The first obvious solution to improving the 
variance performance is to study the use of the optimum MAP 
estimator itself. It may yield superior performance at the SNR's of 
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interest and push the threshold region to the left, i.e., to a lower 
SNR. Further, the barrier to implemention, namely analog 
hardware construction of hyperbolic operators, no longer exists, 
since digital receivers are now commonplace. In point of fact, this 
is exactly the research avenue being pursued by the 
telecommunications department at NMSU. Early results suggest 
that there is a small performance increase in using the optimum 
MAP estimator for 8PSK [42]. 

Another option for possibly increasing the performance of 
the M>4 PSK receiver is to use coding to increase the performance 
of the data estimators. This was suggested for 8PSK by 
Ungerboeck [43], Since the location of the threshold region is 
indirectly traceable to making bad data decisions, if the data 
decision process can be improved by coding gain, then the 
threshold region should move to the left. This idea is particularly 
attractive for systems where coding is already being employed, 
e g., gPSK and 16PSK TCM. Care must be exercised in utilizing this 
option to ensure that loop stability is not compromised by the 
required delay inherent in making estimates with such schemes. 

Another area of interest in these loops that warrants study 
is the acquisition process. Several factors make the M>4 loops 
particularly attractive for study in this regard. First, the number 
of stable phaselock points in the constellation increases with M. In 
fact, there are M phaselock points in an MPSK carrier tracking 
feedback loop. It is also true that there are M unstable points that 
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generate zero phase error. It follows that the maximum phase 
error that can be tolerated drops with increasing M. 

Another important factor in the acquistion process is the 
loop bandwidth. It has a strong relationship to acquistion 
parameters such as mean-time-between-cycle slips and, as was 
demonstrated, the loop bandwidth is narrowed as the SNR drops. 
Other factors, such as lock-in range and hold-in range, would also 
be affected since the loop bandwidth is changing. 
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Appendix A 

SIMULATION COMPARISON BETWEEN HIGH SNR MAP LOOP AND 
LECLERT & VANDAMME VARIATION 


To assist in making the decision as to which type of tracking 
loop to consider for the study of phase error variance two loop 
designs were compared using 8PSK simulations of each. To 
compare the two loops they were simulated at the same loop 
bandwidth over a range of SNR's. A static loop bandwidth was 
maintained in each loop by varying the phase detector gain. These 
early simulations were written in C. The source code for the 
Leclert & Vandamme simulation is: 

/*Leclert & Vandamme 8PSK simulation of earner recovery 
with fixed loop bandwidth-to-symbol rate ratio / 


#include <stdio.h> 
#indude <stdlib.h> 
#include <math.h> 


main() 

{ 


/♦declare system parameters*/ 


int k; 
int m; 
int j; 
int s; 
int n; 
int i; 

int uncore; 
double errk; 

double phsk; 
double refkl; 
double pdifk; 


/♦discrete time variable*/ 

/♦discrete time variable*/ 
/♦discrete time variable*/ 
/♦discrete time variable*/ 
/♦discrete time variable*/ 
/♦discrete time variable*/ 
/♦discrete time variable*/ 

/♦error signal out of phase 
detector at time k*/ 

/♦received signal phase at time k / 
/♦phase of VCO at time k+1*/ 
/♦phase difference between 
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double ik; 

double imk; 
double qk; 

double qmk; 
double ilk; 

double i2k; 

double a; 

double t; 
double kO; 
double ilkl; 

double sum; 

double wn; 

double damp; 
double bw; 
char dptr[10]; 

char pdgainsflO]; 

double endb; 
double pdgain; 
double phsfix; 
double sqrerr; 

double sumerr; 

double meanphs; 
double varphs; 


received signal & VCO at time k*/ 
/^received value for I sample at 
time k*/ 

/★generated I sample at time k*/ 
/*received value for Q. sample at 
time k*/ 

/*generated Q. sample at time k*/ 
/*initial condition of filter 
integrator at time k*/ 

/initial condition of VCO 
integrator at time k*/ 

/^inverse of the one pole filter 
time constant, 3dB break point*/ 
/*symbol period*/ 

/*the VCO slope in Hertz per Volt*/ 
/*hold variable for intial 
condition, II, at time k+1*/ 
/*averaging of phase detector 
error*/ 

/*user entered natural freq. 
of loop*/ 

/*damping factor*/ 

/*loop bandwidth*/ 

/*string pointer for entering 
floating point numbers*/ 

/*string pointer for entering 
floating point numbers*/ 
/*E(s)/N(o) in dB*/ 

/*phase detector gain*/ 

/*fixed initial phase*/ 

/*running sum of squared phase 
error refkl*/ 

/*running sum of phase error 
refkl*/ 

/*mean of phase error*/ 

/*variance of phase error*/ 


static double pdgmx[8] = {1.5425, 2.5, 4.0, 5.7, 7.3, 8.1, 9.0, 10.0}; 


/*declare data generation parameters*/ 


int d; 

double testd; 


/*declare QPSK test parameters*/ 

double xl; 
double x2; 

/*declare phase detector variables. These will be moved into a 
subroutine eventually*/ 
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double rcvdp; 
double ihatk; 


double qhatk; 


/*the phase of the received vector 
whose components are ik and qk*/ 
/*for this phase detector this is the 
output from the circular 
quantizer*/ 

/*for this phase detector this is the 
output from the circular 
quantizer*/ 


double errkhold; 
int count; 

/*declare noise model parameters, sigma is the WON standard 
deviation, nik and nqk are the gaussian noise samples at time 
k with mean zero and standard deviation sigma. / 


double sigma; 
long xconl; 
long xcon2; 
long x; 
long y; 

double sample 1; 
double sampleZ; 
double nik; 
double nqk; 

/*declare output file parameters*/ 


FILE *outfile; 
char *modehold; 
char *nameholdl; 


/initialize file output parameters*/ 


modehold = "w"; 

nameholdl = "Vandam"; 

outfile = fopen(nameholdl, modehold); 

/*initialize uniform RV seed for repeated tests*/ 

srand(2546); 
y = 47; /*seed*/ 

/*initialize phase step*/ 

phsfix = 0.0; 

/*input natural frequency, damping factor, and Es/No for 
loop*/ 
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printf("This is a second order carrier loop synchronizer 
simulator\n"); 

printf("\n"); 

printf( "Please enter the natural frequency of the loop 
transfer functionin''); 

printf("\n"); 

scanf("%s",dptr); 

wn = strtod(dptr,NULL); 

printf("\n"); 

printf("Now please enter the loop damping factorin''); 

printffin"); 

scanf("%s",dptr); 

printf("\n"); 

damp = strtod(dptr,NULL); 
printf("\n"); 

bw = wn*(damp + 1.0/(4.0*damp)); 
a = wn/(2*damp); 

/*determine sample interval that is required to uncorrelate 
phase error variance calculation*/ 

uncore = floor(3.1415927/wn/4.0); 

/initialize phase detector gain counter*/ 

i = -1; 

printf("\n"); 

printf("The loop noise bandwidth is: %f Hz\n",bw); 
printf("\n"); 

printf("The loop filter time constant is: %f sec\n",a); 
printf("\n"); 

fprintf(outfile,"8-PSK reconstruction loop results\n"); 
fprintf(outfile,"\n"); 

fprintf(outfile," Initial phase step = %f\n",phsfix); 
fprintf ( outfile \n " ) ; 

fprintf(outfile, "Noise random seed = %ld\n",y); 
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fprintf ( outf ile, " \n" ) ; 
fprintf(outfile,"wn = %f\n",wn); 
fprintf(outfile,"\n"); „ . 

fprintf (outfile, "damping factor = %f\n ,damp), 

fprintf(outfile,"The loop noise bandwidth is: %f Hz\n",bw); 

fprintf(outfile,"\n"); e „v _ v 

fprintf( outfile, "The loop filter time constant is: ^f sec\n ,a), 

fprintf ( outf ile, " \n" ) ; 

fprintf( outfile, "The number of skipped samples tor 
correlation is: %d \n", uncore); 
fprintf (outfile, "\n"); 

/*begin simulation*/ 


for(j=10;j<=24;j=j+2) 

{ 

i = i+1; 

sigma - ' 1 . 0 /sqrt( 2 . 0 *pow( 10 . 0 ,endb/ 10 . 0 )); 
/*initialize system parameters*/ 


ilk = 0.0; 
i2k = phsfix; 
refkl = i2k; 
t = 1.0; 

errkhold = 0.0; 
phsk = 0.0; 
errk = 0.0; 
xcon2 = 16807; 
xconl = 2147483647; 
sqrerr = 0.0; 
sumerr = 0.0; 
s = 0; 

pdgain = pdgmx[i]; 

k0 = 2 *damp*wn/pdgain; 


/*for RV generation*/ 
/*for RV generation*/ 


for(k=0;k<40;++k) 
for(m = 0;m<500;++m) 

/*construct received phase, noise, and transmitted data / 
/*running variance calculation*/ 

if (k > 19) 

* if (fabs(refkl) < 0.392699) 

{ 

s =s + 1; 
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sqrerr = sqrerr + refkl*refkl; 
sumerr = sumerr + refkl; 


} 

for(n=0;n<uncore;++n) 

{ 

/*received phase*/ 

/*phsk = 0.001*k;*/ 

/*noise*/ 

/* x = rand(); 
y = rand(); 

sample 1 = (x/32767.0)*0.999999; 
sample2 = (y/32767.0)*0.999999; 

*/ 

x =ldiv(y,xconl).rem; 
x=labs(x*xcon2); 

y =ldiv(x,xconl).rem; 
y=labs(y*xcon2); 

samplel = (x/2147483647.0)*0.999999; 
sample2 = (y/2147483647.0)*0.999999; 

nik = sigma*sqrt(-2.0*log(1.0-samplel)) 
*cos(2.0*3.14159265*sample2); 

nqk = sigma*sqrt(-2.0*log(1.0-samplel)) 
*sin(2.0*3.14159265*sample2); 

/transmitted data*/ 

d = rand(); /*d is now divided into 8 regions each 

representing a possible xmitted symbol*/ 


if(d <4096) 

{ 

imk = 0.923879; 
qmk = 0.382683; 

I 

else if((d >= 4096) && (d < 8192)) 

imk = 0.382683; 
qmk = 0.923879; 

} 

else if((d >= 8192) && (d < 12288)) 

{ 

imk = -0.382683; 
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qmk = 0.923879; 

else if((d >= 12288) && (d < 16384)) 

imk = -0.923879; 
qmk = 0.382683; 

else if((d >= 16384) && (d < 20480)) 

imk = -0.923879; 
qmk = -0.382683; 

else if((d >= 20480) && (d < 24576)) 

imk = -0.382683; 
qmk = -0.923879; 

else if((d >= 24576) && (d < 28672)) 

imk = 0.382683; 
qmk = -0.923879; 

} 

else 

{ 

imk = 0.923879; 
qmk = -0.382683; 

} 


/*end generation part of the simulation, this generated 
information will be used for this value of the simulation 
counter only*/ 

/*subtract the reference phase from the received phase*/ 


pdifk = phsk - refkl; 


^calculate the received symbol values for 1 and Qbased on 
the generated data and noise and the phase difference pdifk / 


ik = imk*cos(pdifk) - qmk*sin(pdifk) + nik; 
qk = imk*sin( pdifk) + qmk*cos(pdifk) + nqk; 


/*Now implement a phase detector to generate a phase error 
signal errk. This implementation performs a circular 
quantization on I and QfoUowed by some algebra to develop 
the error signal.*/ 


^calculate the phase of the received vector. The atan2 
function adjusts for the appropriate quadrant.*/ 


rcvdp = atan2(qk,ik); 
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/*quantize rcvdp into ihatk and qhatk*/ 

if ((rcvdp >= 0.0) && (rcvdp < 0.785398)) 

/* 0 <= rcvdp < Pi/4 */ 

I 

ihatk = 0.923879; 
qhatk = 0.382683; 

} 

else if ((rcvdp >= 0.785398) && (rcvdp < 1.570796)) 
/* Pi/4 <= rcvdp < Pi/2 */ 

{ 

ihatk = 0.382683; 
qhatk = 0.923879; 

{ 

else if ((rcvdp >= 1.570796) && (rcvdp < 2.356194)) 
/* Pi/2 <= rcvdp < 3Pi/4 V 

f 

ihatk = -0.382683; 
qhatk = 0.923879; 

I 

else if ((rcvdp >= 2.356194) && (rcvdp < 3.141593)) 
/* 3Pi/4 <= rcvdp < Pi V 
I 

ihatk = -0.923879; 
qhatk = 0.382683; 

1 

else if ((rcvdp < -2.356194) && (rcvdp >= -3.141593)) 
/* -3 Pi/4 < rcvdp >= -Pi V 
1 

ihatk = -0.923879; 
qhatk = -0.382683; 

} 

else if ((rcvdp >= -2.356194) && (rcvdp < -1.570796)) 
/* -Pi/2 > rcvdp >= -3Pi/4 V 


ihatk = -0.382683; 
qhatk = -0.923879; 
1 
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else if ((rcvdp >= -1.570796) && (rcvdp < -0.785398)) 
/* _pi /4 > rcvdp >= -Pi/ 2 V 

ihatk = 0.382683; 
qhatk = -0.923879; 

else if ((rcvdp >= -0.785398) && (rcvdp < 0.0)) 


/* 0 > rcvdp >= -Pi/ 4 */ 


ihatk = 0.923879; 
qhatk = -0.382683; 

} 

^calculate the phase error errk as a function of ik, qk, ihatk, 
and qhatk.*/ 

errk = ((qk-qhatk)/fabs(qk-qhatk))*(ik/fabs(ik)) - 

( ( ik-ihatk) /f abs ( ik-ihatk) )* ( qk/ f abs ( qk) ) ; 

/*This is the end of the phase detector*/ 


/*The next step considers the loop filter and VCO m one 
transfer function, as a difference equation. It calculates the 
phase of the VCO, refkl, for the next simulation iteration, 
refkl is based not only on the error signal, errk, but also on 
the state of the integrators that make up the one pole low-pass 
filter and the VCO. The filter can be represented by the 
function, 1 + a/s. a is one-over-the time constant, RC, of the 
filter The VCO xfer function can be considered as kO/s. kO is 
""describing the output frequency of the VCO m terms 
of the input error voltage from the loopfiltenThe filte 
integrator initial condition at time k is ilk. The VCO integrator 
initial condition at time k is i2k. The integration period 1 j * 1 
discrete time these transfer funcitons can be implemented 
follows*/ 

refkl = i2k + kO*t*errk + ilk*kO*t + a*k0*t*t*errk/2.0; 


/*now reset the integrator intial conditions for the next 
iteration*/ 

i2k = refkl; 

ilkl = ilk + a*t*errk; 

ilk = ilkl; 


} 

} 
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} 

meanphs = (sumerr/s); 

varphs = sqrerr/s - meanphs*meanphs; 

fprintf(outfile,"\n"); 

fprintf(outfile,"E(s)/N(o) = %4.1f, mean of error = %10.7f, 
variance of error = %10.7f\n",endb,meanphs,varphs); 

fprintf(outfile,"The phase detector gain used was: 
%f\n",pdgain); 

fprintf(outfile,"The number of samples used in the previous 
statistical analysis is: %d\n",s); 

fprintf(outfile,"\n"); 

printf("Es/No = %4.1f, mean = %f, variance = %f, samples = 
%d\n " , endb, meanphs, varphs, s ) ; 

} 

fclose(outfile); 

} 


The source code for the high-SNR loop simulation is: 


/*high SNR 8PSK simulation of carrier recovery with fixed 
loop bandwidth-to-symbol rate ratio*/ 

#include <stdio.h> 

#include <stdlib.h> 

#include <math.h> 

main() 

{ 


/*declare system parameters*/ 


int k; 

/*discrete time variable*/ 

int m; 

/*discrete time variable*/ 

int j; 

/*discrete time variable*/ 

int n; 

/*counter*/ 

int s; 

/*number of samples used in 
variance calculation*/ 

int uncore; 

/*number of samples to skip over 
for variance calculation*/ 
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double errk; 

double phsk; 
double refkl; 
double pdifk; 

double ik; 

double imk; 
double qk; 

double qmk; 
double ilk; 

double i2k; 

double a; 

double t; 
double kO; 
double ilkl; 

double sum; 

double wn; 

double damp; 
double bw; 
char dptr[10]; 

char pdgains[10]; 

double endb; 
double pdgain; 
double phsfix; 
double sqrerr; 

double sumerr; 

double meanphs; 
double varphs; 


/♦error signal out of phase detector 
at time k*/ 

/♦received signal phase at time k*/ 
/*phase of VCO at time k+1*/ 
/♦phase difference between 
received signal and VCO at time k*/ 
/♦received value for I sample at 
time k*/ 

/♦generated I sample at time k*/ 
/♦received value for Q. sample at 
time k*/ 

/♦generated Q. sample at time k*/ 
/♦initial condition of filter 
integrator at time k*/ 

/♦initial condition of VCO 
integrator at time k*/ 

/♦inverse of the one pole filter time 
constant (the 3dB break point)*/ 
/♦symbol period*/ 

/♦the VCO slope in Hertz per Volt*/ 
/♦hold variable for intial condition, 
II, at time k+1*/ 

/♦averaging of phase detector 
error*/ 

/♦user entered natural frequency 
of loop*/ 

/♦damping factor*/ 

/♦loop bandwidth*/ 

/♦string pointer for entering 
floating point numbers*/ 

/♦string pointer for entering 
floating point numbers*/ 
/*E(s)/N(o) in dB*/ 

/*phase detector gain*/ 

/*fixed initial phase*/ 

/*running sum of squared phase 
error refkl*/ 

/♦running sum of phase error 
refkl*/ 

/♦mean of phase error*/ 
/♦variance of phase error*/ 


/♦declare data generation parameters*/ 


int d; 

double testd; 

/♦declare QPSK test parameters*/ 

double xl; 
double x2; 
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/*declare phase detector variables. These will be moved into a 
subroutine eventually*/ 

double rcvdp; /*the phase of the received vectort 

whose components are ik and qk*/ 
double ihatk; /*for this phase detector this is the 

output from the circular 
quantizer*/ 

double qhatk; /*for this phase detector this is the 

output from the circular 
quantizer*/ 

double errkhold; 
int count; 

/*declare noise model parameters, sigma is the WGN standard 
deviation, nik and nqk are the gaussian noise samples at time 
k with mean zero and standard deviation sigma.*/ 

double sigma; 
long xconl; 
long xcon2; 
long x; 
long y; 

double samplel; 
double sample2; 
double nik; 
double nqk; 

/*declare output file parameters*/ 

FILE *outfile; 
char *modehold; 
char *nameholdl; 

/initialize file output parameters*/ 

modehold = "w"; 

nameholdl = "kghigh"; 

outfile = fopen( nameholdl, modehold); 

/initialize uniform RV seeds for repeated tests*/ 

srand(2546); 
y = 47; /*seed*/ 

/initialize phase step*/ 

phsfix = 0.0; 
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/*input natural frequency, damping factor for loop*/ 

printf("This is a second order carrier loop synchronizer 
simulator\n"); 

printf( "Please enter the natural frequency of the loop 

transfer functionin' 1 ); 

printf("\n"); 

scanf("%s",dptr); 

wn = strtod( dptr ,NULL) ; 

printf( "Now please enter the loop damping factorin ), 

printf("\n"); 

scanf("%s",dptr); 

printf("\n"); 

damp = strtod( dptr ,NULL) ; 

printfl "And finally, please enter the phase detector gamin ), 

printf("\n"); 

scanf ( "%s" .pdgains) ; 

printf("\n"); , 

pdgain = strtod(pdgains,NULL), 

printf("\n"); 

printf("\n"); 

printf("wn = %f, damp = %f\n",wn,damp); 
printf( "phase detector gain = %f\n",pdgain); 
bw = wn*(damp + 1.0/(4.0*damp)); 
kO = 2*damp*wn/ pdgain; 
a = wn/(2*damp); 

/*determine sample interval that is required to uncorrelate 
phase error variance calculation / 


uncore = floor(3.1415927/wn/4.0); 


printf( "Tim loop noise bandwidth is: %13.10f Hz\n ,bw), 

printf("The gain is; %f rad/sec/Volt\n ,k0), 

printf("The^loop filter time constant is: %f sec\n ,a), 

printfl "The* number of skipped samples for correlation is. %d 

\n", uncore); 

printf("\n"); 


217 


fprintf(outfile,"8-PSK High SNR MAP estimator loop 

results\n"); 

fprintf(outfile,"\n"); 

fprintf(outfile, "Initial phase step = %f\n",phsfix); 
fprintf (outfile, "\n" ) ; 

fprintf(outfile, "Noise random seed = %ld\n",y); 
fprintf ( outfile ," \n ") ; 
fprintf(outfile,"wn = %f\n",wn); 
fprintf (outfile, " \n " ) ; 

fprintf(outfile, "damping factor = %f\n",damp); 
fprintf(outfile,"\n"); 

fprintf (outfile, "The loop noise bandwidth is: %f Hz\n",bw); 
fprintf(outfile,"\n"); 

fp rintf ( outfile , "The gain is: %f rad/sec/Volt\n",kO); 
fp rintf ( outfile , "\n ") ; 

fprintf(outfile,"The loop filter time constant is: %f sec\n",a); 
fprintf(outfile,"\n"); 

fprintf(outfile, "phase detector gain = %f\n",pdgain); 
fprintf(outfile,"\n"); 

fprintf(outfile,"The number of skipped samples for 
correlation is: %d \n", uncore); 
fprintf(outfile,"\n"); 

/*begin simulation*/ 


for( j=30;j<=50;j=j+5) 


endb = j; 

sigma = 1.0/sqrt(2.0*pow(10.0,endb/10.0)); 

/*initialize system parameters*/ 

ilk = 0.0; 
i2k = phsfix; 
refkl = i2k; 
t = 1.0; 

errkhold = 0.0; 
phsk = 0.0; 
errk = 0.0; 

xcon2 = 16807; /*for RV generation*/ 
xconl = 2147483647; /*for RV generation*/ 
sqrerr = 0.0; 
sumerr = 0.0; 
s = 0; 

for(k=0;k<=120;++k) 

I 

/Construct received phase, noise, and transmitted data*/ 
/*received phase*/ 
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/*running variance calculation*/ 


if (k >20) 

if (fabs(refkl) < 0.392699) 

{ 

s =s + 1; 

sqrerr = sqrerr + refkl*refkl; 
sumerr = sumerr + refkl; 

1 


for(n=0;n<uncore;++n) 

{ 


/* phsk = 0.001*k;*/ 

/*noise*/ 

/* x = rand(); 

y = rand(); 

sample 1 = (x/32767.0)*0.999999; 
sample2 = (y/32767.0)*0.999999; 

*/ 

x =ldiv(y,xconl).rem; 
x=labs(x*xcon2); 

y =ldiv(x,xconl).rem; 
y=labs(y*xcon2); 

sample 1 = (x/2147483647.0)*0.999999; 
sample2 = (y/2147483647.0)*0.999999; 

nik = sigma*sqrt(-2.0*log(1.0-samplel)) 
*cos(2.0*3.14159265*sample2); 

nqk = sigma*sqrt(-2.0*log(1.0-samplel)) 
*sin(2.0*3.14159265*sample2); 

/*transmitted data*/ 

d = rand(); /*d is now divided into 8 regions each 

representing a possible xmitted symbol*/ 


if(d <4096) 

\ 

imk = 0.923879; 
qmk = 0.382683; 

else if((d >= 4096) && (d < 8192)) 

{ 
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imk = 0.382683; 
qmk = 0.923879; 

I 

else if((d >= 8192) && (d < 12288)) 

{ 

imk = -0.382683; 
qmk = 0.923879; 

else if((d >= 12288) && (d < 16384)) 

{ 

imk = -0.923879; 
qmk = 0.382683; 

} 

else if((d >= 16384) && (d < 20480)) 

I 

imk = -0.923879; 
qmk = -0.382683; 

1 

else if((d >= 20480) && (d < 24576)) 

{ 

imk = -0.382683; 
qmk = -0.923879; 

I 

else if ((d >= 24576) && (d < 28672)) 

i 

imk = 0.382683; 
qmk = -0.923879; 


else 

{ 

imk = 0.923879; 
qmk = -0.382683; 
I 


/*end generation part of the simulation, this generated 
information will be used for this value of the simulation 
counter only*/ 

/*subtract the reference phase from the received phase*/ 
pdifk = phsk - refkl; 

/Calculate the received symbol values for I and Q. based on 
the generated data and noise and the phase difference pdifk*/ 

ik = imk*cos(pdifk) - qmk*sin (pdifk) + nik; 
qk = imk*sin(pdifk) + qmk*cos(pdifk) + nqk; 

/*Now implement a phase detector to generate a phase error 
signal errk. This implementation performs a circular 
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quantization on I and Q. followed by some algebra to develop 
the error signal.*/ 

/Calculate the phase of the received vector. The atan2 
function adjusts for the appropriate quadrant.*/ 

rcvdp = atan2(qk,ik); 

/*quantize rcvdp into ihatk and qhatk / 

if ((rcvdp >= 0.0) && (rcvdp < 0.785398)) 

/* 0 <= rcvdp < Pi/4 */ 

{ 

ihatk = 0.923879; 
qhatk = 0.382683; 

else if ((rcvdp >= 0.785398) && (rcvdp < 1.570796)) 

/* Pi/4 <= rcvdp < Pi/ 2 */ 

{ 

ihatk = 0.382683; 
qhatk = 0.923879; 

else if /(rcvdp >= 1.570796) && (rcvdp < 2.356194)) 

/* Pi/2 <= rcvdp < 3Pi/ 4 */ 

I 

ihatk = -0.382683; 
qhatk = 0.923879; 

else if /(rcvdp >= 2.356194) && (rcvdp < 3.141593)) 
/* 3Pi/4 <= rcvdp < Pi */ 

{ 

ihatk = -0.923879; 
qhatk = 0.382683; 

else if /(rcvdp < -2.356194) && (rcvdp >= -3.141593)) 
/* -3 Pi/4 < rcvdp >= -Pi V 

ihatk = -0.923879; 
qhatk = -0.382683; 

} 
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else if ((rcvdp >= -2.356194) && (rcvdp < -1.570796)) 
/* -Pi/2 > rcvdp >= -3Pi/4 */ 

I 

ihatk = -0.382683; 
qhatk = -0.923879; 

1 

else if ((rcvdp >= -1.570796) && (rcvdp < -0.785398)) 
/* -Pi/4 > rcvdp >= -Pi/ 2 */ 


1 

ihatk = 0.382683; 
qhatk = -0.923879; 

1 

else if ((rcvdp >= -0.785398) && (rcvdp < 0.0)) 
/* 0 > rcvdp >= -Pi/4 */ 


ihatk = 0.923879; 
qhatk = -0.382683; 

I 

/^calculate the phase error errk as a function of ik, qk, ihatk, 
and qhatk.*/ 

errk = -ik*qhatk + qk*ihatk; 

/*This is the end of the phase detector*/ 

/*The next step considers the loop filter and VCO in one 
transfer function, as a difference equation. It calculates the 
phase of the VCO, refkl, for the next simulation iteration, 
refkl is based not only on the error signal, errk, but also on 
the state of the integrators that make up the one pole low-pass 
filter and the VCO. The filter can be represented by the 
function, 1 + a/s. a is one-over-the time constant, RC, of the 
filter. The VCO xfer function can be considered as kO/s. kO is 
the ratio describing the output frequency of the VCO in terms 
of the input error voltage from the loop filter. The filter 
integrator initial condition at time k is ilk. The VCO integrator 
initial condition at time k is i2k. The integration period is t. In 
discrete time these transfer funcitons can be implemented as 
follows*/ 

refkl = i2k + k0*t*errk + ilk*k0*t + a*k0*t*t*errk/2.0; 

/*now reset the integrator intial conditions for the next 
iteration*/ 
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i2k = refkl; 

ilkl = ilk + a*t*errk; 

ilk = ilkl; 

} 

/*printf(" k = %d\n",k);*/ 

} 

meanphs = (sumerr/s); 

varphs = sqrerr/ s - meanphs*meanphs; 

f printf ( outf ile , "\n " ) ; 

fprintf(outfile,"E(s)/N(o) = %4.1f\n",endb); 

fprintf(outfile,"Mean of error = %10.7f, variance of error = 
% 1 3 . 10 f\n", meanphs, varphs) ; 

fprintf(outfile,"\n"); 

fprintf(outfile,"The number of samples used in the previous 
statistical analysis is: %d\n",s); 

fprintf(outfile,"\n"); 

printf("\n"); 

printf("Es/No = %4.1f, mean = %f, variance = %f, samples = 
%d\n",endb, meanphs, varphs, s); 

i 

fclose(outfile); 

} 


Both simulations were run with a loop bandwidth-to-symbol rate 
ratio of 0.53% and damping factor of 0.7071. The data plotted, in 
dB, in Table A-l indicates the high-SNR loop performs better than 
the Leclert & variation in terms of variance of phase error. 
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Table A-l. The simulation data for comparing loop designs. 



Leclert & 

1 

SNR 

MESEEEBi 


K Jnm 



24 

-44.271 

-49.002 

22 


-48.365 

20 

-41.244 

-45.642 

18 

-39.642 

-42.529 

16 

-37.535 

-40.471 

14 

-34.825 

-39.095 

12 

-30.965 

-32.0768 

1 0 

-25.655 

-26.655 
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Appendix B 


THE 8PSK MAP ESTIMATOR DERIVATION 


Extending the MAP estimator result for 8PSK (M=8) is a 
simple process. Equations (3-1) through (3-13) are applicable for 
any M. The modulation probability density in (3-14) is different 
for 8PSK and is stated as 
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We proceed by inserting (3-12) into (B-2) for all values of 8 m . To 
do this the two quantities r x -s x and r 2 -s 2 must first be analyzed 
for the 8 values of 0 m . Table B-l contains this analysis. 


Table B-l. Two dimensional noise component analysis. 
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Now using Table B-l it is possible to insert (3-12) into (B-2) 
and obtain the likelihood function. The process of obtaining the 8 
PSK likelihood function will only be conducted for the first term of 
(B-2). The procedure is the same for all the terms in (B-2). 
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Inserting the data from Table B-l for the first modulation phase, 

namely when 6 m = the first term in (B-2) can be stated as 
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Completing the square of the exponential arguements in (B-3) and 
then combining, results in 
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By combining and separating like terms (B-4) becomes 
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Applying the appropriate trigonometric identities to expand the 
sinusoidal terms forms 
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Inserting the equations for r x and r 2 , given in (3-9) resolves (B-6) 
into 


(A) = 


1 


16;rcr 2 


r 2 +r\+E 

~ r ~ 


lr(t)J-cos(a> c t)dt 


V2£ 


cos(e, )cos^| j-sin(0 f )sin^j 


xe 


xe 


2a z 


(B-7) 


I r U):M sin(co c r)dr 

0 ' 1 


~J2E 


sin(e, )cos(|) + cos( 01 )sin^jj 


2 a 


By combining the integrals and distributing the carrier 
components within the integrand (B-7) becomes 
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Next the trigonometric identities 


cos(x)cos(y) - sin(x)sin(y) = cos{x + y) 


and 


(B-9) 


cos(x)sin(y) + sin(x)cos(y) = sin(x + y) 


are used to separate the integral in (B-8) into two new integrals. 
This new form is 
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reduces (B-10) to 
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Repeating the process for each of the seven remaining terms in 
the 8PSK likelihood function is accomplished by changing the 
trigonometric coefficients in (B-10) to accomodate the modulation 
angle corresponding to each term. Therefore the second term 
becomes 
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<C> = 


16 no 


(D). 


16 kg 1 


(E) = 


\6no 2 


rj+rj + E 



^ COS; 

2ct 2 e 

l 8 J 

r, 2 +r 2 +£ 

( In > 


^ COS 

2ct‘ e 

( 8 j 

r 2 +r 2 +E i 

' -* y 


<■» cos 

2a 2 e l 

, 8 J 


(B-14) 


(B-15) 


(B-16) 


230 



(F) = 


16 KO 


r*+r\+E (^n\ 


la 


cost — jx 




r^+ri+E (-5n 

1 " J — T — cos "7 

(G) = — - — s-e 211 e ( 8 




\6no 


r 2 +r 2 +£ (- ln\ . (~Tn\ 

i r^A— M—>. 


(H> = 


\6nc 


(B-17) 

(B-18) 

(B-19) 


All of the trigonometric coefficients in(C) through (H) can be 
written in terms of cosjjj, sin(jj, cos(7-)> and sm (irJ' 

Performing the necessary conversions yields 
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the 8PSK likelihood function can be expressed as 
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Expanding the exponents and rearranging results in 
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This reduces to 
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Using the definition of the hyperbolic cosine function, the 8PSK 
likelihood function becomes 
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Next the log likelihood function is obtained and can be stated 
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\6kc 


r 2 +rf+£ 

2a 1 


(B-31) 


— 

ln[cosh(c 1 Jc)cosh(i , 1 >') + cosh (c 2 x ) cosh ( ^ 2 v )] 


Taking the derivative of (B-31) with respect to 6, and noting that 
the first term in the log likelihood function will equate to zero in 
the derivative, the MAP estimator equation is found. The first step 
in computing the derivative is to apply the chain rule yielding 
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(B-32) 


<9[cosh(c! ^)cosh(j 1 y) + cosh(c 2 ;t)cosh(s 2 ;y)] 

= <90, 

d&i cosh(c ] x)cosh(^ 1 }') + cosh(c 2 J:)cosh(j 2 >') 


Now taking the derivative of the products results in 


<9{ln[/?(/g|0,)]} 

<90, 


+ COS! 


h(^)^ c °sh(c,, ) ] + coshM ib^ 
hM a[c^^ +co ^ (c2i) £[co^)l 


cosh(c 1 Jc)cosh(5 1 j) + cosh(c 2 Jt)cosh(j 2 j) 


(B-33) 


Applying the chain rule to each of the derivatives on the right 
side of (B-33) forms 


<90, 


Cj cosh(j 1 y)sinh(c 1 x)-^[^ 

OU l : 

+^! cosh(c ] j:)sinh(s|3')^^ 

OUl 

+c 2 cosh(s 2 ;y)sinh(c 2 x)^l 


+^ 2 cosh(c 2 x)sinh(j 2 y)^^ 

OUi 


cosh(c 1 x)cosh(5 1 >’) + cosh(c 2 *)cosh(.s 2 ;y) 


(B-34) 


Using (B-ll) it is noted that 


Substituting (B-35) into (B-34) yields 




In 


ddj 


-c 1 ycosh(s 1 3 ')sinh(c,jc) + s,xcosh(c 1 x)sinh(s 1 y) 1 
‘ -c 2 y cosh(j 2 y)sinh(c 2 Jc) + j 2 scosh(c 2 x)sinh(j 2 y) j 
cosh(c 1 x)cosh(ji>') + cosh(c 2 x)cosh(.s 2 y) 


(B-36) 


Writing (B-36) as a sum of two fractions yields 


d { 


In 


d6: 


-Ciycosh(5 1 y)sinh(cix) + 5,xcosh(c 1 x)sinh(^ 1 y) 
cosh(c,x)cosh(s 1 y) + cosh(c 2 x)cosh(.s 2 y) 
-c 2 ycosh(s 2 y)sinh(c 2 x) + j 7 scosh(c 2 s)sinh(s 2 y) 
cosh(c 1 x)cosh(5,y) + cosh(c 2 Jc)cosh(j 2 y) 



Multiplying the first fraction by p osh ( c i ) cosh (-yi y) and the second 

cosh(cix)cosh(j 1 >'J 

fraction by cosh M^M forms the sum 

cosn(c 2 Jt)cosh(,s 2 ;y) 


dO t 


+ 


-cjy cosh(5 1 y)sinh(c 1 .x) + ^xcoshfajtjsinh^ y) 
cosh(cj x ) cosh( jjy) 

cosh(c 1 x)cosh(^,y) + cosh(c 2 .r)cosh(j 2 y) 
cosh(c 1 Jc)cosh(j 1 j) 


c 2 y cosh(.y 2 y)sinh(c 2 Jc) + s 2 ;ccosh(c 2 ;c)sinh (,y 2 y) 1 * 
cosh(c 2 *)cosh(.s 2 y) 

cosh(ci^)cosh(j 1 }') + cosh(c 2 *)cosh(.s 2 ;y) 
cosh(c 2 Jc)cosh(s 2 y) 


(B-38) 


This reduces to 


^{ ln [p(^0,)]} 

do. 


.y 1 xtanh(.y 1 >')-c 1 ytanh(c 1 A:) s 2 x tanh(s 2 y) - c 2 y tanh(c 2 x) 
j + cosh(c 2 *)cosh(s 2 y) cosh(c 2 Jt)cosh(s 2 y) 

cosh(c 1 x)cosh(>y ] y) cosh(c 2 x)cosh(j' 2 y) 


(B-39) 


Using the same approach for completing the MAP estimator as 
was done for BPSK and QPSK, the right side of (B-39) is equated to 
zero yielding 
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s ,xtanh( S ,y) - c,y tunhM s,xtanh( J? y)-c 2 ytanh(c 2 *) = 0 . (B-40) 

cos h(c 2 jc)cosh^ 2 y) ~ t + _cosh(c 1 x)cosh(.y 1 y£ 
1+ cosh(c,x)coshM cosh(c 2 x)cosh(s 2 y) 
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Appendix C 


THE 16PSK MAP ESTIMATOR DERIVATION 


Extending the MAP estimator result for 16PSK (M=16) is a 
simple process. Equations (3-1) through (3-13) are applicable for 
any M. The modulation probability density in (3-14) is different 
for 16PSK. Writing p(d m ) for 16PSK as a summation while 

preserving the positive and negative references for the 
modulation points, as was done for the previous PSK modulations 
formats, results in 


P{ 6 m) = I J£ S \ 


n =- 8 


e„ - 

m 


(2 n + l);r 
16 , 


(C-l) 


The 16PSK likelihood function can be expressed in the same 
manner, namely, 


'(«= „?>l 


R 


Oi’dm = 


(2» + l);r 
16 


(C-2) 


We proceed by inserting (3-12) into (C-2) for all values of d m . To 
do this the two quantities r, -s, and r 2 -s 2 must first be analyzed 
for the 16 values of 6 m . Table C-l contains this analysis. 
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Table C-l. Two dimensional noise component analysis. 


ri~s i 





r x - Vficos|^0 z - + — 


/ 3 ^ 

r\ - Vs cosl 0, + — 


rj - Vs cosf 0, + — 


rj - VScos^0,- + ■ 



r- ( n 

r, - Vs cosl 0, + -j-g- 


f 13 tt 
r! - Vs cosl 0, + -^g- 





rj - cos^ 


r x - VS cos^ 


Tj — Vs COS^ 


r, - Vs cos^ 


/ 

rj - Vs cos 


n JL 
6i 16 


01 16 , 


^ 16 


e _Z£ 
0 ‘ 16 


0 — 
'16 


/— („ 1 ITT 

rj -VS cos| 0j --jg- 


/-=■ (. 137T 

rj -VScosj 0,---^g- 


15tt 

rj - VS cosj 0, --gg- 


^ 2-^2 












r 2 + Vs sin^0j + yg 


r 2 + VSsin^0 ( +^g- 


f 57 T 
r 2 + VS sin! 0, + — 


r 2 + Vs sinl 0, + — 


n= . ( _ 97T 

r 2 +VS sin 0.+ — 
z 16 


r 2 + VS sin^0j + gg^- 


y— . f _ 137T 

r 2 + VS sinl 0, + - 


r= . ( n 15 K 
r 2 + VSsinl 0, +-jg- 


r 2 + Vs sin^0,- — — 


r 2 + VS sin^ 0,-^1 


/ 5 ^' 

r 2 + VS sinl 0, - — 


r 2 +VSsin^0,-yg 


/ 9 ^ 

r 2 + VSsinl 0, - — 


r 2 +VSsin^ 0 y-gg^ 


/- • 13/r 

r 2 + VSsml 0, - — 


r- • 15 /r 

r 2 + VS sinl 0, — — 
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Now using Table C-l it is possible to insert (3-12) into (C-2) and 
obtain the likelihood function. The process of obtaining the 16 PSK 
likelihood function will only be conducted for one term of (C-2). 
The procedure is the same for all the terms in (C-2) and the 
results will be presented after the derivation for this first term. 
Inserting the data from Table C-l for the modulation phase 

corresponding to 6 m = — — , the related term in (C-2) can be stated 

16 

as 


< a > 4 * 


R 





- 2 
[r,-V£cos(0 i+ ^]] 

^ 1 

e 2 ° 2 

) 32 KG 2 

2 

(r 2+ V£sin(0 I+ ^]] 


xe 2ff2 


(C-3) 


Completing the square of the exponential arguments in (C-3) and 
then combining, results in 


r\ -2rj 41 cos^ 0, j+£cos 2 ^0, +yl j 
+r\ +2 r 2 V£ sin^ 0, j+Esin 2 ^ 0, j 


<A) = 


2 a 1 


32 KG 


(C-4) 


By combining like terms (C-4) becomes 


<A> = 


l 

32 no 1 


r\ +r\+E-2r x 4E cos^ 0, j+2 r 2 VI sin^ 0, 


2 a 


(C-5) 
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Applying appropriate trigonometric identities to expand the 
sinusoidal terms forms 


(A) = 


1 


r 2 +r 2 +E 

2 a 2 


32k tr 

-2rjV£ 

xe 

2r 2 4E 

xe 


cos(0, )cos(-j| j-sin(0, )sin^ 


2a L 


sin(0, )cos^ j+cos(0, )sin^ j 


2a 


(C-6) 


Inserting the equations for r, and r 2 , given in (3-9) resolves (C-6) 
into 


(A) = 


1 


32 kg 2 


r 2 +r 2 +E 
'~2a 2 ~ 


T [2 

J'-(Oi-cos(£o c r)* 


V2£^cos(0, )c°s[^]-sin(0, )sin^ j j 


X e 


xe 


2er 


(C-7) 


M f )-^r sin (<»cO < * sin( 0, ) cos^ ^ j+c° s ( 0 i ) sin )] 


2 a 


By combining the integrals and distributing the carrier 
components within the integrand (C-7) becomes 
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2a 2 


(A> = 


1 

32 kg 2 


[pT 

c°s( w c /) cos( 0,- ) cosf-^- Vcos(<a c /)sin(0,- )sinf 3i) 



VioJ U6 J 

dt 


-sin(to c f)sin(0 f )cos|^— j-sin(o) c /)cos(^- )sin^-^ j 



x e 2 a 2 


(C-8) 


Next the trigonometric identities 


cos(x)cos( y) - sin(.x)sin(;y) = cos(x + y ) 


and 


(C-9) 


cos(*)sin(;y) + sin(jr)cos(j) = sin(jc + y) 

are used to separate the integral in (C-8) into two new integrals. 
This new form is 


<A> = 


1 

T e 

32 kg 2 


2 7 

r, +r| +E 

2 cos 

l r(t ^ COS ( W c t+0 i) c/t 


-2 sin 

(~ 

U6, 

2<j 2 

jjj: 1 r(t)sin(co c t+d i )dt 

2d 2 e \ 

2a 2 J 


(C-10) 


Using the substitution adopted for QPSK, namely 



X = 




~e t 

- 

To 


r(t)cos(a> c t + di)dt 


and 


(C-ll) 


y- 


-^2 Jr(Osin(o c / + 
(TUo 




reduces (C-10) to 


r?+r?+£ 


(A) = 


1 

32 no 2 


2 a 


- cosl 

e ^ 




(C-12) 


Repeating the process for each of the fifteen remaining terms m 
the likelihood function is accomplished by changing the 
trigonometric coefficients in (C-10) to correspond to the 
modulation angle for that term. Therefore the second term 

becomes 


(B) = 


32 no 


rj+rl+E 

la 1 e 




(C-13) 


Similarly the other fourteen terms are 


r\+r\+E ( 

A k rnsl 


(0 = 


2(7 


.2 C0S l 




(C-14) 


32 no 
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r\ l+r^+E ( ln \ fj n \ 

1 -?> — cos — Lt-sin — y 

i * 2 a 2 P v 16 J 1 16/ 


(D) = 


32 kg 


(C-15) 


r \ +r 2 + E (9n\ . (9n'\ 

l — ^2— cos hrH ,n hr 


(E) = 


32 kg 


2a l P V 16 y v 16 


(C-16) 


r [ +> j+e nijr'i . niw> 


(F) = 


1 — = — =* — COS 

,, 2a 2 „ v 16 


32 KG 


\ ■ (n*\ 

■ x-sin v 

J 16 / 


(C-17) 


(G) = 


rf+rj+E 

1 — ! H — cos pr-sin y 

1 * 2a 2 l 16 7 l 16 V 


32kg 1 


(C-18) 


(H) = 


r\+rl+E 715 ^ 715 ^ 

i — S — * cos bt-sin y 

1 2c 2 p V 16 7 l 16 J 


32 kg 


(C-19) 


(/>= 


32 kg 1 


rf+rj+E 
2 O 2 A 


( Jt\ . ( 7 t\ 

cos x-sm j 

V 16 J l 16/ 


(C-20) 


rJ+rj+E ( 2n \ 7 37 r> 


(J) = 


i — — cos x-sin y 

1 2cr 2 . I 16 J { 16/ 


32 KG 


( 3n\ . ( 

s *-sin - 

l 16 J l 


(C-21) 


r 2 +r|+£ f 5J) . 7 5/r\ 

— 1 — ^ — cosl x-sin > 

i J y i6 / 


(K) = 


32 kg 1 


2o l p \ 16 


(C-22) 


{L) = 


rf+rj+E 7 ln 
j — 5 ^ cosf 


2 a , l 16 


W-T?> 


32 kg 


(C-23) 
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r}+r\+E ( lljr'l . ( 1 l/rA 

i —73- Mir r r 16 7. 


(N) = 


32kg 1 


-e 2o ' e 


r}+r\+E ( 13 ?^ . f 

1 


( 0 ) = 


32 no 2 




^ 32 no 2 € 


+rj+E ( 15 *^ . ( 15 n\ 


2a z 


up 


(C-24) 

(C-25) 

(C-26) 

(C-27) 


All of the trigonometric coefficients in(C) through (P) can be 


written in terms of cos (yg > sin (^)’ cos (kt)’ and sm ( 


3k 

16 


Performing the necessary conversions yields 


rl+rl+E . ( 3 n\ ( 3 ^ 

1 —^Aur^itr. 


(C) = 


32 KO 


(D) = 


r?+r?+£ 

1 s "" 

^ e la e 




32 KO 1 


(E) = 


32no 


rf+rj+E 

2 a 2 




r l 2+r 2 +£ • ( 

1 2 Sin 

^=^? e 20 £ 


3;r^ (3*^ 

— x-cos — > 

16 J U 6 ) 


(C-28) 

(C-29) 

(C-30) 

(C-31) 
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1 


e 


(C-41) 


</>) = 


r 2 +r 2 
'“^ r_ 


-cos(^]x + sin(^) y 


32 kg 


Now making the assignments 


K . K 

c,=cos-, s ,=sm- 


and 


(C-42) 


3 K . 3 K 

c 2 =c° s — , s 2 = sin — , 

the 16PSK likelihood function can be expressed as 

e c\*-s\y + e c 2 x-j 2 y + e s 2 x ~ c 2y + e s \ x ~ c \y 
+e ~ s i x ~ c \y + e ~ s 2 x ~ c 2y + e ~ c 2 x ~ s 2y + e - c \*-*\y 

+e c\X+s x y + g C2 x +S2y + e n x + c 2y + e s \ x+c \y 
+e ~six+ciy + e -s 2 x+c 2 y + e ~c 2 x+S2y + e ~ c l x+s iy 

(C-43) 

Expanding the exponents and rearranging, yields 

e c \ x e s \y + e ClX e~ s ' y + e~ Ct *e Siy + e~ ClX e~* iy 
+e C2X e S2y + e C2X e~ S2y + e~ C2X e Siy + e~ C2X e~ Siy 
+e s i x e c ' y + e s ' x e~ c ' y + <T 51 V 1 * + e -' 1 V c,y 
+e S2X e C2y + e S2X e~ C2y + e~ $2X e C2y + 


rf+rf+E 


pW»,)= 


32 kg 


ic r 


r 2 +r 2 +E 


p(R|e,)= 


32na 


2a 
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(C-44) 


This reduces to 


p(«|0,) = 


r\+r\+E 


2cr z 


32 K( 7 Z 


\e c ' x + e~ c ' x )[e s ' y +e~ s ' y ) 
+(e C2 x + e~ C 2 X )(e S2 y + e~ S2y ) 
+(e s ' x +e~ s ' x )(e c ' y + e- c ' y ) 
\+[e* 2X +e~ S 2 X )[e C2y + e~ C2y ) 


(C-45) 


Using the definition of the hyperbolic cosine function, the 16PSK 
likelihood function becomes 


p(p|e f ) = 


r 2 +r]+E 

L- e 2^~ 

32 kg 2 


cosh(c 1 ^)cosh(5 1 >') + cosh(c2x)cosh(^2> ; ) 
+cosh(5jx)cosh(c!_y) + cosh^Jcjcosh^}') 


(C-46) 


Next the log likelihood function is obtained and can be stated 


ln[p(tf|0 ( .)] = ln 


22 kg 


rf+rj+E 

2a 2 


+ In 


cosh(q x) cosh(sj y) 
+cosh(c2Jc)cosh(5 2 >') 
+ cosh(sj jc) cosh (q y ) 
+cosh(5 2 ^)cosh(c 2 y) 


(C-47) 
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Taking the derivative of (047) with respect to 8, and noting that 
the first term in the log likelihood function will equate to zero in 
the derivative, the MAP estimator equation is found. The first step 
in computing the derivative is to apply the chain rule yielding 


a{ln[p(R|^)]} 

36 , 


COS 
+ COS 


h(c 1 x)cosh(s 1 ;y) + cosh(c 2 *)cosh(s 2 ;y) 
h(^ 1 A:)cosh(c 1 y) + cosh(.S 2 *) cos h( c 2 ) ) ), 


ddi 


cos h(cj^)cosh(5iy) + cosh(c2*)cosh(s2y) 
+cosh(s 1 x)cosh(c,y) + cosh(s 2 ^)cosh(c 2 y) 


. (048) 


Now using the rule for the derivative of a product results in 

. dfcosh(ci*)l , , x d[cosh(sjy)] 
jsh M - 1 + cosh (c,x) ^ 

^[cOSh(coJc)l . a[cOSh(5 2 y)] 

coshM- 1 -^-^ ^ + cosh M — 




In 


do, 


dfcoshfsi*)] / \ <9[cosh(c 1 y)] 

+ cosh(q y) — *■ cosh(sj x) ^ 


+cos 


. dfcosh(^x)l ,a[cosh(c 2 y)] 

h M- L a g. - JA+ cosh M — 


COS h(c 1 x)cosh(s,y) + cosh(c 2 x)cosh(s 2 ;y) 

+ COS h(5!x)cosh(c,y) + cosh(s 2 *)cosh(c 2 y) 


(049) 


Applying the chain rule to each of the derivatives on the right 
side of (049) forms 
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afln[p(«|8,)]} 
dO : 


q cosh(^ 1> y)sinh(c 1 x)^^ 


dO i 

+Sj cosh(c 1 x)sinh( 1 s , 1 j)^|^ 

dOi 

+c 2 cosh(s 2 ;y)sinh(c 2 x)^j^ 

+s 2 cosh(c 2 *)sinh(j 2 3>)^[^ 

dQ t 

+5, cosh(c 1 >')sinh(5 1 J:)^[^ 

dO { 

+q cosh(^x)sinh(c,y)^^ 

dOi 

+s 2 cosh(c 2 >')sinh(s 2 j:)^^ 
+c 2 cosh(5 2 j:)sinh(c 2 y)^lj 


cosh(c 1 A:)cosh( 1 Si>') + cosh(c 2 .x)cosh(s 2 )') 
+cosh(^x)cosh(c, j) + cosh(s 2 *)cosh(c 2 ;y) 


Using (C-ll) it is noted that 


M = _, 

dOi y 


and 


<90, 


= x. 


Substituting (C-51) into (C-50) yields 


(C-50) 


(C-51) 
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d{ln[p(fl|e,-)]} 

do t 


-c 1 >'Cosh(^ 1 }')sinh(c 1 x) + 5 1 Jccosh(c,x)sinh(5 1 >) 
-c 2 ycosh(5 2 }')sinh(c2J:) + S2 j:cosh ( c 2 ;c ) sinh ( s 23') 

-5 1 ycosh(c 1 y)sinh(5 1 Jc) + c 1 xcosh(j 1 x)sinh(c 1 >') 

-s 2 y cosh(c 2 y)sinh(j 2 x) + c 2 xcosh(s 2 x)sinh(c 2 y) 
cosh(c 1 x)cosh(^ 1 y) + cosh(c 2 x ) cosh(s 2 y ) 

+ cosh^x) cosh(cjy) + cosh(j 2 Jc)cosh(c 2 >') 


(C-52) 


Writing (C-52) as a sum of four fractions yields 




In 


-c, 3 ;cosh(^y)sinh(c 1 j:) + SiSCOshfajQsinh^y) 


36 ; 


cos h(c 1 x)cosh(^ 1 )') + cosh(c 2 x)cosh(s 2 y) 

+ cosh( x) cosh(c 1 >') + cosh(.s 2 x)cosh(c 2 y) 
-c 2 y cosh(.s 2 y)sinh(c 2 Jt) + .s 2 -x:cosh(c 2 .x)sinh(.s' 2 y) 


+ 


cosh(c 1 jr)cosh(j 1 y) + cosh(c 2 .x)cosh(s 2 y) 
+cosh(s 1 ;t)cosh(c 1 y) + cosh(.s 2 .x)cosh(c 2 y) 
-s x y cosh(c 1 y)sinh(^Jc) + qxcoshfa x)sinh(ciy) 


cosh(c 1 x)cosh(5 1 y) + cosh(c 2 x)cosh(s 2 y) 
+cosh(s,;t)cosh(c,y) + cosh(s 2 x)cosh(c 2 y)J 
-s ? ycosh(c 2 y)sinh(:? 2 :t) + c 2 scosh(s 2 x)sinh(c 2 y) 


cosh(cjx)cosh(^y) + cosh(c 2 x)cosh(.s 2 y) 
+cosh(s 1 ;t)cosh(c 1 y) + cosh(s 2 x)cosh(c 2 y)J 


(C-53) 
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Multiplying the first fraction by t he second 

cosh^qxjcosh^yj 

„ , coshfcoJtlcoshf^y) A A . j , u coshi^, x ) cosh(q y ) 

fraction by ^ the third fraction by — — tt-K. 

cosh(c 2 ;c)cosh(.y 2 y) cosh(5] x J cosh(c 1 y ) 

and the fourth fraction by cos M* 2 *)cosh(c 2 yj f orms the reduced sum 

cosh[S2X)cosn{c 2 y) 


djl n [p(/?|0 ( )]j _ -c K ytanh(c 1 -y) + 5 1 Jctanh(.y | y) 

d6 t f + cosh(c 2 x)cosh(^ 2 y) + 

cosh(qA:)cosh(^jy) 

coshf^-xjcoshfay) | cosh(s 2 j:)cosh(c 2 y) 
cosh(c 1 Jc)cosh(5 1 y) cosh(c 1 ^)cosh(j ] y) 

-c 2 y tanh( c 2 x) + s 2 x tanh(.s 2 y) 

cosh(cjJc)cosh(5!y) + 
cosh(c 2 ;c)cosh(.s 2 y) 

cosh(5 1 j:)cosh(c 1 y) cosh(s 2 .x)cosh(<: 2 y) 
[ cosh(c 2 x)cosh(y 2 y) cosh(c 2 x)cosh(s 2 y) 

— jjy tanh ^x) + c x x tanh(c 1 y) 

cosh(c 1 x)cosh(5 1 y) cosh(c 2 j) cosh(^ 2 y) 
cosh(5,x)cosh(c ] y) cosh(y 1 ^)cosh(c,y) 

cosh(s 2 ;t)cosh(c 2 y) 
cosh(^ x) cosh(cj y) 

-s 2 y tanh(i 2 x) + c 2 x tanh(c 2 y) 


+ -F 


cosh(c 1 x)cosh(5 1 y) cosh(c 2 Jc)cosh(s 2 y) 
cosh(^ 2 x)cosh(c 2 y) cosh(s 2 *)cosh(c 2 y) 

cosh(s|^)cosh(cjy) 

cosh(s 2 ;t)cosh(c 2 y) 


(C-54) 
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Using the same approach for completing the MAP estimator as 
was done before for BPSK, QPSK, and 8PSK, the right side of (C-54) 
is equated to zero 


-qytanhfax) + ^xtanhfoy) 


r 

coshfc-ixlcoshl v) cosh(jix)cosh(c 1 y) cosh(s 2 x)cosh(c 2 y) 


1 + cosh(c 1 x)cosh(5 1 y) cosh^^cosh^y) cosh(c,x)cosh(^>-) J 



-c 2 y tanh(c 2 x) + s 2 x tanh(s 2 y) 


+ - 

'cosh(c,x)cosh(5 1 y) l cosh{s { x)cosh(c x y) ^ cosh{s 2 x)cosh{c 2 y) 


cosh(c 2 x)cosh(s 2 }’) cosh(c 2 x)cosh(s 2 y) cosh(c 2 x)cosh(.s 2 y) _ 


-s, y tanh( 5, x) + c,x tanh(cj y) 


+| 

cosh(c,x)cosh(5,y) cosh(c 2 x)cosh(s 2 y) ( J ( cosh(s 2 x)cosh(c 2 y) 


cosh(sjx)cosh(c,y) cosh^xjcoshfoy) cosh(^ 1 x)cosh(c 1 y) J 


— s 2 y tanh(s 2 x) + c 2 x tanh(c 2 y) 


+ 

coshfc 1 x)cosh(5 1 y) cosh(c 2 x)cosh(s 2 y) cosh(^,x)cosh(c 1 y) 

— — r H ; ; T7 - / \ w \ 


cosh(s 2 x)cosh(c 2 y) cosh(s 2 x)cosh(c 2 y) cosh(5 2 x)cosh(c 2 y) J 


(C-55) 


and it is noted that 

* = ~T J r(0cos(o> c f -+- 

<y wo 


and 



(C-56) 
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Appendix D 

EQUIVALENT NOISE DATA AND SIMULATION SOURCE CODE 


In Chapter 5 the equivalent noise in the high-SNR loop was 
discussed. To verify the derived numerical approach, it was 
compared to the theoretical result for QPSK (5-18). It was then 
applied to 8PSK and 16PSK. The numerical calculations were 
executed with the use of MATLAB source code. The 8PSK source 
code is presented here for review: 


%This program calculates the variance of the %equivalent 
noise term in the error signal of the 8PSK %HIGH-SNR MAP 
loop. 

clear 

format long 
%Test EsNo 
EsNo=10; 

%sigma2 is the I or Q. channel noise variance out 
%of either of the integrators 

sigma2=l/(2*10 A (EsNo/10)); 

%k is the constant out in front of the exponents in the 
%double integral for the equivalent noise calculation. 

k=l/(2*3.14159*sigma2); 

%ke is the constant in the exponent 

ke=l/(2*sigma2); 

%Clear the variable that tallies up the volume 
volume = 0; 

%To minimize computation time special measures were 
%taken to ensure that the resolution of the integration 
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%was "fine" enough but not more than was necessary 
%Basically, the calculation was repeated, increasing the 
%resolution each time, until the change in volume 
%dropped below 0.1 percent. 


%Set the length of one side of a pixel in the integration 
%grid. This is the integration resolution 


delta = .05; 

%Begin the double loop that integrates across the grid 


for 1=1:281, 
for j=l:281, 


%Convert the matrix indices to bilateral coordinates 
%that are designed to minimize the 
%needed. This minimization was done prior to this 
% co de being run. The size of the grid is designed so as 
%to contain 99.9 percent of the volume of the Gaussian 
%probability mountain. The size of the grid ts 
%therefore dependent on the SNR. Note that the 
%Gaussian probability mountain is centered at the zero 
%phase, modulation point. The ranges over whichx 
%and y are varied are the results of this minimization 

%process. 


%x range from -6 to 8 


x=(j-121)*delta; 

%y range from 7 to -7 


y=(141-i)*delta; 

%As the double loop sequences through each of the 
%locations in the x,y plane (the grid) the ^owmgif- 
%thens assign the data estimates for calculating the 
%equivalent noise. Actually the square of 1 and of CL 
o/nare needed as is the inner product 


if (y>=-0.41421*x)&(y<0.41421*x), 
ihat2=1.0; 
qhat2=0.0; 
ihatqhat=0.0; 

elseif (y>=0.41421*x)&(y<2.41421 x), 
ihat2=0.5; 
qhat2=0.5; 

ihatqhat=0.5; . 

elseif (y>=2.41421*x)&(y>-2.41421 x), 
ihat2=0.0; 
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qhat2=1.0; 

ihatqhat=0.0; 

elseif (y<=-2.41421*x)&(y>-0.41421*x), 
ihat2=0.5; 
qhat2=0.5; 
ihatqhat=-0.5; 

elseif (y<— 0.4 1 42 l*x)&(y>0.41 42 1 *x) , 
ihat2=1.0; 
qhat2=0.0; 
ihatqhat=0.0; 

elseif (y<=0.41421*x)&(y>2.41421*x), 
ihat2=0.5; 
qhat2=0.5; 
ihatqhat=0.5; 

elseif (y<=2.4142 l*x)&(y<-2.4142 l*x), 
ihat2=0.0; 
qhat2=1.0; 
ihatqhat=0.0; 

else, 

ihat2=0.5; 

qhat2=0.5; 

ihatqhat=-0.5; 

end 

%At location x,y calculate the joint noise density %function 

c=exp(-ke*(y) A 2)*exp(-ke*(x-1.0) A 2); 

%Calculate the square of the equivalent noise term 

da(i,j) = ((y) A 2*ihat2)+((x-1.0) A 2*qhat2) 
da(i,j) = da(i,j)-(2*(y)*(x-1.0)*ihatqhat); 

%Multiply the two together and then by the constant k 

d = da(i,j)*k*c; 

%"d" is the magnitude at any point in the grid. It needs 
%to be mulitiplied by the area of pixel to represent a 
%volume element. However, to save computation this 
%final step is withheld until the end 

%Accumulate the "volume" for the variance calculation 

volume = volume + d; 

end 

end 

%Multply by the pixel area and calculate the variance 
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var= 1 O*log 1 0 ( delta*delta*volume ) ; 

%Save data 
save dataout 

The theoretical and numerical data, plotted in Figure 5 1 and 
Figure 5-9 is presented here, in dB, in Table D-l and Table D-2. 
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Table D-l. The first set of MPSK equivalent noise variance data. 
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Table D-2. The second set of MPSK equivalent noise variance data. 



BPSK 

QPSK 

QPSK 

8PSK 

16PSK 


Theoretical 

Theoretical 

Numerical 

Numerical 

Numerical 

SNR 

Data 

Data 

Data 


■SlpMl 


mn 


(variance) 

I®2E1 


10.000 

-13.000 

-13.010 

-13.010 

-13.01 

-13.01 

10.500 

-13.500 

-13.510 




11.000 

-14.000 

-14.010 




11.500 

-14.500 

-14.510 




12.000 

-15.000 

-15.010 




12.500 

-15.500 

-15.510 




13.000 

-16.000 

-16.010 






-16.510 




14.000 

-17.000 

-17.010 




14.500 

-17.500 

-17.510 




15.000 

-18.000 


-18.012 

-18.078 

-18.087 

15.500 

-18.500 





16.000 

-19.000 





16.500 

-19.500 





17.000 

-20.000 





17.500 

-20.500 





18.000 

-21.000 





18.500 

-21.500 





19.000 

-22.000 





19.500 

-22.500 





20.000 

-23.000 


-23.012 

-23.11 

-23.097 
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Appendix E 

PHASE DETECTOR GAIN DATA AND SIMULATION SOURCE CODE 

In Chapter 5 the phase detector gain in the high-SNR loop 
was discussed. To verify the derived numerical approach, it was 
compared to the theoretical result for QPSK (5-64). It was then 
applied to 8PSK and 16PSK. The numerical calculations were 
executed with the use of MATLAB source code. The 8PSK source 
code is presented here for review: 

%This code generates the 8PSK phase detector gains 
%for a high-SNR MAP loop 

clear 

%Set the number of points to integrate the entire 
%phase probability over 

k=20000; 

snrdb 10.0; 

snr= 1 0 A ( snrdb/ 1 0 ) ; 

%Select Phase error that the phase detector is 
"devaluated at 

psi=0.001; 

%Clear storage variables 
popsk=[0 0 0 0 0 0 0 0]; 

%Start counter that integrates across the density 
for thetao=-k:l:k; 

%Create phase variable from counter (+/- pi) 
theta=thetao*pi/k; 

%The next three equation all combine to represent the 
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%density of the phase 

pl=2*sqrt(snr*pi)*cos(theta+psi)*exp(snr*cos(theta+psi)*cos 

(theta+psi)); 

p2=.5*(l-erf(-sqrt(snr)*cos(theta+psi))); 


p=(l/2/pi)*exp(-snr)*(l+pl*p2); 

%Now decide which decision region to accumulate the 
%density component in 


if (theta<0.392699082)&(theta>=-0.392699082), 

elseif ( diSa> =0. 3 9 2^9908 2^& ( theta< 1 . 1 7 809 7 245 ) , 

popsk(2)=popsk(2)+p; 

elseif ( theta>= 1.17 8097245 )&(theta< 1 .963495409) , 
popsk(3)=popsk(3)+p; 

elseif (theta>=1.963495409)&(theta<2.748893572), 
popsk(4)=popsk(4)+p; 

elseif ( theta> =2.748893572)l( theta<-2 .748893572), 
popsk( 5 ) =popsk( 5 ) +p; 

elseif (theta>«-2.748893572)&(theta<-l. 963495409), 
popsk(6)=popsk(6)+p; nnA ^^ r . 

elseif (theta>— 1.963495409)&(theta<-l. 178097245), 
noosk( 7 ) =popsk( 7 ) +p; 


else, 

end 


popsk( 8)=popsk( 8)+p; 


end 

%Having accumulated the probabilities for each of the 
%decision regions they must be scaled so that their 
%total eauals unity 


po 1 =pi/k*popsk( 1 ) ; 
po2=pi/k*popsk(2); 
po3 =pi/k*popsk( 3 ) ; 
po4=pi/k*popsk(4); 
po 5 =pi/k*popsk( 5 ) ; 
po6 = pi/k*popsk( 6 ) ; 
po7 =pi/k*popsk(7 ) ; 
po8=pi/k*popsk(8); 

%These probabilities are now used to compute the 
%phase detector characteristic output 


odo2=sin(psi)*pol+sin(psi+pi/4)*po2; 
pdo2=pdo2+sin(psi+pi/2)*po3+sin(psi+3 pi/4) po4, 
pdo2=pdo2+sin(psi+pi)*po5+sin(psi+5*pi/4) po6; 
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pdo2=pdo2+sin(psi+3*pi/2)*po7+sin(psi+7*pi/4)*po8; 


%The gain is the phase detector characteristic output 
%divided by the phase error. Also note that the square 
%of the gain is what is desired 

alpha2=20*logl0(pdo2/psi); 

%Save the data 

save 8pskdata; 


The theoretical and numerical data, plotted in Figure 5-11 and 
Figure 5-17 is presented here, in dB, in Table E-l and Table E-2. 
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Table E-l. The first set of MPSK phase detector gain data 




BPSK 

QPSK 

QPSK 

8PSK 

16PSK 


Theoretical 

Theoretical 

Numerical 

Numerical 

Numerical 

SISR 

Data 

Data 

Data 

Data 

Data 

(dB) 

(variance) 

(variance) 

(variance) 

(variance) 

(variance) 

n oool 

-1.486 

-14.033 

-14.029 




■smil 






■KEEBI 

mmm 

-11.660 



^n3ST!]| 

-0.846 

-10.533 




mmm 

-0.677 

-9.442 

-9.440 

-42.790 


■aTOII 

■Efim 






BBm 

eesshbi 

-7.399 

-7.397 

-37.067 



■Kssni 

-6.454 

mm 



4 000 

^BssSj 


■R^SI 



4 500 

-0.154 

-4.739 




5 000 

-0.104 

-3.976 

-3.974 

-26.469 


5 500 

-0.068 

-3.281 




6 000 

-0.042 

-2.657 

-2.656 

-21.699 


6.500 

-0.024 

-2.106 




7 000 

-0.014 

-1.628 

-1.627 

-17.350 


7 500 

~~ -0.007 1 

-1.224 




8 000 

-0.003 

-0.891 

-0.89CT 

-13.465 


8 500 

-0.002 

-0.625 




9 000 

-0.001 

-0.420 

-0.420 

-10.080 


9 500 

0.000 

-0.269 




10 000 

0.000 

-0.163 

-0.163 

-7.214 

-41 .304 

10 500 

0.000 

-0.093 




1 1 000 

0.000 

-0.049 

-0.049 

-4.877 

-33.383 

1 1 .500 

0.000 

-0.024 




1 2.000 

0.000 

-0.011 

0.011 

-3.061 

-26.388 

12.500 

0.000 

-0.004 




13.000 

o.ooc 

-0.002 

0.002 

-1.741 

-20.376 

13.50C 

> o.ooc 

-0.001 




14.00C 

) o.ooc 

) O.OOC 

) 

-0.86£ 

5 -15.324 

14.50C 

) o.ooc 

) O.OOC 

) 
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Table E-2. The second set of MPSK phase detector gain data. 



BPSK 

QPSK 

QPSK 

8PSK 

16PSK 







SNR 

Data 

Data 

Data 

Data 

Data 

(dB) 

m 



n 


15.000 

0.000 

0.000 


-0.364 

-1 1.162 

15.500 

0.000 

0.000 




16.000 

0.000 



-0.121 

-7.803 

16.500 

0.000 

0.000 




17.000 

0.000 

0.000 


-0.030 

-5.166 

17.500 

0.000 

0.000 




18.000 

0.000 

0.000 


-0.005 

-3.177 

18.500 

0.000 




. . . . _ 

19.000 

0.000 



0.000 

-1.769 

19.500 

0.000 





20.000 

0.000 

0.000 



-0.861 

20.500 

0.000 

0.000 




21.000 

0.000 

0.000 



-0.350 

21.500 

0.000 

0.000 




22.000 

0.000 

0.000 



-0.112 

22.500 

0.000 

0.000 




23.000 

0.000 

0.000 



-0.026 

23.500 

0.000 

0.000 




24.000 

0.000 




-0.004 

24.500 

0.000 





25.000 

0.000 




0.000 
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Appendix F 

MPSK VARIANCE OF PHASE ERROR RESULTS 


In Chapter 6 the equation that describes the variance of the 
phase error was discussed. Results were presented graphically. 
The components of the equation and the variance itself are 
tabulated here. Equation (6-7) is the equation used for computing 
the variance. Tables F-l and F-2 contain the BPSK data for the 4 
selected high-SNR loop bandwidth-to-symbol rate ratios. Table F- 
3 contains the QPSK variance data. The 8PSK numerical variance 
data is tabulated in Table F-4 and Table F-5 and the 16PSK 
numerical data is tabulated in F-6. All of the data are plotted in 

dB. 
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Table F-l. The first set of BPSK variance data. 


High SNR 
BL/SR=0.125% 
variance) 
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Table F-2. The second set of BPSK variance data 


SMR 

High SNR 

Hiqh SNR 

High SNR 

High SNR 

BL/SR=1% 

BL/SR=0.5% 

BL/SR=0.25% 

BL/SR=0.125% 

(dB) 

(variance) 

(variance) 

(variance) 

(variance) 

1 0 

-33.01 

-36.021 

-39.031 

-42.041 

9.5 

-32.51 

-35.521 

-38.531 

-41.541 

9 

-32.01 

-35.02 

-38.03 

-41.041 

8 

-31.008 

-34.019 

-37.029| 

-40.039 

7 

-30.002 

-33.013 

-36.023 

-39.033 

6 

-28.986 

-31.996 

-35.00o 

-38.017 

5 

-27.948 

-30.958 

-33.968 

-36.979 

4 

-26.878 

-29.888 

-32.899 

-35.909 

3 

-25.765 

-28.776 

-31.786 

-34.796 

2 

-24.602 

-27.612 

-30.623 

-33.633 

1 

-23.383 

-26.393 

-29.404 

T -32.414 

0 

-22.108 

-25.118 

-28.129 

-31.139 




Table F-3. The QPSK variance data. 




IHCS233I 


langsigH 


BL/SR=1% 

BL/SR=0.5% 



\mmm 

(variance) 



WisNM.mi 

40 

-63.01 

-66.021 

-69.031 

-72.041 

30 

-53.01 

-56.021 

-59.031 

-62.041 

25 

-48.01 

-51.021 

-54.031 

-57.041 

24 

-47.01 

-50.021 

-53.031 

-56.041 

23 

-46.01 

-49.021 

-52.031 

-55.041 

22 

-45.01 

-48.021 

-51.031 

-54.041 

21 

-44.01 

-47.021 

-50.031 

-53.041 

20 

-43.01 

-46.021 

-49.031 

-52.041 

1 9 

-42.01 

-45.021 

-48.031 

-51.041 

1 8 

-41.01 

-44.021 

-47.031 

-50.041 

1 7 

-40.01 

-43.021 

-46.031 

-49.041 

1 6 

-39.01 

-42.021 

-45.031 

-48.041 

1 5 

-38.01 

-41.021 

-44.031 

-47.041 

14 

-37.01 

-40.021 

-43.031 

-46.041 

1 3 

-36.009 

-39.019 

-42.03 

-45.04 

12 

-35.004 

-38.014 

-41.024 

-44.035 

1 1 

-33.981 

-36.991 

-40.001 

-43.012 

1 0 

-32.912 

-35.923 

-38.933 

-41.943 

9 

-31.758 

-34.769 

-37.779 

-40.789 

8 

-30.478 

-33.488 

-36.498 

-39.508 

7 

-29.039 

-32.05 

-35.06 

-38.07 

6 

-27.434 

-30.445 

-33.455 

-36.465 

5 

-25.666 

-28.676 

-31.686 

-34.696 

4 

-23.741 

-26.752 

-29.762 

-32.772 

3 

-21.671 

-24.682 

-27.692 

-30.702 

2 

-19.458 

-22.468 

-25.478 

-28.489 

1 

-17.099 

-20.11 

-23.12 

-26.13 

0 

-14.588 

-17.598 

-20.608 

-23.618 
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Table F-4. The first set of 8PSK variance data. 
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Table F-5. The second set of 8PSK variance data. 



memim 










Hffli 


MBHHI 

IHSESSHI 

1 0 

- 28.403 

- 31.414 

- 34.424 

- 37.434 

9.5 

- 26.959 

- 29.969 

- 32.979 

- 35.989 

9 

- 25.393 

- 28.403 

- 31.413 

- 34.424 

8.5 

- 23.7 

- 26.71 

- 29.72 

- 32.731 

8 

- 21.873 

- 24.883 

- 27.894 

- 30.904 

7.5 

- 19.905 

- 22.915 

- 25.926 

- 28.936 

7 

- 17.793 

- 20.803 

- 23.814 

- 26.824 

6.5 

- 15.537 

- 18.548 

- 21.558 

- 24.568 

6 

- 13.139 

- 16.149 

- 19.16 

- 22.17 

5.5 

- 10.606 

- 13.616 

- 16.626 

- 19.637 

5 

- 7.9448 

- 10.955 

- 13.965 

- 16.976 

4 

- 2.286 

- 5.2963 

- 8.3066 

- 11.317 

3 

| 3.7369 

0.72662 

- 2.2837 

- 5.294 

2 i 10.026 

7.0159 

4.0056 

0.9953 
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Table F-6. The 16PSK variance data. 



High SNR 

High SNR 

High SNR 

High SNR 

SNR 

BL/SR=1% 

BL/SR=0.5% 

BL/SR=0.25% 

BL/SR=0.125% 

(dB) 

(variance) 

(variance) 

(variance) 

(variance) 

MEftH 

-63.01 

-66.021 

-69.031 

-72.041 

30 

-53.01 

-56.021 

-59.031 

-62.041 

25 

-48.01 

-51.021 

-54.031 

-57.041 

24 

-47.01 

-50.021 

-53.031 

-56.041 

23 

-46.01 

-49.021 

-52.031 

-55.041 

22 

-45.01 

-48.021 

-51.031 

-54.041 

21 

-44.01 

-47.021 

-50.031 

-53.041 

20 

-42.491 

-45.501 

-48.511 

-51.522 

19.5 

-41.749 

-44.759 

-47.769 

-50.78 

19 

-40.935 

-43.945 

-46.956 

-49.966 

18.5 

-40.04 

-43.051 

-46.061 

-49.071 

1 8 

-39.056 

-42.066) 

-45.076 

-48.087 

17.5 

-37.971 

-40.981 

-43.991 

-47.002 

1 7 

-36.775 

-39.786 

-42.796 

-45.806 

16.5 

-35.457 

-38.467 

-41.478 

-44.488 

1 6 

-34.001 

-37.011 

-40.022 

-43.032 

15.5 

-32.39 

-35.401 

-38.411 

-41.421 

1 5 

-30.605 

-33.615 

-36.625 

-39.635 

14.5 

-28.62 

-31.63 

-34.641 

-37.651 

14 

-26.411 

-29.422 

-32.432 

-35.442 

13.5 

' -23.952 



-32.983 

1 3 

-21.218 

El 

-27.239 

-30.249 

12.5 

-18.19 

-21.2 

-24.211 

-27.221 

1 2 

^ -14.857 

-17.868 

-20.878 


■BE 

■HHIRfr 

-14.229 

-17.24 


EH 


-10.296 

-13.306 

-16.316 

Bee 


-6.0892 

-9.0995 

-12.11 
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Appendix G 

MPSK HIGH-SNR MAP CARRIER TRACKING SIMULATION RESULTS 

In Chapter 6 the simulation results were presented. 
Simulations were conducted to verify that the numerical results 
for the 8PSK and 16PSK variance of the phase error were correct. 
These simulations were conducted using the nonlinear baseband 
model of Figure 4-2. All of the baseband simulations assume 
phaselock at the modulation point corresponding to the zero- 
phase angle modulation symbol and assume that only the zero- 
angle modulation symbol is being transmitted. This is exactly 
what was done for the variance of the phase error component 
calculations of Chapter 5. Further, it speeds up the simulation time 
since random data is not required. 

One sample per symbol period is used in these time- 
simulations where the output at the end of a simulation run is the 
statistical variance of 35 samples of the VCO phase in the 
baseband loop. Since the "received" signal phase is held at a 
constant zero radians, the VCO output phase is indeed the phase 
error. Only 35 samples were used to save computational time. 

Each simulation actually executes several thousand 
iterations to obtain the 35 statistical samples. This is required 
because adjacent samples in time are highly correlated with each 
other. It is the nature of the phaselock loop to average out the 
effects of noise and thus vary the phase slowly in response. To 
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obtain statisically uncorrelated samples required analysis of the 

transfer function of the loop. 

The tracking loop is a filter whose transfer function, 
translated down to baseband, is H(s). The response of this filter to 

white Gaussian noise is a noise process whose normalized power 
spectral density magnitude response is | H(sf. The inverse Fourier 

transform of this response yields the normalized autocorrelation 
function of the process. What is of interest in this autocorrelation 
function are the points at which the function equates to zero. 
These points correspond to the time difference required between 
noise samples so that they have zero correlation with each other. 
The response, | \H(sf, is a lowpass function and its inverse Fourier 

transform yields a function that has a damped oscillation in both 
directions as the time difference increases. The first zero crossing 
represents the shortest amount of time required between 
uncorrelated samples. Therefore for all 4 PSK schemes, at 3 
different loop bandwidths, and for all the desired SNR's, the first 
zero crossing was obtained. This minimized the required run time 
of the simulations. 

The source code for the simulation of the 8PSK baseband 
carrier tracking loop is: 

%8PSK high-SNR MAP carrier tracking simulator 

clear 
tl=clock; 
format long 
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%Initialization parameters: natural freq, damping 
%factor, symbol rate 

wn=1171.88; 

zeta=1.2; 

sr=5e5; 


%Calculate filter integrator gain: a, integrator parmeter 
%c=a*T, high-SNR gain: k,vco integrator parameter; 

a=wn/(2*zeta); 

c=a/sr; 

k=2*wn*zeta; 

k2=k/sr; 


%Initialize carrier phase, vco initial phase, and loop 
%filter integrator 

sigphs=0; 

vcophs=0; 

fltint=0; 


%set random number distribution 
rand ('normal'); 

%Set statcnt size for snrs. This is the number of 
%samples used to compute the variance of the phase 
%error. It sets up the size of the outer loop. The size of 
%the inner loop is found from the auto correlation 
%function associated with white noise that has passed 
%through a lowpass function that has the same 
%transfer function as the loop, namely H(s). 

stattot=35; 

% For snr's from 5.5 to 15 dB in .5 dB steps and then 
%16 to 25 in 1 db steps the zero crossings occur at the 
%following simulation time sample increments 

statsize=[6280 

5495 

4867 

4239 

3768 

3454 

3140 

2826 

2512 

2355 

2198 
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2041 

1884 

1727 

1727 

1570 

1570 

1570 

1570 

1413 

1413 

1413 

1413 

1413 

1413 

1413 

1413 

1413 

1413 

1413 

1413]; 

%The snr's that the simulations will be run at 

snrv=[5.5 

6 

6.5 

7 

7.5 

8 

8.5 

9 

9.5 

10 

10.5 
11 

11.5 
12 

12.5 

13 

13.5 

14 

14.5 

15 

16 

17 

18 

19 

20 
21 
22 
23 
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24 

25]; 


%outer loop to count through the snr array 

for snrcnt = 1:1:30; 

snrdb = snrv(snrcnt); 
snr= 10 A (snrdb/10); 

%noise STDEV out of correlators 

sigma= sqrt(l/(2*snr)); 

%set carrier phase, vco initial phase, loop filter 
%integrator variables and phase error 

sigphs=0; 

vcophs=0; 

fltint=0; 

phserr=0; 

fltlast=0; 

%Reset stats for calculating variance 

sumerr=0; 

sqrerr=0; 

maxerr=0; 

%Initialize random number generator 
rand('seed',650); 

%Statistics loop. This loop iterates once for each sample 
%used in the variance calculation 

for statcnt = l:l:stattot; 

%The loop that sets the number of samples used 
%between the statistical samples 

for simcnt = l:l:statsize(snrcnt); 

%Fix data 

im=1.0; 

qm=0.0; 


%Generate noise 
ni=sigma*rand( 1 ); 
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nq=sigma*rand( 1 ) ; 

%Calculate correlator outputs 

i=im*cos(phserr) +qm*sin(phserr ) +ni; 
q=-im*sin(phserr ) +qm*cos(phserr)+nq; 

%Make data estimates 

rcvdphs= atan2 ( q,i) ; 

if (rcvdphs<0.392699082)&(rcvdphs>=-0.392699082), 

ihat=1.0; 

elseif (rcvdphs<l . 1 7 8097 245)&(rcvdphs>=0.3 92699082 ) , 
ihat=0.707 106781; 

qhat=0.707 106781; . 

elseif (rcvdphs<1.963495409)&(rcvdphs>=l. 178097245), 
ihat=0.0; 

elseif (rcvdphs<2. 748893 572)&(rcvdphs>=1.963495409), 
ihat=-0.707 1 067 8 1 ; 

elseif ( rcv^phs<>2. 7488935 7 2 ) I ( rcvdphs>=2 .748893 572), 
ihat=-1.0; 

elseif (n^dphs< , - 1 - 963495409)&(rcvdphs>= ' 2 ‘ 748893572) ’ 

ihat=-0.707106781; 

elseif (^^phs<-^ ? l 7^097245 )&(rcvdphs>=-1.963495409), 
ihat=0.0; 
qhat=-1.0; 

else 

ihat=0. 707 106781; 
qhat=-0.707 106781; 

end 

%Generate error signal 
error=i*qhat-q*ihat; 

%Calculate filter output 
filtout= fltint+error; 

%Calculate filter integrator's new value for next time 
fltint=fltint+error*c; 

%Calculate next value of the vco phase 


vcophs=vcophs+fltlast*k2 ; 

%Set next vco input value; 
fltlast=filtout; 

%Calculate next value of the phase error 

phserr=-vcophs; 

%Store max phase encountered 

if ( abs ( phserr) >maxerr ) , 

maxerr=abs(phserr); 

end 

end 

%Iterate statistical parameters 

sqrerr=sqrerr+phserr A 2 ; 
sumerr=sumerr+phserr; 

end 

meanphs(snrcnt)=sumerr/stattot; 

vp=sqrerr/ (stattot-l)-sumerr A 2/( stattot* ( stattot- 1 ) ) ; 

varphs ( snrcnt ) = vp ; 

maxsa ve ( snrcnt ) =m axerr ; 

end 

%check time and save data 

t2=clock; 

tdiff=t2-tl; 

save 8pkl650b wn zeta sr meanphs varphs maxsave 


The simulation data is tabulated in the following tables for BPSK, 
QPSK, 8PSK, and 16PSK. The simulations were conducted using 
three specific high-SNR loop bandwidths (the same three used in 
the hardware) and a symbol rate of 500 kilo-samples-per-second. 
The ratios of the three specific bandwidths to the above symbol 
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rate are listed on the simulation data graphs of Figures 6-6 
through 6-17 and at the tops of the following data tables. 

The BPSK simulation data are tabulated in Table G-l. The 
QPSK data is contained in Table G-2 and Table G-3 contains the 
8PSK variance data. The 16PSK data are tabulated in Table G-4. 
All the data are plotted in dB. 


Table G-l. The BPSK simulation variance data. 



High SNR 

High SNR 

High SNR 

gsR 

BL/SR=0.19% 

BL/SR=0.33% 

BL/SR=0.66% 

(dB) 

(variance) 

(variance) 

(variance) 

25 



-49.67 

20 

-50.31 

-50.31 

-44.67 

1 5 

-45.31 

-45.31 

-39.67 

10 

-40.31 

-40.31 

-34.67 

5 

-35.74 

-35.74 

-29.98 

4 

-34.67 

-34.67 

-28.99 

3 

-33.14 

-33.14 

-27.26 

2 

-32.06 

-32.06 

-26.25 

1 

1 -30.89 

-30.89 

-25.39 

0 

-29.15 

-29.15 

-24.82 
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Table G-2. The QPSK simulation variance data. 



mill, iii ii ‘,i 

WEMUSiM 

1— 



P:lBdsB«ycfctli 




vnlnfRffini 

■■zuffERnni 

25 



- 49.67 

20 

- 50.11 

- 47.98 

- 44.67 

1 5 

- 45.31 

- 42.99 

- 39.67 

1 0 

- 40.44 

- 37.74 

- 34.74 

9 

- 39.70 

- 37.12 

- 33.77 

8 

- 36.88 

- 35.75 

- 32.73 

7 

- 35.06 

- 33.64 

- 30.76 

6 

- 35.34 

- 32.74 

- 30.27 

5 

- 34.45 

- 30.91 

- 29.00 

4 

- 31.85 

- 27.77 

- 27.31 

3 

- 29.16 

- 26.13 

- 22.89 

2 

- 28.88 

- 26.42 

- 22.21 

1 

- 24.01 

- 22.88 

- 19.11 

0 

- 20.75 

- 21.14 

- 16.87 
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Table G-3. The 8PSK simulation variance data. 
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Table G-4. The 16PSK simulation variance data 






SNR 


BL/SR=0.33% 

BL/SR=0.66% 

(dB) 

HBS53B 

(variance) 


25 



-49.67 

24 



-48.67 

23 



-47.7 

22 


-50.3 

-46.61 

21 


-48.87 

-45.52 

20 

-49.93 

-46.81 

-44.48 

19.5 

-49.2 

-46.63 

-44.23 

1 9 

-48.43 

-46.69 

-44 

18.5 

-47.34 

-45.14 

-40.66 

18 

-46.76 

-44.7 

-39.92 

17.5 

-45.5 

-43.57 

-39.76 

1 7 

-44.33 

-43.37 

-38.18 

16.5 

-44.13 

-42.02 

-38.15 

1 6 

-41.72 

-41.3 

-37.4 

15.5 

-41.83 

-39.36 

-36.81 

1 5 

-39.14 

-37.39 

-33.75 

14.5 

-37.11 

-35.41 

-33 

1 4 

-35.85 

-33.51 

-30.71 

13.5 

-31.57 

-30.03 

-29.17 
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Appendix H 

THE HARDWARE TEST CONFIGURATION AND RESULTS 

As discussed in Chapter 6 an MPSK high-SNR carrier 
tracking loop was constructed in hardware to verify the variance 
results. A simulated MPSK transmission was created and used as 
an input to the loop. In fact, just as was done with the numerical 
anaylsis and the computer simulations, the MPSK transmission 
represented the transmission of the same symbol continuously. 
This made the MPSK signal a continuous wave (CW) carrier signal 
that could be obtained from a standard frequency synthesizer and 
it removed the impact of intersymbol interference as a dependent 
variable in the testing. White noise was added to the CW carrier 
signal and the result passed to the demodulator front-end of the 
carrier tracking loop. This demodulator front end is a quadrature 
structure that power-divides the CW carrier-with-noise signal into 
two signals that are mixed in quadrature with the "VCO' signal. 
Actually a numerically controlled oscillator (NCO) was used in 
place of a VCO since much of the loop was constructed with digital 
hardware. The outputs of the quadrature mixers are fed through a 
set of lowpass filters. The use of lowpass filters instead of 
integrate-and-dump circuits does not cause a problem since the 
test set-up is calibrated at the output of the filters. Analog-to- 
digital conversion takes place after the lowpass filters and the 
generated 8 bits of I channel baseband data and 8 bits of Q_ 
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channel baseband data are passed to a set of electronically 
programmable read only memory (EPROM) chips. These EPROM 
chips generate the error signal of the high-SNR loop using 16 bits. 
This error signal is processed by a digital filter that creates a 
second order transfer function for the loop. The output from 
digital filter, a 16-bit word, is added to a 32-bit word that 
represents the quiescent frequency of the output sinusoid from 
the loop NCO. This 32-bit sum is the frequency that the NCO will 
output during the next symbol period and it is transferred byte- 
serially to a Stanford Telecomm 1172/1272 NCO/ evaluation board 
that generates both the sine and cosine of the desired output 
frequency. These two sinusoids are analog when they leave the 
evaluation board and return to the demodulation front-end where 
they are mixed in quadrature with the "received" signal. This 
circuit is shown in Figure H-l. 

The I and Q_ baseband data at the output of the lowpass 
filters are corrupted by the input noise and thus variance 
measurements cannot be made using these data. To facilitate 
making variance measurements the outputs from the NCO are 
mixed in quadrature with the "clean" CW carrier from the 
synthesizer. The mixer outputs are lowpass filtered to remove the 
double frequency term and create I and Q. baseband data that are 
corrupted only with phase jitter from the NCO. 

This I and Qdata are displayed on a modulation analyzer 
that can plot the I and Qdata in the constellation space. If there is 
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no noise and if the carrier tracking loop is phaselocked the display 
will show a single lock point. If there is noise present and 
therefore jitter on the NCO signals then the lock point will take the 
shape of an arc. The arc is a plot of past samples of the phase 
error. It is brightest in the center, at the lock point, indicating that 
the mean of the phase error is zero. The modulation analyzer 
samples the I and Q. components 4096 times over a period of 
approximately 30 seconds and calculates the statistical variance of 
the samples. By taking such a large sample set it was hoped that 
any correlation between adjacent samples would be averaged out. 
The results seem to indicate this is indeed the case. All the 

tabulated data that follows is plotted in dB. 

Table H-l and Table H-2 contain the hardware results for 
BPSK that were plotted in Figure 6-18. The BPSK data plotted in 
Figure 6-19, for the second loop bandwidth is tabulated in Table 
H-3 and H-4. The third bandwidth that was tested for BPSK 
yielded the data plotted in Figure 6-20 and tabulated in Tables H- 

5 and H-6. 

For QPSK, H-7 and H-8 contain the hardware results for the 
smallest high-SNR loop bandwidth-to-symbol rate ratio. This is 
the data shown in Figure 6-21. The data in Figure 6-22 is 
tabulated in Tables H-9 and H-10, and the data in Figure 6-23 can 
be found in Tables H-l 1 and H-l 2, corresponding to the widest 
high-SNR loop-to-symbol rate ratio tested. 
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Table H-13 and Table H-14 contain the hardware results for 
8PSK that were plotted in Figure 6-24. Table H-15 is the tabulated 
data for the results in Figure 6-25 and the data corresponding to 
the largest high-SNR loop bandwidth-to-symbol rate ratio, plotted 
in Figure 6-26, is tabulated in Tables H-16 and H-17. 

For 16PSK the data shown in Figure 6-27 is tabulated in 
Table H-18 and the data of Figure 6-28 is contained in Table H-19. 
The final two tables, Table H-20 and FI-21, contain the the data for 
the largest high-SNR loop bandwidth-to-symbol rate ratio. This is 
the data shown in Figure 6-29. 
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H-l. The MPSK high-SNR carrier tracking hardware test configuration. 

















Table H-l. The first set of BPSK hardware variance data at 0.19% 
high-SNR loop bandwidth-to-symbol rate ratio. 


LAB 



SNR for 

iKiilHSffnffiHi 


\WsS 

Measured 


■■"raBIlSW 


Theoretical 

■sun 

Variance 


BB21031 


variance 

32.9 

-43.775 

111 ' 

35 



27.9 

-44.223 


30 

-43.898 

-60.232 

25.4 

-43.736 

27.5 




-43.59 

&|te; 

25 

-43.685 

-55.232 




24.5 



21.9 

-43.673 


24 

-43.607 

-54.232 

21.4 

-43.644 


23.5 






20.9 

-43.623 

toil 

23 

-43.51 

-53.232 

20.4 

-43.004 


22.5 



19.9 

-43.405 

22 

-43.392 

-52.232 

19.4 

-42.997 

21.5 



18.9 

-43.401 


CVJ 

-43.248 

-51.232 

18.4 

-42.93 


20.5 



17.9 

-42.679 

0W^£':i , 

20 

-43.072 

-50.232 

17.4 

-42.519 


19.5 



16.9 

-42.461 


1 9 

-42.861 

-49.232 

16.4 

-42.471 


18.5 



15.9 

-42.169 


18 

-42.609 

-48.232 

15.4 

-41.869 

■ - V k 

17.5 



14.9 

-41.604 

\ ; . 

1 7 

-42.312 

-47.232 

14.4 

-41.538 

:V: , 

16.5 



13.9 

-41.201 


1 6 

-41.964 

-46.232 

13.4 

-40.599 


15.5 



12.9 

-40.321 1 


15 

-41.562 

-45.232 

12.4 

-39.876 


14.5 

-41.34 

-44.732 

11.9 

-39.865 


14 

-41.104 

-44.232 

11.4 

-39.644 1 

% - 

13.5 

-40.853 

-43.732 

10.9 

-39.137 

f[| II : 

13 

-40.589 

-43.232 

10.4 

-38.874 j 

V ' i t> ■ 

12.5 

-40.309 

-42.732 

9.9 

-38.481 j 

f 

12 

-40.016 

-42.232 

9.4 

-38.033 


11.5 

-39.709 

-41.732 
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Table H-2. The second set of BPSK hardware variance data at 
0.19% high-SNR loop bandwidth-to-symbol rate ratio. 


LAB 

SNR for 

Theoretical 


SN=t 

Measured 

Theoretical 

variance with 

Theoretical 

(dB) 

Variance 

Data (dB) < 

added noise floor 

variance 

ft Q 

-37.538 

1 1 

-39.389 

-41.232 


-37.276 

10.5 

-39.055 

-40.732 



10 

-38.709 

-40.232 


- ft ft 474 

9.5 

-38.351 

-39.732 



9 

-37.981 

wm&m 



8.5 

-37.601 

-38.732 


-35.1 48 

8 

-37.209 

-38.23 



7.5 

-36.811 

-37.732 


-33.937 

7 

-36.395 

-37.223 


-33.854 lillli 

6.5 

-35.985 

-36.732 


-33.308 

6 

-35.538 1 

-36.206 


-32.568 J 

5.5 

^ -35.128 

-35.732 

2.9 

-31.852BBH 

5 

-34.632 

-35.166 

2.4 

-3 1.857 >3$$^ 

4.5 


SgB 

1 .9 


4 

-33.671 

I^EZXjUJI 

1 .4 

-30.865 

3.5 




-30 

3 

-32.645 

-32.975 

mm 

-29.315f:|||Il 

2.5 


_ 

BsB 

-28.551 lllSlf 

2 

-31.549 

-31.804 


-28.171' {\~ 

1.5 




-27.231 ttfillg! 

1 1 

-30.381 

-30.574 

1 -1 6 

-26.453 

| 0.5 

t 



-26.101 

1 0 

i -29.141 

-29.285 
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Table H-3. The first set of BPSK hardware variance data at 0.33% 
high-SNR loop bandwidth-to-symbol rate ratio. 


LAB 



SNR for 

Theoretical 


SNR 

Measured 


Theoretical 

variance with 

Theoretical 

(dB) 

Variance 


Data (dB) 

added noise floor 

variance 

30.4 

-43.952 

lliiilllll 

32.5 



27.9 

-43.693 


30 

-43.776 

-57.824 

25.4 

-43.976 


27.5 



22.9 

-43.878 


25 

-43.421 

-52.824 

20.4 

-42.771 


22.5 



19.4 

~ 4- 2 . 4 6 

21.5 



18.9 

-41.89 


21 

-42.726 

-48.824 

18.4 

-41.844 


20.5 

-42.458 

-47.824 

17.9 

-41.376 


20 



17.4 

-41.262 


19.5 

-42.143 

-46.824 

16.9 

-40.806 

l • ^ ■ 

1 9 



16.4 

-40.844 


18.5 

-41.776 

-45.824 

15.9 

-40.543 

IllSlf:;: 

1 8 



15.4 

-39.92 


17.5 

-41.355 

-44.824 

14.9 

-39.651 

fiillll 

17 



14.4 

-39.41 9 

16.5 

-40.876 

-43.824 

13.9 

-39.066 

. . . 

1 6 



13.4 

-38.712 


15.5 

-40.34 

-42.824 

12.9 

-38.463 

i ; , ; ; : 

15 

-40.051 

-42.324 

12.4 

-37.94 

: : :-y 

14.5 

-39.748 

-41.824 

11.9 

-37.376 

:jT', 

1 4 

-39.431 

-41.324 

11.4 

-37.271 

: ;.y. " 

13.5 

-39.101 

-40.824 

10.9 

-36.902 


13 

-38.759 

-40.324 

10.4 

-36.069 

:: .. ... ; ; 

12.5 

-38.404 

-39.824 

9.9 

-35.815 


12 

-38.038 

-39.324 

9.4 

-35.528 


11.5 

-37.661 

-38.824 

8.9 

-35.014 


1 1 

-37.273 

-38.324 

8.4 

-34.561 

: - y 'i 

10.5 

-36.876 

-37.824 

7.9 

-34.156 


10 

-36.47 

-37.324 

7.4 

-33.601 

: : 

9.5 

-36.054 

-36.824 

6.9 

-33.033 


9 

-35.632 

-36.324 
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Table H-4. The second set of BPSK hardware variance data at 
0.33% high-SNR loop bandwidth-to-symbol rate ratio. 


LAB 

SNR for 

Theoretical 


SNR 

Measured 

Theoretical 

variance with 

Theoretical 

(dB) 

Variance 

Data (dB) 

added noise floor 

variance 

6.4 

-32.483 

8.5 

-35.201 

-35.822 

5.91 

-32.081 ■> fat* 

8 

-34.766 

-35.324 

5.4 

-31.793l^®r ; 

7.5 

-34.316 

-34.817 

4.9 

-31.512 

7 

-33.875 


4.4 

-30.824 

6.5 

-33.4 

wmm\ 

3 9 

-30.217®^^ 

6 

-32.964 

-33.324 

3.4 

-29.643 j|l|8 §|j§Si 

5.5 

-32.446 

-32.764 

~ 

2.9 

-29.231^^^1 

5 



2.4 

-29.029 jMmjKfc 

4.5 

-31.447 

-31.698 

1 .9 

-27.975 1 

4 



1 .4 

-27.667 

3.5 

-30.394 

-30.59 

0.9 

-26.965 

3 



0.4 

-26.345|l|li;; c 

2.5 

-29.283 

-29.434 

-0.1 

-25.66 '*- 

2 



-0.6 

-25.077 

1.5 

-28.11 

-28.225 

-1.1 

-24.51 

1 



-1.6 

-23.864 

0.5 

-26.877 

-26.963 
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Table H-5. The first set of BPSK hardware variance data at 0.66% 
high-SNR loop bandwidth-to-symbol rate ratio. 



Theoretical 



SNR 

variance with 

Theoretical 

Measured 

(dB) 

added noise floor 

variance 

variance 

30 

-43.469 

-54.814 

-43.464 

29 



-43.529 

28 



-43.372 

27 



-43.193 

26 



-43.036 

25 

-42.83 

-49.814 

-42.849 

24.5 



-42.711 

24 

-42.61 

-48.814 

-42.573 

23.5 



-42.658 

23 

-42.349 

-47.814 

-42.619 

22.5 



-42.148 

22 

-42.04 

-46.814 

-41.995 

21.5 



-41.988 

21 

-41.681 

-45.814 

-41.672 

20.5 



-41.395 

20 

-41.267 

-44.814 

-40.893 

19.5 



-40.909 

1 9 

-40.797 

-43.814 

-40.456 

18.5 



-40.325 

1 8 

-40.269 

-42.814 

-40.005 

17.5 



-39.652 

1 7 

-39.684 

-41.814 

-39.384 

16.5 



-38.788 

1 6 

-39.045 

-40.814 

-38.664 

15.5 



-38.093 

1 5 

-38.354 

-39.814 

-37.781 

14.5 

-37.991 

-39.314 

-37.435 

14 

-37.617 

-38.814 

-37.036 

13.5 

-37.232 

-38.314 

-36.779 

1 3 

-36.838 

-37.814 

-36.477 

12.5 

-36.434 

-37.314 

-35.788 
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Table H-6. The second set of BPSK hardware variance data at 
0.66% high-SNR loop bandwidth-to-symbol rate ratio. 



Theoretical 



SNR 

variance with 

Theoretical 

Measured 

(dB) 

added noise floor 

variance 

variance 

12 

-36.021 

-36.814 

-35.537 

11.5 

-35.601 

-36.314 

-35.061 

1 1 

-35.173 

-35.814 

-34.53 

10.5 

-34.738 

-35.314 

-34.015 

10 

-34.297 

-34.814 

-33.715 

9.5 

-33.851 

-34.314 

-33.163 

9 

-33.398 

-33.813 

-32.832 

8.5 

-32.942 

-33.314 

-32.384 

8 

-32.479 

-32.812 

-31.836 

7.5 

-32.016 

-32.314 

-31.341 

7 

-31.54 

-31.806 

-30.721 

6.5 

-31.075 

-31.314 

-30.628 

6 

-30.578 

-30.79 

-29.845 

5.5 

-30.123 

-30.314 

-29.298 

5 

-29.586 

-29.754 

-28.683 

4.5 



-28.53 

4 

-28.556 

-28.687 

-27.837 

3.5 



-27.167 

3 

-27.477 

-27.58 

-26.747 

2.5 



-26.076 

2 

-26.345 

-26.424 

-25.76 

1.5 



-25.243 

1 

-25.155 

-25.215 

-24.479 

0.5 



-23.781 

0 

-23.908 

1 -23.953 

1 -23.279 



Table H-7. The first set of QPSK hardware variance data at 0.19% 
high-SNR loop bandwidth-to-symbol rate ratio. 



Theoretical 



STvR 

variance with 

Theoretical 

Measured 

(dB) 

added noise floor 

variance 

variance 

30 

-41.047 

-60.232 

-41.002 

29 



-41.101 

28 



-41.184 

27 



-41.045 

26 



-40.962 

25 

-40.935 

-55.232 

-40.657 

24 

-40.894 

-54.232 

-40.472 

23 

-40.842 

-53.232 

-40.504 

22.5 



-40.83 

22 

-40.778 

-52.232 

-40.667 

21.5 



-40.477 

21 

-40.698 

-51.232 

-40.593 

20.5 



-40.492 

20 

-40.6 

-50.232 

-40.478 

19.5 



-40.342 

1 9 

-40.479 

-49.232 

-40.178 

18.5 



-40.281 

1 8 

-40.332 

-48.232 

-40.153 

17.5 



-39.854 

17 

-40.153 

-47.232 

-39.925 

16.5 



-39.922 

1 6 

-39.938 

-46.232 

-39.628 

15.5 



-39.322 

1 5 

-39.682 

-45.232 

-39.302 

14.5 



-39.165 

1 4 

-39.379 

-44.232 

-39.104 

13.5 



-38.726 

1 3 

-39.026 

-43.23 

-38.656 

12.5 



-38.394 

1 2 

-38.616 

-42.225 

-38.063 

11.5 



-37.943 
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Table H-8. The second set of QPSK hardware variance data at 
0.19% high-SNR loop bandwidth-to-symbol rate ratio. 



Theoretical 



SNR 

variance with 

Theoretical 

Measured 

(dB) 

added noise floor 

variance 

variance 

1 1 

-38.14 

-41.201 

-37.666 

10.5 




10 

-37.661 

-40.129 







-36.894 

-38.968 

-36.303 

8.5 



-35.746 

8 

-36.046 

-37.673 

-35.393 




-35.095 


_ -34.992 

-36.212 

-34.354 

6.5 



-33.77 

6 

-33.701 

-34.574 

-32.996 

5.5 



-32.469 

5 

-32.166 

-32.76 

-31.553 

4.5 



-31.134 

4 

-30.392 

-30.778 

-30.111 

3.5 



-29.217 

3 

-28.398 

-28.638 

-27.976 

2.5 



-26.731 

2 

-26.2 

-26.343 

- 25.37 

1.5 



-24.211 

1 

-23.813 

-23.895 

-23.505 

0 

-21.243 

-21.288 
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Table H-9. The first set of QPSK hardware variance data at 0.33% 
high-SNR loop bandwidth-to-symbol rate ratio. 


LAB 



SNR for 



SfsF! 

Measured 




Theoretical 

(dB) 

Variance 




variance 

42.9 

-41.589 

fill,. 

30 

-41.107 

-57.824 

32.9 

-41.234 


25 

-40.911 

-52.824 

27.9 

-41.291 


24 

-40.839 

-51.824 

22.9 

-41.115 


23 

-40.751 

-50.824 

21.9 

-40.665 


22 

-40.641 

-49.824 

20.9 

-40.555 


21 

-40.508 

-48.824 

19.9 

-40.491 111:11. 

20 

-40.345 

-47.824 

18.9 

-40.136 


19 

-40.149 

-46.824 

17.9 

-39.788 


18 

-39.913 

-45.824 

17.4 

-39.737 

w 

17 

-39.634 

-44.824 

16.9 

-39.619 

1 6 

-39.307 

-43.824 

16.4 

-39.256 

j; ;’;V '1 ' 1; 

1 5 

-38.926 

-42.824 

15.9 

-39.346 


1 4 

-38.491 

-41.824 

15.4 

-38.803 

IMIil 

1 3 

-37.997 

-40.823 

14.9 

-38.703 


1 2 

- 37.444 

-39.818 

14.4 

-38.449- 

1 1 

-36.823 

-38.796 

13.9 

-38.341 

■ 

1 0 

-36.169 

-37.73 

13.4 

-38.013 


9 

-35.295 

-36.583 

12.9 

-37.369 


8 

-34.319 

-35.315 

12.4 

-37.549 

7.5 



11.9 

-37 

■ 

7 

-33.157 

-33.898 

11.4 

-36.501 

•<l . ■ l \ 

6.5 



10.9 

-36.291 


6 

-31.795 

-32.324 

10.4 

-35.905 


5.5 



9.9 

-35.418 

: . 

5 

-30.236 

-30.599 

9.4 

-35.243 


4.5 



8.9 

-34.63 

< : 1 : 

4 

-28.491 

-28.731 

8.4 

-34.217 


3.5 



7.9 

-33.549 


3 

-26.579 

-26.731 

7.4 

-33.167 


2.5 



6.9 

-32.637 


2 

-24.509 

-24.603 
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Table H-10. The second set of QPSK hardware variance data at 
0.33% high-SNR loop bandwidth-to-symbol rate ratio. 







Table H-ll. The first set of QPSK hardware variance data at 0.66% 
high-SNR loop bandwidth-to-symbol rate ratio. 



Theoretical 



SNR 

variance with 

Theoretical 

IHSSS1 

(dB) 

added noise floor 

variance 

| 

30 

-40.823 

-54.814 

-40.821 

29.5 



-40.909 

29 



-40.754 

28.5 



-40.593 

28 



-40.576 

27.5 



-40.735 

27 



-40.703 

26.5 



-40.51 

26 



-40.528 

25.5 



-40.48 

25 

-40.464 

-49.814 

-40.118 

24.5 



-40.548 

24 

-40.335 

-48.814 

-40.365 

23.5 



-40.038 

23 

-40.178 

-47.814 

-40.215 

22.5 



-39.852 

22 

-39.989 

-46.814 

-39.989 

21.5 



-39.826 

21 

-39.761 

-45.814 

-39.597 

20.5 



-39.52 

20 

-39.491 

-44.814 

-39.002 

19.5 



-39.021 

1 9 

-39.173 

-43.814 

-38.84 

18.5 



-38.631 

1 8 

-38.803 

-42.814 

-38.466 

17.5 



-38.287 

1 7 

-38.378 

-41.814 

-37.868 

16.5 



-37.699 

1 6 

-37.896 

-40.814 

-37.459 

15.5 



-37.39 

1 5 

-37.356 

-39.814 

-36.56 
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Table H-12. The second set of QPSK hardware variance data at 
0.66% high-SNR loop bandwidth-to-symbol rate ratio. 



Theoretical 



SNR 

variance with 

Theoretical 

Measured 

(dB) 

added noise floor 

variance 

variance 

1 4.5j 



-36.672 

14 

-36.76 

-38.814 

-36.104 

13.5 



-35.698 

1 3 

-36.11 

-37.813 

-35.61 

12.5 



-35.101 

12 

-35.406 

-36.808 

-34.917 

11.5 



-34.376 

1 1 

-34.643 

-35.786 

-33.969 

10.5 



-33.533 

1 0 

-33.834 

-34.72 

■EiB 

9.5 

■mmm 


-32.782 

9 

liMMBMTTTI 

-33.573 

-32.225 

8.5 

— jMfBM 



8 

-31.755 

-32.30o 

-31.064 

B 



-30.551 

B 

-30.484 

-30.888 

-29.88 

6.5 



-29.554 

6 

-29.029 

-29.314 

-28.639 

5.5 



-27.87 

5 

-27.395 

-27.588 

-27.069 

4.5 



-26.656 

4 

-25.593 

-25.72 

-25.52 

3.5 



-24.494 

3 

j -23.641 

-23.721 

-23.368 

2.5 



-22.224 

2 

-21.543 

i -21.593 

-21.034 

1.5 



-19.804 

1 

-19.302 

! -19.332 


0 

-16.91 

T -16.927 
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Table H-13. The first set of 8PSK hardware variance data at 0.19% 
high-SNR loop bandwidth-to-symbol rate ratio. 


LAB 



SNR for 

Theoretical 


SNR 

Measured 



Bom 

Theoretical 

(dB) 

Variance 




variance 

42.9 

-42.266 iililll 30 

-42.329 

-60.232 

32.9 

-42.532 

1 c i''* 

I 3 !??? • ' -p: : ' i .■ ■ 
til f , $ 

25 

-42.179 

-55.232 

27.9 

-42.383 

24 

-42.124 

-54.232 

26.9 

-42.293 

23 

-42.055 

-53.232 

25.9 

-42.085 iiiiii* 22 

-41.97 

-52.232 

24.9 

-41.556 

i-> : * '' 

21 

-41.866 

-51.232 

23.9 

-41.973 

20 

-41.738 

-50.232 

22.9 

-41.486 

1 9 

-41.581 

-49.232 

22.4 

-41.068 

1 8 

-41.392 

-48.228 

21.9 

-41.668 filiSll 17 

-41.161 

-47.213 

21.4 

-41.632 


16 

-40.874 

-46.156 

20.9 

-41.943 

r V: -yy 

15 

-40.499 

-45.002 

20.4 

-41.42 


14.5 

-40.263 

-44.368 

19.9 

-41.154 

■ ■' 

14 

-39.984 

-43.682 

19.4 

-41.414 

13.5 

-39.649 

-42.935 

18.9 

-41.278 

1 3 

-39.248 

-42.12 

18.4 

-41.084 

12.5 

-38.765 

-41.229 

17.9 

-40.961 

1 2 

-38.185 

-40.253 

17.4 

-40.827 

11.5 

-37.493 

-39.188 

16.9 

-40.932 

■ ■; : - r v : - : . : 

1 1 

-36.673 

-38.024 

16.4 

-40.587 

10.5 

-35.709 

-36.757 

15.9 

-41.133 

1 0 

-34.591 

-35.377 
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Table H-14. The second set of 8PSK hardware variance data at 
0.19% high-SNR loop bandwidth-to-symbol rate ratio. 


LAB 

SNR for 

Theoretical 


SNR 

Measured Theoretical 

variance with 

Theoretical 

■ESI 

Variance Data (dB) 

added noise floor 

variance 

■121 

mrtn 

-40.573 9.5 

-33.307 

-33.878 

-40.31 9 

-31.851 

-32.252 

mm 

■ran 

-40.038 MSP 8.5 

-30.222 

-30.493 

-39.918^^» : 8 

-28.419 


-39.414^^® 7.5 

-26.444 


1 IP 0 


-24.303 

-24.37 

mtm 

mnm 

-38.66 iilaijii 6.5 

-22.004 

-22.044 

-37.924 MM* 6 

-19.5571 

-19.579 

mm 

-37.293 5.5 

-16.973 

-16.986 

10.9 

- 3 R Q 1 2 



10.4 

-35.193^PW5f 



9.9 
9.4 

8.9 

-34.191 'MPlS-w 



-33.189^^^1 



-31.169|Mk 



8.4 

7.9 

-29.823Mifel?! — 



-27.1 21 IfSllllJ 


J 
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Table H-15. 8PSK hardware variance data at 0.33% high-SNR loop 
bandwidth-to-symbol rate ratio. 


LAB 



SNR for 

Theoretical 


SNR 

Measured 


Theoretical 

variance with 

Theoretical 

■0D1 

Variance 

— 


added noise flooi 

variance 


-44.041 


40 

-43.982 

-67.824 


-43.395 

| | 

30 

-43.824 

-57.824 



mm 

25 

-43.465 

-52.824 



24 

-43.337 

-51.824 

| 22.9 

-43.713 


23 

-43.337 

-51.824 


■fr%*^| 

uBlg j 

22 

-42.991 

-49.824 


-42.506 feiiaBl 

21 

-42.764 

-48.824 

BiffilEl 

-42.375 


20 

-42.494 

-47.824 


-41.532 


1 9 

-42.176 

-46.824 

■ill 

-41.343 

BMiP 

18 

-41.806 

-45.821 


-41.116 


1 7 

-41.374 

-44.807 

MiM 

-40.953 

; . . 

16 

-40.865 

-43.755 


-40.559 


15 

-40.242 

-42.615 

| 15.9 

-40.144 


14.5 

-39.871 

-41.992 


-39.933 ggm 

14 

-39.448 

-41.323 


-39.591 


13.5 

-38.965 

-40.599 

■00 

-39.38 

' . : 

13 

-38.41 

-39.813 

■ESI 

-38.847 


12.5 

-37.776 

-38.959 

13.4 

-38.432 


12 

-37.051 

-38.03 

mz&il 

-37.892 


11.5 

-36.228 

-37.022 

■00! 

-37.122 

foil Is 

1 1 

-35.298 

-35.927 

■mi 

-36.557 

wmm io.5 

-34.254 

-34.74 

■mi 

-36.17 


10 

-33.089 

-33.456 

| 10.9 

-35.095 


9.5 

-31.798 

-32.067 

■ESI 

-34.4s 

1 - 

9 

-30.373 

-30.566 

KH 

-33.1071 

■ 13 - - 1 1 

8.5 

-28.811 

-28.945 

ms 

-32.425 

v.V : !‘ " :.C 

8 

-27.108 

-27.198 

■OJi 

-30.8641 


7.5 

-25.257 

-25.315 

mo 

-29.8121 


7 

-23.255 

-23.292 

IHlS 

-27.931 | 

. i; . - ' . 

6.5 

-21.103 

-21.125 


1 


6 

-18.8 

-18.814 
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Table H-16. The first set of 8PSK hardware variance data at 0.66% 
high-SNR loop bandwidth-to-symbol rate ratio. 



Theoretical 




SNR 

variance with 

Theoretical 

Measured 

(dB) 

added noise floor 

variance 

variance 

1 L 

30 

-40.535 

-54.814 

-40.56 

29 



-40.664 

28 



-40.466 

27 



-40.27 

26 



-40.36 

25 

-40.198 

-49.814 

-40.153 

24.5 



-40.106 

24 

-40.076 

-48.814 

-40.041 

23.5 



-39.795 

23 

-39.929 

-47.814 

-39.597 

22.5 



-39.685 

22 

-39.749 

-46.814 

-39.48 

21.5 



-39.375 

21 

-39.534 

-45.814 

-39.617 

20.5 



-38.955 

20 

-39.277 

-44.814 

-39.309] 

19.5 



-38.694 

1 9 

-38.973 

-43.814 

-38.632 

18.5 



-38.436 


1 8 

-38.618 

-42.811 

-38.104 

17.5 



-38.236 

17 

-38.203 

-41.797 

-37.637 


16.5 



-37.357 


1 6 

-37.712 

-40.744 

-37.306 


15.5 



-37.063 








































_J 


303 




Table H-17. The second set of 8PSK hardware variance data at 
0.66% high-SNR loop bandwidth-to-symbol rate ratio. 



Theoretical 



SNR 

variance with 

Theoretical 

Measured 

m 

added noise floor 

variance 

variance 

15 

-37.107 

-39.604 

-36.568 

14.5 

-36.746 

-38.982 

-36.278 

14 

-36.334 

-38.313 

-35.785 

13.5 

-35.861 

-37.589 

-35.392 

1 3 

-35.318 

-36.803 

-34.756 

12.5 

-34.695 

-35.949 

-34.096 

1 2 

-33.981 

-35.02 

HBESEZS 

11.5 

-33.168 

-34.011 


1 1 


-32.917 

-31.882 

10.5 

-31.212 


-30.813 

10 

-30.054 


-29.454 

9.5 

-28.769 


-28.191 


-27.35 

-27.555 



-25.792 



8 

-24.091 

-24.187 


7.5 

-22.242 

-22.305 


7 

-20.242 

-20.281 


6.5 

-18.091 

-18.115 


6 

-15.789 

-15.803 
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Table H-18. 16PSK hardware variance data at 0.19% high-SNR loop 
bandwidth-to-symbol rate ratio. 


LAB 


SNR for 

Theoretical 


SNR 

Measured 

Theoretical 

variance with 

Theoretical 

(d 

Variance 

Data (dB) £ 

idded noise floor 

variance 


-43.233 


30 

-43.213 

-60.232 


-43.496 


25 

-43.03 

-55.232 

Bjiw 

-43.003 


24 

-42.963 

-54.232 


-43.587 

'inpaiii" 

23 

-42.88 

-53.232 

KL4n 

-43.277 

lirajl' 

22 

-42.777 

-52.232 

EH 

-43.136 


21 

-42.652 

-51.232 


-43.125l^feSr 

20 

-42.402 

-49.687 

ran 

-42.966 


19.5 

-42.25 

-48.933 

mm 

-42.898 

19 

-42.059 

-48.103 

ran 

-42.776 

18.5 

-41.812 

-47.188 

■m 

-42.799 

18 

-41.494 

-46.178 

20 9 

-42.653 llplii 

17.5 

-41.082 

-45.061 

20 4 

-42.535 

17 

-40.545 

-43.827 

IQ 9 

-42.191 

16.51 

-39.851 

-42.462 

IBhfi 

-42.313fc'M%&i 

16 

-38.959 

-40.952 

iBn*i 

-42 125 \ 

15.5 

-37.829 

-39.278 

18 4 


15 

-36 422 

-37.42 

1 7 Q 


14.5 

-34.709 

-35.356 

15m 

-40.591 11%!^! 

14 

-32.669 

-33.062 

IBSB 

-40.048' IS SB 

13.5 

-30.293 

-30.516 



1 13 

-27.578 

-27.696 


-38.354(111;# 

i 12.5 

1 -24.531 

-24.589 

Ira 

-37.427 

1 12 

; -21.16 

i -21.187 

■rag 


i 11.5 

i -17.46 

I -17.491 

Iras 

-33.335 

S 1 1 

-13.506 

( -13.513 

Iras 

> -33 859 

10. £ 

> -9.271£ 

5 -9.2733 

Ibex 

IHZE 

i 
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Table H-19. 16PSK hardware variance data at 0.33% high-SNR loop 
bandwidth-to-symbol rate ratio. 


■sa 



SNR for 

Theoretical 


SNR 

Measured 


■BEH 

variance with 

Theoretical 

m 

Variance 



added noise floor 

variance 

42.9 

-43.466 


40 

-44.231 

-67.824 

37.9 

-43.967 


30 

-44.063 

-57.824 

32.9 

-44.254 


25 

-43.685 

-52.824 

27.9 

-44.104 


24 

-43.55 

-51.824 

22.9 

-43.473 

iifiiii 

23 

-43.386 

-50.824 

22.4 

-43.406 


22 

-43.188 

-49.824 

21.9 

-43.129 


21 

-42.95 

-48.824 

21.4 

-42.939 


20 

-42.512 

-47.328 

20.9 

-43.135 

y : .;; 

19.5 

-42.256 

-46.597 

20.4 

-42.722 

y,' 

19 

-41.945 

-45.798 

19.9 

-42.21 1 


18.5 

-41.563 

-44.923 

19.4 

-42.081 


18 

-41.093 

-43.962 

18.9 

-41.888 


17.5 

-40.517 

-42.908 

18.4 

-41.398; 

• .v 

1 7 

-39.812 

-41.75 

17.9 

-41.029 


16.5 

-38.956 

-40.478 

17.4 

-40.326 


16 

-37.925 

-39.078 

16.9 

-39.703 


15.5 

-36.695 

-37.533 

16.4 

-39.023 


15 

-35.242 

-35.825 

15.9 

-38.276! 

- ■ r: •; 

• •• wgi f f _ 

14.5 

-33.544 

-33.93 

15.4 

-37.185 


14 

-31.579 

-31.82 

14.9 

-36.303 


13.5 

-29.325 

-29.467 

14.4 

-34.703 

■ : ~ 

1 3 

-26.764 

-26.842 

13.9 

-33.224 

: ■ *■' ' ^ 

12.5 

-23.88 

-23.92 

13.4 

-31.945 


12 

-20.665 

-20.684 


| 


11.5 

-17.12 

-17.128 


f 


1 1 

-13.257 

-13.261 


10.5 

-9.1027 

-9.104 
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Table H-20. The first set of 16PSK hardware variance data at 
0.66% high-SNR loop bandwidth-to-symbol rate ratio. 
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Table H-21. The second set of 16PSK hardware variance data at 
0.66% high-SNR loop bandwidth-to-symbol rate ratio. 



Theoretical 



SNR 

variance with 

Theoretical 

Measured 

(dB) 

added noise floor 

variance 

variance 

20 

-39.739 

-44.318 


19.5 

-39.47 

-43.587 

HEMS 

1 9 

-39.143 

-42.788 


18.5 

-38.743 

-41.912 

HEED 

1 8 

-38.254 

-40.952 

HBHH3 

17.5 

-37.656 

-39.898 

heed 

1 7 

-36.928 

-38.74 


16.5 

-36.05 

-37.468 


16 

-34.996 

-36.067 

-35.075 

15.5 

-33.746 

-34.523 

-33.843 

1 5 

-32.275 

-32.81 5 I 

-32.491 

14.5 

-30.563 

-30.919 

-31.217 

14 

-28.587 

-28.81 

-29.658 

13.5 

-26.326 

-26.457 


13 

-23.76 

-23.832 


12.5 

-20.873 

-20.91 


12 

-17.656 

-17.674 
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